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Observation  of  new  type  resonances  in  tripie  barrier  resonant 
tunneling  diodes 

J.  Jo®) 

Department  of  Electronics  Engineering,  Ajou  University,  Suwon,  442-749,  Korea 

K.  Alt  and  K.  L  Wang 

Department  of  Electrical  Engineering,  University  of  California  Los  Angeles,  California  90095 
(Received  21  February  1997;  accepted  for  publication  14  June  1997) 

Current-voltage  characteristics  measured  in  triple  barrier  resonant  tunneling  diodes  reveal  new 
features  about  energy  level  alignment  in  a  superlattice.  Our  data  indicate  that  energy  levels  in  the 
two  quantum  wells  are  not  aligned  at  current  peaks.  Current  peaks  are  observed  when  one  of  the 
energy  levels  in  the  two  wells  becomes  resonant  with  the  emitter  level.  One  of  the  current  peak  was 
thermally  activated,  and  showed  inverted  bistability  at  77  K  temperature.  We  explain  that  this  is  due 
to  X  state  assisted  tunneling.  ©  1997  American  Institute  of  Physics.  [80021-8979(97)04818-4] 


I.  INTRODUCTION 

Since  the  original  proposal  of  a  superlattice  structure  by 
Esaki  and  TsuJ  there  has  been  considerable  interest  in  the 
transport  properties  of  this  structure.  In  a  superlattice, 
current-voltage  characteristics  show  multiple  peaks  due  to 
nonuniform  electric  field  distribution  in  the  structure,  which 
is  called  high  field  domain  formation.^^  It  was  previously 
thought  that  the  electric  field  in  a  superlattice  always  aligns 
the  energy  levels  of  the  neighboring  wells  at  the  same  height. 
However,  Kwok  et  al  recently  reported  that  the  electric  field 
in  the  high  field  domain  is  lower  than  the  field  for  resonant 
alignment."^  Their  photoluminescence  experiments  showed 
that  resonant  alignment  of  energy  levels  occurred  only  for  a 
narrow  region  of  bias.  It  was  explained  that  the  nonresonant 
alignment  is  due  to  current  conservation  across  the  structure. 
Since  the  upper  level  current  is  much  larger  than  the  lower 
level  current,  resonant  alignment  is  not  allowed  because  of 
the  current  difference. 

The  details  of  energy  level  alignment  by  the  applied 
electric  field  is  not  yet  thoroughly  understood.  It  was  not 
easy  to  interpret  superlattice  current-voltage  (I-V)  data 
since  there  are  many  wells  involved,  where  each  well  can 
give  current  peaks  of  a  different  nature.  In  order  to  simplify 
this  problem,  we  studied  the  transport  properties  of  triple 
barrier  resonant  tunneling  diodes.^  Although  these  diodes  are 
a  short  version  of  a  superlattice  structure,  they  still  demon¬ 
strate  the  interactions  between  quantum  wells.  Our  data  re¬ 
veal  new  features  about  energy  level  alignment.  Our  data 
show  that  a  current  peak  is  observed  when  any  energy  level 
in  the  two  wells  is  resonant  with  the  emitter  level,  without 
energy  level  alignment  between  neighboring  wells.  In  addi¬ 
tion,  our  data  show  that  tunneling  between  next-to-nearest 
wells  can  be  as  strong  as  tunneling  between  nearest  neigh¬ 
bors.  One  sample  showed  inverted  bistability,  where  upward 
sweep  shows  lower  current  than  downward  sweep.  We  ex¬ 
plain  that  the  bistability  is  related  to  the  X  state  in  AlAs 
barrier. 
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II.  EXPERIMENTS 

The  structures  used  in  this  study  are  made  of 
Ino.53Gao.47As  wells  and  AlAs  barriers,  grown  by  molecular 
beam  epitaxy  (MBE)  on  InP  substrates.  The  triple  barrier 
structure  is  made  of  (starting  from  the  substrate  side)  AlAs 
barrier,  narrow  InGaAs  well,  AlAs  barrier,  wide  InGaAs 
well,  and  AlAs  barrier.  Three  different  structures  A,  B,  and  C 
were  grown.  In  these  structures,  the  width  of  the  narrow  well 
was  fixed  to  58  A.  The  widths  of  the  wide  well  and  the 
barrier  are  7 1  and  25  A,  respectively,  for  sample  A;  65  and 
25  A  for  sample  B;  65  and  35  A  for  sample  C.  The  contact¬ 
ing  layers  on  the  top  and  the  substrate  side  are  Si  doped 
InGaAs  layers  with  7X  10*^  cm”^  doping  density.  An  un¬ 
doped  InGaAs  spacer  (64  A)  was  inserted  between  the  con¬ 
tacting  layers  and  the  triple  barrier  structure.  An  Au  disk  of 
50  or  100  fjim  diam  was  deposited  for  an  ohmic  contact, 
which  also  serves  as  a  mesa  etch  mask.  Current-voltage 
characteristics  were  measured  at  room  temperature  and  77  K 
using  a  Hewlett-Packard  4142B  dc  source/monitor. 

We  calculated  energy  levels  and  wavefunctions  of  triple 
barrier  resonant  tunneling  diode  (RTD)  structures.  The  en¬ 
ergy  levels  and  wavefunctions  were  obtained  by  matching 
wavefunctions  [exp(±/kx),  exp(±^j:)]  and  their  derivatives 
at  the  six  interfaces  between  InGaAs  and  AlAs.  Nonparabo- 
licity  was  included  by  using  the  method  in  Ref.  6.^’^  Figure  1 
shows  squares  of  electron  wavefunctions  in  sample  B.  The 
calculated  wavefunctions  are  localized  in  a  single  well  with 
only  small  penetration  into  the  other  well,  due  to  the  thick 
barriers.  The  four  energy  levels  in  Fig.  1  are  equivalent  to  the 
first  and  second  energy  levels  in  the  two  wells,  and  we  call 
these  n  i  and  ni  for  the  levels  in  the  narrow  well,  and  and 
W2  for  those  in  the  wide  well.  The  calculated  energies  for 
sample  B  are  Wi  =  l34,  W2  =  389,  ni=159,  and 

^2  =  435  meV.  For  sample  C  (35-A-thick  barriers),  the  cal¬ 
culated  energies  are  within  ±2  meV  of  those  in  sample  B. 

Figures  2  and  3  show  room  temperature  and  77  K  I-V 
characteristics  measured  in  the  triple  barrier  diodes.  At  a 
positive  bias,  the  top  side  is  positively  biased,  and  the  sub¬ 
strate  side  becomes  the  emitter.  Although  the  current  scale  is 
different  for  the  three  samples,  the  peak  voltages  are  nearly 
the  same,  which  are  marked  by  the  vertical  lines.  In  the  fig- 
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RG.  I.  Squares  of  electron  wavefunctions  calculated  in  sample  B.  Substrate 
is  on  the  left  side  of  the  triple  barrier  diode  structure. 


ures,  the  peak  positions  are  named  by  the  energy  levels.  The 
relation  between  the  peaks  and  the  energy  levels  will  be  dis¬ 
cussed  later.  Sample  C  has  thicker  barriers,  so  that  the  peak 
positions  are  slightly  higher  than  those  of  sample  B. 
(Samples  A  and  B  are  different  in  the  wide  well  width,  and 
samples  B  and  C  are  different  in  the  barrier  thickness).  In 
Fig.  2,  two  current  peaks  are  observed  on  the  positive  side 
(0.65  and  1.10  V),  and  two  peaks  on  the  negative  side 
(—0.25  and  —0.82  V).  The  peak  at  —0.82  V  becomes 
stronger  in  the  77  K  data  shown  in  Fig.  3.  In  addition.  Fig.  3 
shows  another  peak  at  0.32  V.  This  peak  becomes  clear  in 
the  dUdV  curve,  and  the  position  is  marked  as  a  dashed  line 
in  Fig.  2.  In  Fig.  3,  sample  B  shows  inverted  bistability 
around  1.2  V. 

In  Fig.  4,  we  show  I-V  data  for  a  larger  bias  range  up  to 
±4  V.  The  upper  curve  is  from  a  50  jim  diam  device  of 
sample  B,  and  the  two  lower  curves  are  from  a  100  /xm  diam 
device  of  sample  C.  The  solid  curves  were  measured  at  room 
temperature,  and  the  dashed  curve  was  measured  at  77  K. 
Note  that  four  peaks  are  observed  on  the  positive  and  nega¬ 
tive  sides  of  sample  C.  Sample  B  (50  /xm  device)  also  shows 
four  peaks  in  the  negative  side.  In  sample  B,  the  small  peak 
positions  (wj  and  wj)  are  the  same  as  those  of  sample  C,  but 
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FIG.  3.  I-V  characteristics  of  samples  B  and  C  (100  fim  diam)  measured  at 
77  K,  and  derivative  of  sample  B’s  I-V.  In  the  I-V  curve  of  sample  B,  the 
arrows  around  1.2  V  indicate  voltage  sweep  direction. 


the  two  large  peak  positions  on  the  negative  side  {rii  and 
are  different  from  those  of  sample  C.  The  ^2 1  ^2  Po¬ 
sitions  of  sample  B  were  not  the  same  when  different  pieces 
from  the  same  wafer  were  measured,  and  we  think  that  this 
variation  is  due  to  peak  voltage  shift  from  the  large  current. 

HI.  DISCUSSIONS 

It  is  important  to  note  that  four  current  peaks  are  ob¬ 
served  in  a  two-well  RTD  structure,  in  which  each  well  has 
two  energy  levels.  Our  structures  have  five  energy  levels: 
one  is  the  emitter  level,  and  the  other  four  are  the  quantum 
well  levels.  If  the  quantum  well  energy  levels  have  to  align 
to  make  current  peaks,  a  maximum  of  three  peaks  should  be 
observed  at  both  polarities  of  bias.  The  three  peaks  at  nega¬ 
tive  bias  would  be  ^  — Wi  — /Z|,  e  — W|  — n2,  and  ^  — W2 
— /i2,  where  e  is  the  emitter  level.  The  alignment  e  — W2 
—  nj  is  not  allowed  because  the  W2  level  is  much  higher  than 
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RG.  2.  I-V  characteristics  of  100  fim  diam  devices  of  samples  A-C,  mea¬ 
sured  at  room  temperature.  The  sequence  of  electron  flow  for  positive  bias  is 
emitter — narrow  well — ^wide  well. 


FIG.  4.  I-V  characteristics  of  50  fim  diam  sample  B  and  100  /tim  diam 
sample  C  for  a  larger  bias  range.  Solid  curves  were  measured  at  room 
temperature,  and  dashed  curve  was  measured  at  77  K. 
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TABLE  L  Energy  levels  obtained  from  the  I-V  peak  positions  of  sample  C, 
and  calculated  from  the  Schrodinger  equation.  Wide  well  width  is  65  A,  and 
narrow  well  width  is  58  A.  X  state  energy  was  obtained  from  Wj  peak  of 
sample  B. 


Energy  level 

Positive  bias 

Negative  bias 

Calculation 

w, 

0.22eV 

0.33eV 

I34meV 

W2 

0.72eV 

0.89eV 

389meV 

ni 

0.25eV 

0.24eV 

159meV 

0.80eV 

0.98eV 

435meV 

X 

0.54eV 

level.  At  a  positive  bias,  the  three  peaks  would  be  e 
—  rii  —  wi,  e  — rti-“W2,  and  — W2.  In  the  level  align¬ 
ment  picture,  a  total  of  six  peaks  should  be  observed,  but  our 
data  show  eight  peaks.  We  believe  that  the  peaks  observed  in 
our  data  are  from  resonant  tunneling  between  the  emitter 
level  and  one  of  the  quantum  well  levels,  not  from  the  level 
alignment  between  neighboring  wells. 

Instead  of  using  a  level  alignment  picture,  we  can  relate 
the  four  current  peaks  to  the  four  individual  energy  levels  in 
Fig.  1.  In  this  picture,  the  current  peaks  show  up  when  any 
quantum  well  level  and  the  emitter  level  are  resonant.  On  the 
positive  side,  the  four  peaks  are  Wi,  nj,  W2,  and  ^2  reso¬ 
nances.  On  the  negative  side,  they  are  nj ,  Wj ,  n2,  and  W2. 
This  is  shown  in  Figs.  2-4.  We  tried  to  calculate  energy 
levels  from  the  peak  positions  in  the  I-V  data.  In  the  calcu¬ 
lation,  we  assumed  that  the  Fermi  level  at  the  end  of  the 
spacer  is  in  resonance  with  the  quantum  well  energy  level. 
Depletion  was  not  included,  since  at  high  doping  density  it 
does  not  make  much  difference  in  peak  voltage.  The  energy 
level  of  a  double  barrier  RTD  (55  A  InGaAs  well,.  25  A  AlAs 
barrier,  2xl0**cm“^  doping  density)  calculated  by  this 
method  was  0.20  eV,  which  is  half  of  the  peak  voltage.  The 
results  for  triple  barrier  diodes  are  listed  in  Table  I.  They  are 
much  larger  than  those  from  the  Schrodinger  equation,  but 
compared  to  0,20  eV  of  the  double  barrier  structure  they  are 
acceptable.  The  values  on  the  negative  bias  are  larger  than 
those  on  the  positive  side,  and  this  could  be  due  to  asymmet¬ 
ric  dopant  distribution.  Although  the  energy  levels  obtained 
this  way  is  larger  than  those  from  the  calculations,  the  over¬ 
all  behavior  supports  that  the  current  peaks  in  our  data  are 
from  the  four  energy  levels. 

The  W2  peak  (l.l  V)  of  sample  B  showed  a  thermally 
activated  behavior.  When  the  temperature  was  lowered,  the 
peak  current  decreased,  while  other  peaks  remained  almost 
unchanged.  This  means  that  around  1,1  V  bias  there  is  an 
energy  level  which  is  slightly  higher  than  W2  level.  But  the 
n2  level  is  too  high  to  be  thermally  activated,  and  we  think 
that  this  is  due  to  the  X  state  in  the  middle  AlAs  barrier.  By 
assuming  that  this  is  X  assisted  W2  peak,  we  obtained  0.54 
eV  for  the  X  state  energy.  The  AlAs  X  state  is  located  at  0. 14 
eV  from  GaAs  conduction  band  edge,  and  the  conduction 
band  energy  difference  between  Ino.53Gao.47As  and  GaAs  is 
estimated  at  about  0.4  eV,  so  that  this  is  in  good  agreement 
with  0.54  eV  obtained  above. 

In  Fig.  3,  sample  B  shows  inverted  bistability  at  77  K  in 
the  region  around  1.2  V.  The  bistability  disappeared  below 
50  K  when  the  peak  current  was  small.  The  inverted  bista¬ 


bility  has  been  reported  before  in  a  double  barrier  structure,^ 
and  it  was  explained  by  charge  build  up  in  the  well.  In  our 
picture,  the  peak  at  1,27  V  is  due  to  the  resonance  between 
the  emitter  level  and  the  second  level  (^2)  in  the  wide  well, 
and  assisted  by  the  X  state  in  the  middle  barrier.  We  think 
that  the  bistability  is  related  to  the  charge  build  up  in  the 
AlAs  barrier.  In  sample  A,  the  wide  well  width  is  71  A,  and 
the  energy  separation  between  W2  level  and  the  X  state  is 
larger,  so  that  the  bistability  is  not  observed. 

One  of  the  strong  peaks  in  the  negative  side  (Fig.  4)  is  a 
resonance  between  the  emitter  and  the  narrow  well  (rt2)*  I^ 
this  bias  polarity,  the  emitter  and  the  narrow  well  are  not  the 
nearest  neighbors,  since  the  wide  well  is  in  between  them.  It 
means  that  a  resonance  between  next-to- nearest  wells  can  be 
as  strong  as  the  resonances  between  the  nearest  wells.  This 
kind  of  resonance  in  a  superlattice  has  not  been  considered 
so  far.  In  a  superlattice,  each  well  can  be  the  emitter  for  the 
wells  in  the  low  energy  side,  and  it  was  assumed  that  only 
the  nearest  neighbor  can  be  a  current  source.  In  contrast,  our 
data  indicate  that  a  well  can  inject  electrons  farther  than  to  its 
nearest  neighbor.  This  can  be  explained  by  considering  the 
length  of  the  localized  electron  wavefiinctions  in  a 
superlattice.^  If  the  wavefiinction  covers  more  than  one  pe¬ 
riod  of  a  superlattice,  it  will  be  possible  to  observe  this  type 
of  resonance. 

There  can  be  other  mechanisms  to  explain  the  four 
peaks,  such  as  phonon  assisted  tunneling.  In  I-V  data,  pho¬ 
non  assisted  tunneling  peak  follows  the  main  peak,  and  is 
very  weak  at  room  temperature.*®  From  the  characteristics  of 
our  peaks,  we  think  that  they  are  not  originated  from  phonon 
assisted  tunneling. 

IV.  CONCLUSION 

We  present  transport  data  measured  in  triple  barrier  reso¬ 
nant  tunneling  diodes.  Current  peaks  are  observed  when  any 
energy  level  in  the  wells  becomes  resonant  with  the  emitter 
level.  Our  data  indicate  that  the  energy  levels  in  the  neigh¬ 
boring  wells  are  not  aligned,  and  the  resonance  condition  can 
also  be  satisfied  between  next-to-nearest  neighboring  wells. 
One  of  the  peak  showed  thermally  activated  behavior  and 
bistability.  We  explain  that  this  is  due  to  X  assisted  tunneling 
when  the  X  state  and  one  of  the  energy  levels  is  energetically 
close. 
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We  studied  capacitance  and  parallel  resistance  in  resonant  tunneling  diodes  as  a  function  of  the 
doping  density  in  the  emitter  and  the  collector  regions.  Capacitance  was  obtained  by  analyzing 
resonance  in  the  admittance  measured.  Our  data  show  that  the  capacitance  varies  with  the  doping 
density,  and  that  the  capacitance  is  smaller  than  the  value  expected  from  the  growth  parameters. 
Electron  density  modulation  exists  around  the  barriers,  and  capacitance  has  doping  density 
dependence  as  a  result  of  the  modulation.  ©  1997  American  Institute  of  Physics. 
[50021-8979(97)05422-4] 


INTRODUCTION 

In  a  resonant  tunneling  diode  (RTD),  capacitance  is  an 
important  parameter  in  determining  its  high  frequency  char¬ 
acteristics.  There  have  been  reports  about  capacitance  mea¬ 
surements  in  RTD  structures, but  the  effect  of  doping  den¬ 
sity  in  the  emitter  and  the  collector  has  not  yet  been 
analyzed.  In  this  work  we  investigated  the  device  parameters 
of  a  RTD  as  a  function  of  the  doping  density.  The  RTD 
capacitance  was  usually  obtained  from  the  slope  of  the 
imaginary  part  of  the  admittance,  which  was  measured  by  an 
impedance  meter  at  a  fixed  frequency.  This  method, 
however,  is  not  accurate  because  a  RTD  device  is  a  combi¬ 
nation  of  resistance,  capacitance,  and  wire  inductance. 

In  order  to  obtain  an  accurate  capacitance  value,  we  used 
resonance  of  the  RTD.  It  was  determined  that  the  Au  wire 
(18  fxm  diameter,  4  mm  long)  used  in  the  RTD  connection 
had  about  8  nH  of  inductance.  This  is  larger  than  the  theo¬ 
retically  expected  value  of  3  nH,^  because  of  the  magnetic 
material  used  in  the  header  where  the  RTD  is  mounted. 
When  a  RTD  is  connected  by  the  Au  wire,  the  whole  circuit 
becomes  a  series  inductance-capacitance  (LC)  circuit  due  to 
the  wire  inductance  and  the  RTD  capacitance,  and  a  LC  reso¬ 
nance  is  observed  in  the  admittance.  For  a  100  fim  diameter 
RTD,  the  resonance  frequency  is  in  the  range  of  300  MHz- 
500  MHz,  depending  on  the  RTD  capacitance.  We  measured 
the  admittance  by  sweeping  the  frequency  up  to  1  GHz,  us¬ 
ing  a  Hewlett-Packard  4 191 A  RF  impedance  analyzer.  The 
RTD  capacitance  and  the  parallel  resistance  were  obtained 
by  fitting  the  calculated  admittance  to  the  measured  data. 
The  admittance  changed  drastically  around  the  LC  resonance 
frequency,  and  the  device  parameters  could  be  easily  ex¬ 
tracted  because  the  resonance  behavior  was  very  sensitive  to 
the  parameter  changes. 


EXPERIMENTS 

The  structures  used  in  this  study  were  made  of 
Ino.53Gao  47AS  layers  grown  by  molecular  beam  epitaxy 
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(MBE)  on  n+  InP  substrates.  Four  samples  (A,  B,  C,  and  D) 
were  investigated.  The  double  barrier  structure  is  composed 
of  an  InGaAs  well  with  AlAs  barriers  attached  on  both  sides 
of  the  well,  all  undoped.  The  growth  parameters  of  the  four 
samples  are  listed  in  Table  I.  In  all  samples  the  well  width 
was  fixed  at  51  A.  The  R,  L,  and  C  values  in  Table  I  were 
obtained  from  the  admittance  data,  the  details  of  which  will 
be  discussed  later.  The  contacting  layers  outside  the  double 
barrier  structure  are  1000  A  thick  Si-doped  InGaAs,  with  the 
doping  density  ranging  from  8.6X10^^  cm""^  to  2.0 
X  10**  cm"^.  An  undoped  InGaAs  spacer  layer  was  inserted 
between  the  contacting  layer  and  the  double  barrier  structure. 
A  100  jjLxn  diameter  Au  disk  was  deposited  as  an  ohmic 
contact,  and  also  served  as  a  mesa  etch  mask.  An  18  /4m 
diameter  Au  wire  of  3-5  mm  length  was  bonded  to  the  Au 
disk  for  the  connection. 

Figure  1  shows  the  current-voltage  characteristics  of 
samples  A  and  C,  measured  at  room  temperature.  The  data 
were  measured  by  a  Hewlett-Packard  4142B  dc  source/ 
monitor.  Since  sample  C  has  a  lower  doping  density,  its  peak 
voltage  is  higher  than  that  of  sample  A.  The  lower  peak 
current  of  sample  A  is  due  to  the  thicker  barriers.  Current 
oscillations  were  observed  in  sample  C  when  it  was  biased  in 
the  negative  differential  resistance  region. 

The  admittance  (T)  of  sample  A  measured  by  the  RF 
impedance  analyzer  is  shown  in  Fig.  2.  Figure  2(a)  is  without 
any  bias  voltage,  and  Fig.  2(b)  is  with  0.45  V  bias  applied. 
The  dots  in  Fig.  2  represent  experimental  data,  and  the  solid 
curves  represent  calculated  results  using  the  circuit  model  of 
the  Fig.  2,  The  real  part  (Fq)  shows  a  sharp  peak  at  the 
resonance  frequency,  while  the  imaginary  part  (Tp^)  changes 
from  a  positive  to  a  negative  value  as  the  frequency  (/)  is 
increased  beyond  the  resonance  frequency.  To  clarify  the 
origin  of  the  inductance  in  our  device,  we  compared  the  in¬ 
ductances  of  the  wire  connected  between  the  terminals  of  the 
header,  one  with  a  RTD  mounted  to  it,  and  one  without  a 
RTD.  The  two  inductances  measured  this  way  were  nearly 
equal  for  the  same  wire  lengths.  In  addition,  the  measured 
inductances  were  proportional  to  the  wire  length.  Thus  it  is 
obvious  that  the  inductance  in  the  RTD  system  originates 
from  the  Au  wire. 
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TABLE  I,  Growth  parameters  and  measured  R,  L,  C  values  for  the  four 
samples. 


A 

B 

C 

D 

AlAs  barrier  (A) 

32 

25 

25 

25 

InGaAs  spacer  (A) 

64 

64 

64 

16 

Doping  density  (10*®  cm) 

2.0 

2.0 

0.86 

46 

^EXT  (O) 

1.8 

1.9 

2.4 

9.9 

^EXT  (nH) 

8.3 

7.0 

8.3 

8.9 

Rd  (O) 

1250 

350 

260 

630 

Cd  (pF) 

16.5 

17.8 

15.8 

30.3 

EQUIVALENT  CIRCUITS  OF  THE  RTD 

We  modeled  the  RTD  as  a  parallel  circuit  oi  Cq  and 
/?£,,  with  an  external  inductor  and  resistor  ^ext) 

connected  in  series  to  the  RTD.  This  model  is  shown  in  Fig. 
2(a).  The  /?,  L,  C  values  listed  in  Table  I  were  obtained  by 
fitting  the  calculated  frequency  characteristics  to  the  /—  Y 
resonance  data,  measured  at  zero  bias  voltage. 

When  a  bias  voltage  is  applied  to  a  RTD,  the  thermionic 
current  increases,  and  the  parallel  resistance  Rq  decreases. 
The  height  of  the  Tq  is  reduced  by  the  applied  bias  as 
shown  in  Fig.  2(b),  since  any  decrease  in  /?^  lowers  the  Q 
factor  of  the  resonance  circuit.  Figure  2(a)  shows  the  admit¬ 
tance  of  sample  A  at  zero  bias  voltage.  In  Fig.  2(a),  the  slope 
oi  Y go  the  low  frequency  region  is  16.5  pF,  which  is  the 
same  as  the  value  obtained  from  the  resonance  fitting.  In  Fig. 
2(b),  the  slope  is  decreased  by  the  applied  bias,  but  the 
slope  is  different  from  the  capacitance  obtained  by  the  fit¬ 
ting.  (The  slope  of  Y go  is  7.2  pF,  while  the  capacitance  from 
the  fitting  is  11.3  pF.)  We  calculated  the  admittance  of  the 


FIG.  1.  Current-voltage  characteristics  of  samples  A  and  C  measured  at 
room  temperature. 


FIG.  2.  Frequency  dependent  admittance  of  sample  A  at  (a)  0  V  and  (b)  0.45 
V  bias  voltage.  The  solid  curves  are  calculated  results,  and  the  dots  are 
measured  data.  The  straight  line  in  (b)  is  for  11.3  pF. 


circuit  in  Fig.  2  for  different  /?£>  and  L^xt  values  to  see  their 
effect  on  the  Ygo  slope. 

Our  calculations  show  that  when  Rq  is  small,  the  slope 
of  K90  is  lower  than  C^.  In  addition,  for  /?£,<60  H  [other 
parameters  being  fixed  to  those  values  in  Fig.  2(b)],  the  pres¬ 
ence  of  Lext  also  lowers  the  slope  of  Ygo-  This  means  that 
the  slope  of  Y go  gives  lower  capacitance  values  when  either 
Lext  or  exists  in  the  circuit.  Since  a  RTD  is  made  of  R 
and  C,  and  is  connected  by  an  Au  wire  which  has  non- 
negligible  inductance,  only  frequency  sweeping  as  shown  in 
Fig.  2  can  give  correct  values. 

In  Fig.  3  we  show  Co  (circles)  and  Ro  (triangles)  mea¬ 
sured  in  sample  A,  as  a  function  of  applied  bias.  The  solid 
line  represents  calculated  capacitance  using  the  depletion 
model  in  Ref.  2,  with  an  initial  separation  of  590  A.  Al¬ 
though  there  is  a  discrepancy  between  the  data  and  the  cal¬ 
culation,  the  measured  and  calculated  capacitances  behave  in 
similar  ways.  Due  to  the  internal  resistance  of  the  impedance 
analyzer,  the  sample  could  not  be  biased  in  the  negative  dif¬ 
ferential  resistance  (NDR)  region  (0.8  V~i.3  V),  where  a 
capacitance  peak  was  observed  before.^  Ro  would  have  large 
negative  and  positive  peaks  in  the  NDR  region,  and  the 
larger  Ro  shown  at  the  higher  bias  is  a  decreasing  part  of  the 
resistance  peak.  In  Fig.  3  Rext  obtained  from  the  fitting  in¬ 
creased  with  applied  bias,  from  1.8  fl  to  4.0  H.  We  think 
that  this  increase  is  due  to  the  series  resistance  in  the  deple¬ 
tion  region  of  the  collector.  At  a  higher  bias  Lext  decreased 
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FIG.  3.  Capacitance  and  parallel  resistance  of  sample  A  as  a  function  of  bias 
voltage.  The  circles  represent  capacitance,  and  the  triangles  represent  paral¬ 
lel  resistance.  The  solid  line  was  calculated  from  the  depletion  model. 


from  8.3  nH  to  7.4  nH.  An  inductance  change  by  applied 
bias  was  recently  reported  in  Ref.  3,  and  it  was  suggested 
that  the  current  delay  time  change  in  the  quantum  well  was 
the  reason  for  the  change.  In  our  structures  the  inductance 
from  the  Au  wire  is  much  larger,  and  at  this  point  we  do  not 
understand  the  reason  for  the  inductance  decrease. 

DISCUSSION 

In  Fig.  4  we  plotted  zero  voltage  capacitance  as  a  func¬ 
tion  of  the  doping  density.  The  data  labelled  U  in  the  figure 
is  from  Ref.  3.  The  inverse  of  the  capacitance  is  shown  by 
circles,  and  the  emitter-collector  distance  from  the  growth 
parameters  is  shown  by  triangles.  Here,  the  inverse  of  the 
capacitance  is  AefC^^  where  A  is  the  area  of  a  RTD,  and  e 
is  the  dielectric  constant  (14  for  InGaAs).  It  was  expected 


FIG.  4.  The  inverse  of  the  measured  capacitance  (circles)  and  the  emitter- 
collector  distance  (triangles)  from  the  growth  parameters.  The  solid  line  is 
the  calculated  distance  between  the  emitter  and  the  collector.  The  data  la¬ 
belled  as  U  are  from  Ref.  3. 


that  the  inverse  of  the  capacitance  at  zero  bias  voltage  would 
directly  reflect  the  emitter-collector  distance,  independent  of 
the  doping  density  Nj^  in  the  emitter.  However,  Fig.  4  shows 
that  Cq  depends  on  ,  and  the  inverse  of  the  capacitance  is 
much  longer  than  the  growth  parameter.  When  yv^=2.0 
X  10*^  cm”^,  the  of  a  100  fMm  diameter  RTD  are  16.5 
pF  (sample  A)  and  15.8  pF  (sample  B),  while  the  of 
4.6X  10*^  cm“^  (sample  D)  is  30.3  pF  under  the  same  con¬ 
ditions.  The  inverse  of  the  capacitance  is  590  A  for  16.5  pF 
and  321  A  for  30.3  pF.  In  comparison,  the  distances  between 
the  emitter  and  the  collector  from  the  growth  parameters  (in¬ 
cluding  spacer  layer)  are  243  A  for  sample  A  and  133  A  for 
sample  D. 

We  discuss  the  doping  dependence  of  the  capacitance  by 
considering  the  shape  of  the  wave  functions  in  the  emitter 
and  the  collector.  The  wave  functions  in  the  emitter  and  in 
the  collector  can  be  expressed  as  (^f(z)exp(z/:xx4-/jtyy),  where 
z  is  along  the  growth  direction.  Without  a  barrier,  ^(z)  will 
be  When  a  barrier  is  introduced  at  z  =  0,  ^(z) 

changes  to  exp(/A:2Z)+R  exp(-//:2^)  for  the  z<0  region, 
where  R  is  the  reflection  coefficient.  If  the  barrier  is  high, 
there  is  full  reflection  (/?=-!,  ^(z)«sin(/:23))»  and  the  elec¬ 
tron  density  has  a  sine  wave  shape  starting  from  zero  at  the 
barrier  interface.  In  the  RTD  case,  the  AlAs  barrier  height  is 
1  eV,  and  the  electron  energy  is  around  100  meV,  so  we  can 
expect  large  reflection  and  electron  density  modulation  near 
the  barriers.  Because  of  the  charge  modulation,  the  distance 
between  the  electron  systems  in  the  emitter  and  the  collector 
is  larger  than  the  growth  parameter.  Since  the  doping  density 
N [)  determines  the  Fermi  energy  level  in  the  doped  region, 
capacitance  will  have  doping  density  dependence.  When  Nq 
increases,  higher  kz  states  are  occupied,  and  electron  density 
will  reach  in  a  shorter  distance. 

When  an  electron  profile  has  a  shape  distributed  over  a 
finite  width,  the  capacitance  should  be  calculated  by  consid¬ 
ering  the  change  in  the  electron  profile  due  to  the  applied 
voltage.^  For  a  RTD  structure  we  can  approximate  the  dis¬ 
tance  between  the  points  where  electron  density  reaches  No 
to  be  the  emitter-collector  distance.  The  capacitance  in  a 
two-dimensional  electron  system  was  calculated  in  a  similar 
way  with  a  numerical  constant  multiplied.^ 

We  plotted  the  distance  calculated  for  sample  A  as  a 
function  of  No  (solid  line)  in  Fig.  4.  The  electron  density 
was  calculated  by  solving  the  Schrodinger  equation  numeri¬ 
cally.  Although  there  is  about  a  100  A  discrepancy,  the  cal¬ 
culated  distance  shows  similar  doping  density  dependence  as 
the  inverse  of  the  capacitance  (circles).  It  is  evident  that  the 
distance  from  the  growth  parameters  (triangles)  is  too  small 
to  explain  the  measured  capacitances.  We  believe  that  the 
general  agreement  between  the  data  and  the  calculation  dem¬ 
onstrates  that  RTD  capacitance  is  determined  by  the  charge 
modulation  in  the  region  around  the  barriers. 

We  reported  previously  that  a  capacitance  peak  was  ob¬ 
served  in  the  NDR  region,^  where  the  peak  was  explained  by 
a  charge  buildup  in  the  quantum  well.  The  charge  buildup  in 
the  well  makes  the  emitter-collector  distance  shorter,  and  the 
capacitance  shows  a  peak  structure.  In  our  previous  work  the 
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adriiittance  was  measured  at  a  single  frequency,  and  the  slope 
was  taken  as  the  capacitance.  As  mentioned  in  the  beginning, 
this  method  gives  lower  capacitance  values,  and  the  actual 
capacitance  value  is  larger  than  the  values  shown  in  Ref.  2. 
Therefore,  our  observation  that  the  measured  capacitance 
was  much  larger  than  the  depletion  capacitance  is  still  valid. 

The  doping  density  dependence  can  be  explained  by  an 
electrostatic  model,  where  electrons  spill  into  the  undoped 
region  in  a  n-i-n  structure.^  However,  in  a  small  dimension 
structure,  the  exact  nature  of  charge  distribution  near  the 
barriers  should  be  described  by  the  Schrodinger  equation, 
and  the  electrostatic  model  can  give  wrong  answers. 

When  N[)  was  increased  to  2X  10^^  cm”^,  Rp  at  zero 
bias  voltage  dropped  to  5  H,  and  the  LC  resonance  behavior 
disappeared  (not  shown  here).  LC  resonance  also  disap¬ 
peared  at  high  bias  when  Rp  became  small.  It  is  believed  that 
the  small  Rp  shorted  the  capacitance  to  change  the  whole 
system  into  a  L-R  circuit. 

In  our  calculations  the  electron  density  in  the  quantum 
well  increases  with  Np .  This  electron  density  is  a  tail  of  the 
wave  functions  in  the  emitter  and  the  collector.  The  electron 
density  in  the  well  increases  with  Np ,  since  higher  kz  is 
occupied  at  higher  Np .  We  can  roughly  estimate  the  order  of 
magnitude  of  the  electron  density  to  be  Np\i{/o/il/p\^.  Here, 
tj/o  ts  ijf  at  the  center  of  the  well,  and  is  the  peak  value  in 
the  doped  region.  In  our  calculations  changed  by  a 

large  amount  when  Np  changed.  The  ratio  at 

Np—ly.  10*^  cm”^  is  about  0.006,  and  the  electron  density 
in  the  well  could  be  as  high  as  10^^  cm""^.  At  this  high 
density  we  can  expect  strong  conduction  between  the  emitter 
and  the  collector.  This  is  in  agreement  with  the  observation 
of  low  Rp  (5  ft)  at  Np-lX  10^^  cm"^. 


CONCLUSIONS 

We  obtained  the  capacitance  and  parallel  resistance  of 
resonant  tunneling  diodes  by  measuring  frequency  dependent 
admittance  around  a  LC  resonance  frequency.  The  measured 
capacitance  showed  doping  density  dependence.  We  ex¬ 
plained  that  electrons  in  the  emitter  and  the  collector  were 
separated  from  the  barriers  by  a  finite  distance,  since  the 
wave  functions  were  a  sine  wave  shape  starting  from  the 
barrier.  Our  calculations  explained  the  behavior  of  Cp  and 
Rp  at  different  Np's.  The  information  obtained  by  this 
method  will  be  important  in  designing  a  RTD  with  improved 
high  frequency  performance. 
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Current- voltage  characteristics  in  InGaAs/AIAs  triple-barrier  resonant  tunneling  diodes  were  investigated.  Our  data  show 
that  a  current  peak  is  observed  when  any  of  the  quantum  well  energy  levels  is  resonant  with  the  emitter  level,  without  level 
alignment  between  neighboring  wells.  When  one  of  the  well  widths  was  increased,  some  peak  voltages  shifted,  while  others 
remained  unchanged.  This  indicates  that  the  peak  positions  are  associated  to  a  specific  well.  When  the  temperature  was 
lowered,  some  peaks  became  stronger,  while  others  became  weaker.  Our  analysis  shows  that  the  weaker  peaks  are  caused  by 
the  resonances  between  the  emitter  and  the  far  side  well.  At  low  temperatures,  current  injection  to  the  far  side  well  is  more 
difficult,  and  the  current  peak  becomes  weaker.  These  results  indicate  the  conculsion  that  neighboring  energy  levels  are  not 
aligned  in  these  structures. 

KEYWORDS:  level  alignment,  resonant  tunneling,  charge  build-up,  temperature  dependent  current,  electric  field  distribution, 
triple-barrier  diode 


1.  Introduction 

Electric  field  distribution  in  a  superlatiice  structure  is  im¬ 
portant  for  detemaining  the  performance  of  vertical  transport 
devices.  At  low-bias  voltage,  the  electric  field  is  uniformly 
distributed  within  the  suructure.  At  higher  voltages,  the  elec¬ 
tric  field  is  no  longer  uniform  because  a  high  field  domain  is 
formed.*^  It  was  previously  thought  that  the  electric  field  in  a 
superlatiice  always  aligns  the  energy  levels  of  the  neighboring 
wells  at  the  same  height.  Recently,  there  was  an  experimen¬ 
tal  report  that  the  electric  field  in  a  superlattice  is  lower  than 
the  field  for  resonant  alignment.-^  The  report  showed  that  res¬ 
onant  alignment  of  energy  levels  occurred  only  for  a  narrow 
region  of  bias.  In  order  to  investigate  electric  field  distribu¬ 
tion,  we  studied  the  current-voltage  (/-V)  characteristics  of 
triple-barrier  resonant  tunneling  diodes.^^  Although  this  struc¬ 
ture  has  only  two  wells,  it  demonstrates  the  interactions  be¬ 
tween  quantum  wells.  Our  data  show  that  a  current  peak  is 
observed  when  any  energy  level  in  the  two  wells  is  resonant 
with  the  emitter  level,  but  shows  no  energy  level  alignment 
between  neighboring  wells.  Peak  voltage  shift  could  be  re¬ 
lated  to  the  well  width  change,  and  current  peaks  showed  two 
different  types  of  temperature  dependence.  These  results  in¬ 
dicate  the  conculsion  that  energy  levels  in  neighboring  wells 
are  not  aligned. 

2.  Experiments 

The  structures  used  in  this  study  are  made  of  Ino.53Gao.47As 
wells  and  A1  As  barriers,  grown  by  MBE  on  n"*"  InP  substrates. 
The  triple  barrier  structure  consists  of  (starting  firom  the  sub¬ 
strate  side)  an  AlAs  barrier,  narrow  InGaAs  well,  AlAs  bar¬ 
rier,  wide  InGaAs  well,  and  AlAs  barrier.  Three  different 
structures;  A,  B  and  C  were  studied.  In  these  structures, 
the  width  of  the  narrow  well  is  fixed  to  58  A.  The  wide-well 
width  is  71  A  for  sample  A,  and  65  A  for  samples  B  and  C. 
The  barrier  thickness  is  25  A  for  A  and  B,  and  35  A  for  C. 
The  contacting  layers  on  the  top  and  the  substrate  side  are  In¬ 
GaAs  layers  with  7  x  lO^^crn*"^  Si  doping  density.  An  Au 
disk  of  lOO'/tm  diameter  was  deposited  for  an  ohmic  con¬ 
tact.  Current-voltage  characteristics  were  measured  using  a 


Hewlett-Packard  41 42B  dc  source- monitor. 

We  calculated  the  energy  levels  and  wavefunctions  by  nu¬ 
merically  solving  the  Schrodinger  equation.  The  calculation 
results  show  two  energy  levels  in  each  well,  and  the  calcu¬ 
lated  energies  for  sample  B  are  134  and  389  meV  fof  the  wide 
well,  and  159  and  435  meV  for  the  narrow  well. 

Figure  1  shows  the  l-V  characteristics  measured  at  room 
temperature.  At  positive  bias  the  top  side  is  positively  biased, 
and  electrons  flow  from  the  substrate  through  the  narrow  and 
wide  wells.  In  this  figure,  the  peak  voltages  are  similar  for 
the  three  samples,  although  the  current  scale  is  very  different. 
At  low  temperature,  one  more  current  peak  was  resolved  at 
0.3  V  (shown  in  Fig.  4  inset),  and  the  position  is  indicated 
by  a  dashed  line  in  Fig.  1.  Sample  C  has  thicker  barriers, 
so  that  the  peak  positions  are  slightly  higher  than  those  of 
sample  B.  Figure  2  shows  the  l-V  data  of  sample  C  taken 
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Fig.  1,  Currcnt-voliagc  characieristics  of  the  three  samples  measured  at 
room  temperature. 
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Voltage  (V) 

fig.  2.  Cuireni -voltage  data  for  a  larger  iHas  range,  measured  at  room  tem¬ 
perature. 

at  higher  bias  range.  Samples  A  and  B  did  not  exhibit  any 
peaks  in  the  higher  voltage  region.  In  this  figure,  four  peaks 
are  observed  in  one  polarity  of  bias  for  sample  C.  If  the  levels 
were  aligned,  the  maximum  number  of  peaks  would  be  three. 
Therefore,  we  can  infer  that  the  four  peaks  in  our  data  are 
from  the  resonances  of  the  four  energy  levels  with  the  emitter 
level,  and  the  peaks  in  the  figure  are  named  by  the  associated 
energy,  levels.  They  are  ni,  tt/i,  n2  and  xv2  on  the  negative 
side,  and  tz;i,  ni,  W2  and  n2  on  the  positive  side. 

Figure  3  shows  how  the  peak  positions  respond  to  well- 
width  change.  In  this  figure,  the  hV  data  of  samples  A  and  B 
are  compared.  Samples  A  and  B  are  different  only  in  wide- 
well  width.  The  width  of  the  wide  well  is  71 A  for  sample 
A  and  65  A  for  sample  B.  In  Fig.  3(a)  the  ni  position  is 
nearly  fixed,  but  the  W2  position  shifted  to  lower  voltage  when 
the  wide-well  width  was  increased  in  sample  A.  Figure  3(b) 
shows  negative-side  peak  positions.  In  this  figure  the  wi  peak 
also  shifted  to  lower  voltage  in  sample  A.  The  Fig.  3(c)  11 K 
/-V  data  clearly  show  that  the  wi  peak  position  is  shifted  to 
lower  voltage.  Figure  3(c)  also  shows  that  the  two  peaks  on 
the  negative  side  (sample  A)  are  merged  into  one  peak  with 
one  wide  negative  differential  resistance  region.  The  shift  of 
the  peaks  in  sample  A  indicates  that  they  are  related  to  the 
wide-well  width  change.  When  the  wide-well  width  was  in¬ 
creased,  the  energy  levels  of  the  wide  well  were  lowered,  and 
the  peak  voltages  were  shifted  to  lower  values.  If  the  neigh¬ 
boring  well  levels  are  aligned,  the  change  of  well  width  in  one 
of  the  wells  should  affect  all  of  the  peak  positions.  Our  data 
show  that  the  level  alignment  does  not  occur  in  this  case. 

In  Figs.  4  and  5  we  show  the  temperature  dependence  of 
/-V  data.  In  Fig.  4,  the  W2  peak  at  1.2  V  became  weaker  at 
low  temperature,  while  the  n\  peak  at  0.7  V  became  stronger. 
On  the  negative  side,  n\  became  weaker,  while  u?i  became 
stronger.  Similar  behavior  is  observed  in  the  I-V  data  of  sam- 


Voltage  (V) 

Fig.  3.  Comparison  of  peak  positions  for  samples  A  and  B.  (a)  is  for  pos¬ 
itive  bias  and  (b)  and  (c)  are  for  negative  bias.  The  wide- well  resonances 
(tui  and  W2)  clearly  show  shift  when  the  width  was  increased  from  65  A 
to  71  A. 

pie  C  shown  in  Fig.  5.  We  believe  that  these  results  also  sup¬ 
port  our  individual  resonance  model.  If  current  peaks  show  up 
as  a  result  of  individual  level  resonance,  the  peaks  can  be  di¬ 
vided  into  two  groups.  One  group  contains  the  resonances  be¬ 
tween  the  emitter  and  the  near  side  well,  and  the  other  group 
contains  the  resonances  between  the  emitter  and  the  far  side 
well.  For  positive  bias,  the  wide  well  is  the  far  side  well  in  our 
structures,  and  in  Figs.  4  and  5  we  can  see  that  the  W2  peak  is 
weaker  at  low  temperature. 

When  the  structure  is  at  the  far  resonance,  the  first  well  can 
be  considered  to  be  a  barrier  to  the  far  side  well.  The  trans¬ 
mission  through  the  first  well  can  be  altered  by  temperature, 
since  the  transmission  depends  on  the  energy  level  position 
and  thermal  activation.  The  enersv  levels  in  the  first  well  are 
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Fig.  4.  Temperaiure'dependence  of  the  resonances  in  sample  B.  m  on  the 
negative  side  and  ti;2  on  the  positive  side  became  weaker  at  low  tempera¬ 
tures.  Inset:  /-V  characteristics  of  sample  B  on  an  expanded  scale,  between 
0  and  0.6  V. 
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Fig.  5.  Temperamre  dependence  in  sample  C.  The  behavior  is  similar  to 
those  in  Fig.  4. 


not  resonant  to  the  emitter,  and  at  room  temperature,  current 
can  flow  via  thermal  excitation.  At  low  temperature,  ther¬ 
mal  activation  becomes  difficult,  and  transmission  through  the 
first  well  decreases.  In  other  words,  the  effective  barrier  be¬ 
comes  thicker  at  low  temperature.  This  explains  why  the  far 
resonance  shows  weaker  resonance  at  low  temperature. 

3.  Conclusion 

We  presented  experimental  evidence  which  supports  our 
individual  resonance  model.  The  number  of  peaks,  well- 
width  dependence  and  temperature  dependence  indicate  that 
the  peaks  in  our  data  are  due  to  resonances  between  the  emit¬ 
ter  and  independent  energy  levels,  without  level  alignment  be¬ 
tween  neighboring  wells. 

We  would  like  to  note  that  this  model  does  not  deny  the  oc¬ 
currence  of  charge  accumulation  in  quantum  wells.  Charge 
accumulation  in  the  wells  could  exist,  but  apparently  the 
charge  accumulation  is  not  strong  enough  to  align  neighbor¬ 
ing  energy  levels  to  the  same  height.  In  some  of  our  low  tem¬ 
perature  I-V  data  we  observed  bistability,  which  we  think  is 
due  to  charge  accumulation. 

We  demonstrated  with  our  data  that  the  resonance  between 
wells  with  one  well  in  between  can  show  peaks,  especially  at 
room  temperature.  So  far.  only  resonances  between  adjoining 
neighbors  were  considered  in  a  superlattice  structure.  How¬ 


ever,  our  data  show  that  a  well  can  inject  electrons  farther  than 
to  its  joining  neighbor.  This  can  be  explained  by  considering 
the  length  of  the  localized  electron  wavefunctions  in  a  super- 
lattice."*^  If  the  wavefunction  covers  more  than  one  period  of 
a  superlattice,  it  will  be  possible  to  observe  this  type  of  reso¬ 
nance. 

In  Fig.  3(c),  two  peaks  were  merged  into  one  wide  peak  at 
77  K  (sample  A).  This  behavior  has  not  been  previously  re¬ 
ported.  In  a  multi-barrier  structure,  it  was  not  clear  whether 
two  closely  spaced  levels  merge  into  one  peak  or  remain  sepa¬ 
rated  by  anti-crossing.  Our  data  show  that  merging  can  occur 
in  such  a  case. 
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Abstract 

We  have  proposed  and  analyzed  a  new  submillimeter  and  deep  millimeter  wave  transit 
time  device  named  BIPOLITT  (Bipolar  Injection  and  Transit  Time).  The  device  structure 
is  similar  to  a  Heterojunction  Bipolar  Transistor  (HBT)  with  the  base  terminal  open 
circuited  to  AC  signal.  The  ambipolar  diffusion  process  in  the  relatively  long  base  is  used 
to  achieve  a  phase  delay  of  the  injection  current  close  to  180°  (relative  to  the  total 
current).  The  collector  depletion  region  is  used  as  the  drift  zone.  As  an  example,  a 
BIPOLITT  diode  for  300-400  GHz  is  designed  and  analyzed.  The  small  signal  specific 
negative  resistance  is  calculated  to  be  about  -3x  10“’Qcm^ .  A  detailed  large  signal 
calculation  at  300  GHz  has  shown  an  efficiency  of  at  least  12%  with  about  2  mW  output 
power  without  spurious  oscillation.  The  BIPOLITT  is  expected  to  work  up  to  at  1  THz. 
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Semiconductor  device  oscillators  in  the  submillimeter  and  deep  millimeter  wave 
frequencies  are  recently  being  pursued  by  using  Resonant  Tunneling  Diode  (RTD)  [1,2]. 
Currently,  RTD  sets  the  record  of  the  highest  oscillation  at  712  GHz  for  semiconductor 
devices.  In  comparison,  the  highest  frequency  transistor  oscillator  to  date  is  213  GHz.  [3] 
However,  the  power  of  RTD  is  very  low,  usually  in  the  jJW  range  [1,2].  The  recent  result 
of  a  16  element  RTD  oscillator  arrays  designed  with  Schottky  collectors  and  special  bias 
circuits  has  reported  28  yW  at  290  GHz  [2].  The  spurious  oscillation  [4,5]  due  to  the 
inherent  dc  negative  resistance  of  the  RTD  is  the  major  difficulty  in  effectively  applying 
the  device,  despite  clever  remedies  proposed  [5]. 

In  this  letter,  we  report  for  the  first  time  the  simulation  results  of  a  proposed  new 
transit  time  device  [6,7]  called  BIPOLITT  (Bipolar  Injection  and  Transit  Time)  for 
applications  in  submillimeter  and  deep  millimeter  frequencies.  It  uses  a  bipolar  injection 
into  the  drift  zone  in  a  structure  like  HBT  to  generate  RF  powers.  Detailed  simulations  in 
submillimeter  and  deep  millimeter  wave  frequencies  show  good  negative  resistance  and 
decent  efficiency  comparable  to  RTD  but  without  the  problem  of  spurious  oscillation. 

The  transit  time  effects  of  the  bipolar  junction  transistor  (BIT)  have  been  known 
in  the  early  development  of  BJTs  [8].  The  effect  has  recently  been  analyzed  in  much  more 
detail  for  applications  using  modem  HBTs  [9].  Even  though  negative  resistance  is 
expected,  its  has  not  been  clearly  demonstrated  in  experiment.  One  severe  difficulty  of 
applying  the  transit  time  effects  in  transistors  at  high  frequencies  is  the  detrimental  RF 
signal  loss  which  must  come  in  through  the  base  terminal.  The  other  thing  hindering  the 
effective  use  of  the  transit  time  effects  in  transistors  so  far  is  the  limited  optimization  of 
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the  effects,  i.e.,  only  collector  transit  time  is  optimized  [8,9].  We  approach  the  problem 
from  the  point  of  view  of  a  two  terminal  device  and  show  for  the  first  time:  (1)  the 
complete  circumvention  of  the  RF  signal  loss  through  the  base  by  letting  the  base 
terminal  RF  open;  (2)  the  importance  of  optimizing  the  base  transit  time  to  drastically 
enhance  performance. 

The  basic  structure  of  the  proposed  BIPOLITT  is  illustrated  in  Fig.  1(a). 
Electrons  are  injected  from  the  emitter-base  junction  into  the  base.  The  neutral  base 
region  will  provide  a  phase  delay  by  the  diffusion  process.  The  delayed  electrons  are  then 
injected  into  the  collector  depletion  region  which  acts  as  the  drift  zone  for  producing  the 
transit  time  negative  resistance.  The  HBT  technology  allows  ftie  base  doping  to  become 
higher  than  10'’cw"^[10]  with  a  conductivity  therefore  larger  than  10^Q”'caw"'[7].  With 
this  high  conductivity,  the  dielectric  relaxation  time  [7]  is  less  than  10  femtosecond  and 
the  displacement  current  is  less  than  10%  of  the  conduction  current  up  to  1  THz.  This 
means  the  charge  neutrality  implicitly  assumed  in  the  ambipolar  diffusion  equation  [7,1 1] 
in  the  base  is  valid  up  to  1  THz. 

In  order  to  fully  utilize  the  transit  time,  the  injected  current  needs  to  have  a  phase 
delay  (relative  to  the  total  current)  close  to  180°  but  still  with  a  significant  magnitude 
[6,7].  Is  this  feasible  by  the  difftision  process  in  base?  A  quick  way  to  the  answer  can  be 

obtained  by  investigating  the  complete  AC  base  transport  factor  which  is  the  ratio  of 
electron  diffusion  current  at  the  end  of  the  neutral  base  region  to  that  at  the  beginning.  It 
is  well  known  that  [11] 
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a  j*  =  sec  =  sqc  h{j2A'io)  I  a  =  sec  /i(r  )  ,  (1) 

^nb 

where  is  the  thickness  of  the  neutral  base,  I.„j  =  is  the  dc  diffusion  length, 

(o  is  the  angular  frequency,  co^  =  2A3D„^  /  is  the  base-cutoff  angular  frequency,  is 
the  base  minority  carrier  life  time,  and  we  define  T  =  ^  j2A3a>  /  co^  .  Since  we  are 
interested  in  the  deep  millimeter  and  submillimeter  wave  frequencies,  » 1  and  the 
term  1  in  Eq.(l)  can  be  neglected.  At  the  base-cutoff  n)*,  it  is  known  that  the  magnitude 

of  the  aY  is  reduced  to  about  71%  with  a  corresponding  phase  delay  of  58°  [11].  To 
obtain  a  phase  delay  of  180°  for  a^-*,  it  simply  requires  T  =  n+ jn .  The  requirement 
demands  0  =  8.66)4  and  | <3:^*1  =0.09  .  Then  the  collector  injection  current  have  a  phase 
delay  of  180°  and  a  magnitude  of  9  %  relative  to  the  conduction  current  at  the  emitter 
base  junction.  The  base  thickness  should  be  designed  to  have  a  base  cutoff  frequency  of 
12%  of  the  operating  frequency.  Therefore,  it  is  possible  to  achieve  a  good  injection 
phase  delay  to  generate  RF  power  with  an  efficiency  of  order  of  10  %  ,  which  is  very 
good  at  deep  millimeter  and  submillimeter  wave  frequencies. 

The  accurate  account  of  injection  current  phase  delay  should  refer  to  the  total 
current  of  the  device.  The  ratio  y  of  the  injection  current  with  respect  to  the  total  current 

can  be  shown  as 

(l/rJTcscKr) 

(l/);)rcoth(r)  +  76)Q, 
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where  the  emitter  junction  dc  resistance  =  kT  I  ql^ ,  with  as  the  dc  emitter  current. 
The  is  the  emitter  junction  depletion  capacitance.  The  emitter  efficiency  is  counted  as 
one  in  Eq.(2)  since  a  heterojunction  at  the  emitter  junction  is  used. 

The  numerator  in  Eq.(2)  is  the  trans-admittance  of  the  conduction  current  at  the 
beginning  of  the  collector  drift  region.  We  call  it  a  trans-admittance  because  the  voltage  is 
referenced  to  the  emitter-base  junction.  The  denominator  in  Eq.(2)  is  the  total  admittance 
of  the  emitter-base  jimction.  The  equivalent  circuit  for  the  emitter  junction  can  be 
constructed  accordingly.  Note  that  Eq.(2)  is  reduced  to  Eq.(l)  if  the  displacement 
admittance  y<yQ^is  negligible  when  compared  with  the  admittance  (l/?;)rcoth(r)of 
the  emitter  junction  conduction  current.  This  is  usually  the  case.  Therefore,  the  earlier 
discussion  of  the  co  dependence  of  aj.*provides  a  general  and  simple  trend  for  y  as  well. 


The  equivalent  circuit  elements  for  the  emitter-base  junction  have  therefore  three 
parallel  elements:  depletion  capacitance,  AC  jimction  conductance  and  AC  diffusion 
capacitance.  The  exact  equivalent  circuit  in  the  base  region  is  very  complicated  due  to  the 
ambipolar  diffusion  process.  However,  since  the  base  doping  is  sufficiently  high,  it  is  a 
very  good  approximation  to  represent  it  by  the  base  resistance  parallel  with  the  trans¬ 
admittance. 

The  equivalent  circuit  element  for  the  collector  depletion  region  is  the  well  known 
transit  time  impedance  of  the  drift  region.  It  is  determined  [6,7]  by  the  injection  ratio 
y  and  transit  angle  6^  =  o)W^  /  , 


7^. 


(3) 
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where  C^=sJWj^\%  the  specific  capacitance  of  the  drift  region,  with and the 
corresponding  length  and  the  saturation  velocity  ftierein.  Both  diffusion  and  velocity 
overshoot  are  neglected  in  Eq.(3).  We  have  carried  out  more  detailed  sunulations 
including  both  effects.  Since  Eq.(3)  is  obtained  by  an  integration  of  the  carriers  in  the 
drift  region  [6,7],  we  find  the  neglect  of  the  diffusion  is  not  serious  for  negative 
resistance.  But  we  find  as  expected  the  diffusion  effects  are  important  for  the  accurate 
calculation  of  efficiency.  The  velocity  overshoot  effects  is  found  to  be  negligible  as  long 
as  the  electric  field  in  the  drift  region  is  sufficiently  high. 

The  resultant  small  signal  equivalent  circuit  including  both  the  injector  and  the 
drift  zone  of  the  BEPOLITT  is  shown  in  Fig.l(b).  We  have  accordingly  designed  a 
BIPOLITT  for  300  to  400  GHz  operation  and  have  carried  out  extensive  small-signal  and 
large  signal  simulations.  The  details  will  be  published  elsewhere.  In  the  following,  we 
summarize  our  design  and  simulation  results.  The  thickness  of  the  base  region  (  GaAs 
doped  to  10‘^cw"^)  is  designed  to  be  0.1 5  5//m,  corresponding  to  a  base-cutoff  fi-equency 
of  65  GHz.  The  injection  phase  delay  is  about  170°  at  400GHz  and  about  150°  at  300 
GHz.  They  are  slightly  larger  than  180°  to  take  advantage  of  a  larger  magnitude  for  y  ■ 

The  thickness  of  the  collector  region  (doped  to  W' cnf^)  is  0.14//W  corresponding  to  a 
transit  angle  of  about  180°  at  400  GHz.  The  calculated  specific  resistance  of  the  device 
shows  a  very  broad  negative  resistance  fi:om  250  to  450  GHz.  Compared  to  RTD,  this 
negative  resistance  band  is  sufficient  to  circumvent  the  spurious  oscillation  problem.  At 
300  GHz  with  a  bias  current  of  W  A!  crn^ ,  the  specific  resistance  of  the  collector  is 
-  5  X 10”’  Q  •  cm^ .  The  specific  resistance  of  the  base  is  1.3  x  10”’  Q  •  cm^ .  The  specific 
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contact  resistance  can  now  be  made  smaller  dian  10  ’Q-c/w^[10].  With  an  emitter 
(AlGaAs)  doped  to  lO'^c/n"^  with  a  thickness  of  0.05/^,  the  total  specific  resistance  of 

the  device  is  -  3  x  10"’  Q  •  cnP'  at  300  GHz ,  comparable  to  that  of  RTD  [13]. 

From  the  large  signal  simulation,  we  find  the  efficiency  is  at  least  12%  at  300 
GHz.  This  efficiency  is  also  comparable  to  that  of  RTD  and  much  better  than  IMP  ATT. 
The  corresponding  efficiency  of  the  HBT  structure  designed  in  [9]  is  found  to  be  3%. 
(This  low  efficiency  is  due  to  a  phase  delay  too  small  in  the  base,  less  than  40° .)  Our 
device  can  be  biased  by  a  crurent  source  through  the  emitter.  When  matched  to  a  load  of  3 
Q,  an  emitter  area  of  about  10  /jn?  can  be  used.  A  resistor  of  about  10  ^  across  the 
base  collector  terminal  will  provide  a  DC  bias  voltage  of  3  volt  if  the  common  emitter 
current  gain  p  is  about  40.  This  bias  resistor  is  effectively  an  open  circuit  for  the  AC 

signal.  Therefore,  the  foregoing  two  terminal  analysis  is  valid  and  accurate.  The 
maximum  total  ou^ut  power  of  a  single  device  is  about  2mW  at  300  GHz,  limited  by  the 
breakdown  voltage  and  the  carrier  space  charge  effects  in  the  collector.  Further 
calculations  show  that  BIPOLITT  can  work  up  to  1  THz  with  an  efficiency  of  about  5%. 
We  have  not  optimized  the  performance  by  choosing  other  material  systems  and  have  not 
explored  the  effects  of  a  drift  field  in  the  base  region.  The  BIPOLITT  can  also  be 
integrated  in  series  without  the  dc  instability  problem  of  RTD  [14]. 

This  work  was  supported  by  the  UCLA  JSEP  program. 
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CAPTIONS 


Fig.  1.  (a)  A  schematic  device  structure  for  BIPOLITT  diode,  (b)  The  corresponding 
small  signal  equivalent  circuit  for  the  BIPOLITT  diode.  Note  that  AC  junction  resistance 
and  AC  diffusion  capacitance  are  used. 

Fig.  2.  The  calculated  specific  AC  negative  resistance  of  a  BIPOLITT  diode  designed 
for  300  ~  400  GHz  at  three  bias  current  densities  .  The  key  parameters  of  the  device 


structure  are  also  shown. 


(b) 


Figure  1 
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Small  Signal  Analysis  of  Bipolar  Injection 
Transit  Time  Diode  (BIPOLITT) 


L.  Chen  and  D.  S.  Pan 


Abstract — We  have  proposed  and  analyzed  a  new  millimeter 
and  submillimeter  wave  transit  time  device  named  BIPOLITT 
(Bipolar  Injection  and  Transit  Time).  The  device  structure  is 
similar  to  a  Heterojunction  Bipolar  Transistor  (HBT)  with  a 
floating  base.  The  relatively  longer  base  is  utilized  to  achieve 
an  injection-phase  delay  close  to  180°.  The  collector  depletion 
region  is  used  as  drift  region.  As  an  example,  a  BIPOLITT 
diode  for  300-400  GHz  operation  is  designed,  and  its  small 
signal  specific  negative  resistance  is  calculated  to  be  about 

-3xl0’’Q-cm^ 

Index  Terms — BIPOLITT,  HBT,  Transit  Time 
I.  Introduction 

The  negative  resistance  arising  from  transit  time  in  the 
semiconductor  diodes  was  first  considered  by  Shockely 
[1,2]  in  1954.  It  was  shown  that  when  a  conduction  current 
pulse  is  injected  into  a  drift  region  with  a  phase  delay  with 
respect  to  the  total  current,  we  would  have  a  negative 
resistance  at  certain  transit  angles.  Due  to  the  rapid 
development  of  FETs  and  HBTs,  the  tamsit  time  devices 
such  as  IMPATT  [2,3]  are  now  used  in  deep  millimeter 
wave  frequencies.  In  submillimeter  wave  frequencies,  there 
is  currently  no  good  solid  state  source  except  perhaps  the 
Resonant  Tunneling  Diode  (RTD)  [3].  The  performance  of 
IMPATT  is  not  satisfactory  in  deep  millimeter  wave 
frequencies  and  above  (with  efficiencies  lower  than  a  few 
percent  [3]).  This  low  efficiency  is  primarily  caused  by  the 
saturation  of  the  ionization  coefficients  when  the  electric 
field  is  very  high  [2,3].  As  a  result,  the  impact  ionization 
injection  cannot  be  well  localized  and  renders  the  IMPATT 
as  a  low  efficiency  Misawa  diode  [2,3]. 

In  this  paper,  we  investigated  the  possibility  to  employ  a 
HBT  structure  with  a  relatively  longer  floating  base  to 
realize  the  current  injection.  The  desired  phase  delay  is 
introduced  when  the  electrons  move  across  the  so-called 
base  (injection  region)  in  our  designed  diode,  which  is 
similar  to  the  delay  introduced  when  electrons  flow  through 
the  base  in  the  Bipolar  transistor.  That  is  the  reason  why  we 
call  our  diode  as  the  Bipolar  Injection  Transit  Time  diode 
(BffOLITT).  Since  the  injection  mechanism  is  well 
localized  at  the  beginning  of  the  drift  region,  the 
performance  of  the  BIPOLITT  diode  can  be  made  much 
better  than  IMPATT  at  deep  millimeter  wave  frequencies 


The  authors  are  with  the  Department  of  Electrical  Engineering,  University 
of  California,  Los  Angeles,  CA  90024.  This  work  is  supported  by  the 
UCVLA  JSEP  program 


and  above.  Meanwhile,  with  the  mature  technology  for 
HBTs,  it  would  be  easy  to  fabricate  the  BIPOLITT  diodes. 

Actually,  the  transit  time  effects  have  been  long  known 
for  the  Bipolar  transistors.  In  1966,  J.  F.  Gibbons  [4] 
pointed  out  that  there  may  exist  a  transit  time  mode 
oscillation  in  BJTs  because  of  the  excess  phase  delay  due  to 
the  transit  time  effect.  In  1988,  Sheila  Prasad  etc.  [5,6]  have 
analyzed  the  High-Frequency  characteristics  of  HBTs  and 
pointed  out  that  the  transit-time  effect  plays  a  significant 
role  at  high  frequencies.  Their  numerical  calculation  shows 
that  Re(/t2j)  assumes  negative  values  in  some  portions  of 
the  frequency,  and  the  unilateral  gain  would  peak  when 
Re(/jj2)  =  0.  However,  those  results  are  not  observed  in 
the  commercially  designed  BJTs  or  HBTs.  Based  on  our 
analysis,  we  found  out  that  the  main  reason  is  because  in 
those  BJTs  or  HBTs,  the  transit  time  effect  was  not 
optimized  and  not  sufficient  to  overcome  the  parasitic  to 
obtain  an  overall  negative  resistance.  The  detail  will  be 
discussed  later. 
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Fig.  1  A  schematic  device  structure  for  the  BIPOLITT 
diode 


IL  DEVICE  STRUCTURE 

The  basic  structure  of  the  proposed  BIPOLITT  diode  is 
illustrated  in  Fig.  1.  Basically,  it  looks  like  a  HBT  with  the 
base  floating.  Compared  to  the  regular  HBTs,  the 
BIPOLITT  has  a  relatively  longer  so-called  base  and  a 
shorter  collector.  Similar  to  all  transit  time  devices  [2],  a 
BIPOLITT  diode  consists  of  an  injection  region  and  an 
adjacent  drift  region.  Electrons  are  injected  from  the 
emitter-base  junction  into  the  base  region.  The 
AlGaAs/GaAs  Heterostructure  here  is  utilized  to  improve 


the  current  injection  efficiency  by  reducing  the  hole  current 
significantly  due  to  the  presence  of  the  valence  band 
discontinuity.  The  neutral  base  region  will  provide  a  phase 
delay  by  the  diffiision  process.  The  delayed  electrons  are 
then  injected  into  the  collector  depletion  region  which  acts 
as  the  drift  zone  for  producing  the  transit  time  negative 
resistance.  The  HBT  technology  allows  the  base  doping  to 
become  higher  than  [7]  with  conductivity 

therefore  larger  than  00Q”^*cm  \  With  this  high 
conductivity,  the  dielectric  relaxation  time  [3]  is  less  than 
10  fs,  and  &e  displacement  current  is  less  than  10%  of  the 
conduction  current  up  to  1  THz.  This  means  the  charge 
neutrality  implicitly  assumed  in  the  ambipolar  diffusion 
equation  [3]  in  the  base  is  valid  up  to  1  THz. 


III.  SMALL  SIGNAL  ANALYSIS 

It  is  well  known  that  the  ac  impedance  of  the  drift  region 
is  given  as  [2,3] 


1  r€-e“’‘h 

— [1--^ - ] 

jcoCj  je 


where  Cd=  Es/Wd  is  the  capacitance  per  unit  area  of  the  drift 
region,  and  the  transit  angle  0d  equals  coWd  /Vg,  with  Wd  the 
drift  region  length  and  Vg  the  carrier  saturation  velocity  in 
the  drift  region. 

Both  the  diffusion  and  velocity  overshoot  effects  are 
neglected  in  Eq.  (1).  Since  Eq.  (1)  is  obtained  by  an 
integration  of  the  carriers  in  the  drift  region  [2,3],  the 
neglect  of  the  diffusion  will  not  be  too  serious.  The 
inclusion  of  the  velocity  overshoot  is  expected  to  improve 
the  negative  conductance. 

The  term  y  is  defined  as  the  ratio  of  the  injection  current 
at  the  beginning  of  the  drift  region  to  the  total  ac  current 

and  can  be  written  as  /  =  |  ,  where  ^  is  the  so  called 

injection  phase  delay.  For  non-zero  (|),  the  real  part  of  Zd 
will  show  negative  resistance  and  reach  the  optimum  value 
when  the  sum  of  ^  and  0d  is  equal  to  2n. 

The  calculation  of  y  depends  on  the  current  injection 
mechanism.  As  mentioned  above,  we  employ  the  HBT 
structure  to  realize  the  current  injection  in  the  BIPOLITT 
diode.  Assuming  a  uniform  doping  and  by  solving  the 
diffusion  equation  of  electrons  in  the  injection  region 
(neutral  base  region)  at  high  fi'equency  [8],  we  can  obtain 


^  (i/r.)rcsch(r) 

(i/r,)rcoth(r)+;<«:», 

where  the  emitter  junction  dc  resistance  =  kTlql^ ,  with 

as  the  dc  emitter  current.  The  is  the  emitter  junction 
capacitance  and  o  is  the  angular  firequency.  The  emitter 


efficiency  is  counted  as  1  in  Eq.  (2)  since  a  heterojunction 
at  the  emitter  junction  is  used  and  the  hole  current  is  highly 
compressed.  The  term  T  is  defined  as 


W 

r  =  -f^(l  +j(orJ'^  =  U2A3col0),r  0) 

^nb 

where  Wb  is  the  thickness  of  the  neutral  base, 
^nb  ”  yl^ub'^nb  diffusion  length;  Xnb  is  the  base 

electron  life  time,  and  o)^  =2A?>  is  the  base- 

cutoff  angular  frequency.  Since  we  are  interested  in  the 
deep  millimeter  and  submillimeter  wave  fi-equencies, 
(DXnb>l  and  the  term  1  in  Eq.  (3)  can  therefore  be  neglected. 


Fig.  2  Frequency  variation  of  p* 


Note  that  if  the  displacement  admittance  jcoC^^  is 
negligible  when  compared  with  the  admittance 
(l/r^)r  coth(r),  Eq.  (2)  is  reduced  to  the  base  transport 

factor  =  SQch(J2A3  coj 0) at  high  fi*equencies.  p* 
is  the  ratio  of  the  electron  diffusion  current  at  the  end  to  that 
at  the  beginning  of  the  neutral  base  region,  and  the 
fi'equency  response  characteristic  of  P  is  shown  in  Fig.  2. 
As  we  can  see,  at  the  base-cutoff  fi'equency  the 

magnitude  of  the  P*  is  reduced  to  about  71%  with  a 
corresponding  phase  delay  of  58°.  However,  this  phase 
delay  is  not  sufficient  for  a  good  transit  time  device.  More 

phase  delay  is  required.  To  obtain  a  80°  phase  delay  of  P*, 
it  demands  0)  =  S.6C0fj  and  makes  |/?*|=0,09.  This 

means  that  the  collector  injection  current  can  have  a  phase 
delay  of  80°  with  a  magnitude  of  9%  with  reference  to  the 
current  at  the  emitter  base  junction,  when  the  base  thickness 
has  a  base  cutoff  firequency  of  12%  of  the  operating 
frequency. 

The  numerator  in  Eq.(2)  is  the  trans-admittance  of  the 
conduction  current  at  the  end  of  the  neutral  base  region.  We 


call  it  a  trans-admittance  because  the  voltage  is  referenced 
to  the  emitter-base  junction.  The  denominator  in  Eq.(2)  is 
the  total  admittance  of  the  emitter-base  junction.  The 
equivalent  circuit  for  the  emitter  junction  can  be 
constructed  accordingly.  The  effect  of  Ac  diJEfusion 
capacitance  is  included  in  (l/r^)r  coth(r)  since  the  Ac 

difiusion  length  is  used  in  T  [8].  It  turns  out  that  the  is 
usually  small  when  compared  with  the  Ac  diffusion 
capacitance  when  the  DC  biased  current  is  high  enough, 
and  only  makes  the  phase  delay  slightly  larger  and  the 
magnitude  of  slightly  smaller.  Therefore,  the  earlier 

discussion  of  the  0)  dependence  of  j3*  provides  a  general 
and  simple  trend  for  /  as  well.  When  the  injection  phase 

delay  is  about  80°,  we  expect  an  injection  efficiency  yof 
about  10%,  which  is  sufficient  for  a  transit  time  device. 

The  equivalent  circuit  elements  for  the  emitter-base 
junction  compose  three  parallel  elements:  depletion 
capacitance  Ac  jxmetion  resistance  ^be  and  Ac  diffiision 
capacitance  C^.  The  exact  equivalent  circuit  in  the  base 
region  is  very  complicated  due  to  the  ambipolar  diffusion 
process.  However,  since  the  base  doping  is  sufficiently 
high,  it  is  a  very  good  approximation  Represent  it  by  the 
base  resistance  parallel  with  a  voltage-controlled  current 

source  where  is  the  Ac  voltage  across  the 

emitter-base  junction,  and  a  is  the  trans-admittance  as 
shown  in  the  numerator  in  Eq.  (2). 


Fig.  3  The  small  signal  equivalent  circuit  for  BIPOLITT 

The  resultant  small  signal  equivalent  circuit  including 
both  injector  and  the  drift  zone  of  the  BIPOLITT  is  shown 
in  Fig.3.  is  the  impedance  of  the  drift  region,  and  is 
the  impedance  for  the  Ohmic  contact  of  the  collector.  Z^  is 
the  sum  of  the  Ohmic  contact  impedance  of  the  emitter  and 
the  extrinsic  emitter  series  resistance.  It  should  be  pointed 
out  that  even  though  the  physics  of  fbe  and  is  clear,  we 
do  not  give  the  simplified  formula  for  them,  like  in  the 
regular  analysis  in  BJTs  or  HBTs  [4-6,8].  In  fact,  both 
and  are  dependent  of  fi’equencies.  However,  since  both 
and  are  much  smaller  than  the  drift  region 


impedance  Z^,  the  replacement  of  and  by  their  BIT 
first  order  simplified  formulae  will  not  cause  serious  errors 
in  the  calculation  of  the  total  device  impedance. 

IV.  NUMERICAL  RESULTS 

According  to  the  device  operating  principle  and  the 
equivalent  circuit,  we  have  designed  a  BIPOLITT  for  300 
to  400  GHz  operation.  The  parameters  for  the  device  is 
shown  in  Fig.  4.  The  thickness  of  the  base  region  is 
designed  to  be  OlSSjum,  corresponding  to  a  base-cutoff 
frequency  of  about  65  GHz.  The  injection  phase  delay  is 
about  70°  at  400GHz  and  about  50°  at  300  GHz.  They 
are  slightly  larger  than  80°  to  take  advantage  of  a  larger 
magnitude  for  y  .  The  thickness  of  the  collector  region  is 

014jLm  corresponding  to  a  transit  angle  of  about  80  at 
400  GHz.  Without  including  the  ohmic  contact  resistance, 
the  calculated  specific  resistances  of  the  device  are  depicted 
fi*om  200  to  800  GHz  for  the  bias  current  density  of 

0^  A/cm^  as  shown  in  Fig.  2.  The  solid  line  represents 
the  total  Ac  resistance  of  the  device,  while  the  dotted  line 
represents  the  resistance  of  the  drift  region  only.  In  the 
designed  fi‘equency  range,  they  are  approximately 

-3  X  10”^  Q  •  cm^  and  x  10"^  Q  •  cm^,  respectively. 
These  values  are  good  and  are  comparable  to  that  of  the 
RTD  [9]. 
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Fig.  4  The  calculated  spedfic  Ac  negative  resistance 
of  a  BIPOLITT  diode 

Note  that  the  specific  resistance  of  the  base  of 
0.1 5 5//m  thickness  with  doping  is  about 

.3  X  10"^Q  •  cm^,  which  decreases  the  total  negative 
resistance  by  the  same  value.  It  is  therefore  desirable  to 

dope  the  base  more  heavily  than  O^^cm’^,  for  example 
5  X  which  is  possible  in  today’s  technology,  to 

decrease  the  base  resistance  to  be  lower  than 

4  X  10"*Ocm^. 


BIPOUTTDiode: 

Emrtter(AIGaAs):  eloping  10^®/cm^;  length  0.05  pm 
Base  (GaAs,  floating):  doping  10^®/cm^;  length  0.155  pnj 
i  Collector  ( GaAs ):  doping  10^^/cm^;  length  0.14  \sm 


400  GHz) 


R  in  drift  region 
■  Rtotal 


At  the  same  time,  we  must  point  out  that  the  DC  biased 
current  plays  a  critical  role  in  the  BIPOLITT 
performance.  affects  the  term  y  through  by  changing 
the  weight  of  the  displacement  current  in  the  total  Ac 
current.  When  is  high  enough,  =  kTjqlg  becomes 
very  small  so  that  the  depletion  capacitance  Cbe  can  be 
ignored  when  compared  to  the  Ac  diffusion  capacitance 
since  is  proportional  to  /a;  at  the  first  order 

approximation.  Since  the  is  about  —5  x  10"^  Q  •  cm^ 
and  the  total  negative  resistance  is  decreased  by  the  emitter- 
base  junction  resistance  he  which  is  equal  to  at  the  first 
order  approximation,  therefore  in  order  to  have  a  relative 
large  total  specific  negative  resistance,  7^  should  be  larger 

than  6  X 10"^  v4/cm^,  which  gives 

=  4.4  X  10"^Q  •  cm^.  Taking  the  parasitic  Ohmic 
contact  resistance  into  consideration,  it  is  desired  to  have 
larger  than  8x10"^  AjovcP' ^  which  gives 

“  3  X  cm^.  Meanwhile,  a  0^  Ajzvci  Dc 

biased  current  will  give  =  2.6  x  10~^Q  *  cm^  and  we 
can  safely  ignore  the  he  compared  with  Z^.  However, 
when  the  Dc  current  is  high,  the  high  current  effects 
happened  in  Bipolar  transistor  will  occur,  which  means  that 
a  good  DC  biased  current  density  range  in  our  sample 

design  is  fi’om  0^  ^/cm^  to  0^  AjcYO^  . 

Since  the  total  specific  resistance  is  at  about 
-3  X  order,  the  contact  resistance  is  a  big 

issue  in  order  for  the  BIPOLITT  diode  has  a  good 
performance.  Note  that  the  DC  Ohmic  contact  on  GaAS 
with  value  less  than  3  x  lO^^Q-cm^  is  possible  to  be 
made  in  today’s  technology  [10-12].  In  the  meantime,  since 
our  BIPOLITT  diode  is  designed  at  high  jfrequencies,  the 
shimt  capacitance  becomes  important  and  decreases  the 
total  resistance  of  the  ohmic  contact.  In  our  design  example, 

the  capacitance  is  about  30  x  10”^  7^/cm^,  and  at  300 

GHz  the  impedance  is  about  .8  x  10”^Q  •  cm^. 
Therefore,  the  total  resistance  of  the  Ohmic  contact  at  300 
GHz  will  be  about  0.8  x  10"^Q  •  cm^,  which  allows  our 
BIPOLITT  diode  to  perform  with  sufficiently  specific 
negative  resistance  for  many  applications.  Due  to  the 

jfrequency  dependence  of  y  (similar  to  that  of  P* ),  the 
BIPOLITT  is  a  wide  band  transit  time  device  as  we  can  see 
&om  Fig.  4. 


V.  THE  BL\SING  SCHEMES 

In  this  section,  we  would  like  to  discuss  two  possible 
biasing  circuits  for  our  BIPOLITT  diode  as  shown  in  Fig.  5. 
It  is  easier  to  explain  by  employing  the  BIT  symbol. 
Basically,  both  can  be  view  as  a  Dc  current  biased  BIT 
operating  at  base  Ac  floating  mode,  the  Ac  signal  drops 
across  the  emitter  and  collector. 


(b) 

Fig.  5  The  biasing  circuits  for  BIPOLITT  diode,  (a):  The 
voltage  across  the  drifl  region  is  adjusted  by  a  variable 
resistor  to  control  the  efficiency  and  power  of  the  diode; 
(b):  Without  the  variable  resistor,  the  voltage  across  the 
drift  region  is  self-sustained. 


In  Fig.  5(a),  the  BIPOLITT  diode  is  shown  as  a  BIT.  An 
ideal  current  source  Is  is  used  to  bias  the  BIPOLITT  diode 
at  DC  condition.  A  RF  choke  (RFC)  is  used  to  allow  the 
DC  bias,  while  preventing  the  current  source  fi’om 
degrading  the  Ac  signal. 

A  variable  resistor  is  used  to  allow  voltage  dropping 
across  the  collector-base  junction  to  be  changeable, 
therefore,  we  can  adjust  the  efficiency  of  the  BIPOLITT 
diode. 

Taking  our  designed  BIPOLITT  diode  as  an  example.  The 
specific  negative  resistance  is  about  -3  x  10”^  Q  •  cm^ .  In 
order  to  match  a  3Q  load,  the  area  for  our  BIPOLIT  diode 
is  about  X  1 0  ^  cm^ .  Biased  at  Aj cm^ ,  the  current 
source  will  be  0nL4.  Meanwhile,  according  the 
parameters  shown  in  Fig.  4,  the  DC  common-emitter 


current  gain  p  can  be  made  to  be  about  100  400. 

Therefore,  the  current  flowing  through  the  variable 
resistance  R  is  0.1  ~  0.025  11l4  .  If  we  required  the  DC 
voltage  across  the  collector-base  junction  as  3  volts  to 
have  a  good  ejBficiency,  the  R  should  be  adjusted  to  be 
30~120kQ.  At  the  same  time,  the  specific  negative 
resistance  of  the  drift  region  (the  collector-base  junction)  is 
about  —5  X  10"^Q*  cm^/l  X  10”^  cm^  =  5Q,  which  is 
negligible  when  compared  to  the  much  larger  value  R. 
Therefore,  the  Ac  signal  will  not  go  through  the  DC  biased 
R  branch,  but  totally  flow  through  the  BIPOLITT  diode, 
which  makes  the  BIPOLITT  diode  as  a  two-terminal 
device. 

The  difference  between  Fig.  5(a)  and  (b)  is  that  in  Fig. 
5(b),  the  variable  resistor  R  is  eliminated.  Actually,  the 
effect  of  the  dropout  of  the  R  is  significant.  First,  it  greatly 
simplifies  the  fabrication  of  diode  and  circuit.  Second, 
without  R,  the  multi-layer  structure  (series  of  BIPOLITT 
diodes)  will  be  possible  for  BIPOLITT  integration  to  realize 
a  high  power  device.  However,  there  exists  a  big  problem 
for  our  AlGaAs/GaAs  structure:  the  break  down  voltage  of 

the  device  (similar  to  the  BVceo  of  BJTs)  is  as  low  as 
about  2  votls  since  our  drift  region  is  at  0.2  fjm  order. 

Biased  at  0^  to  0^  Ajc^VCi  DC  current  density, 

the  voltage  across  the  AlGaAs/GaAs  heterojunction  (the 
emitter-base  junction)  is  about  1.4  votls;  therefore,  there  is 
only  about  0.6  volts  voltage  left  for  the  drift  region  (the 
base-collector  junction),  which  is  not  enough  for  a 
reasonable  efficiency. 

Fortunately,  with  today’s  mature  technology  of  HBTs,  we 
can  figure  out  several  ways  to  overcome  this  drawback  by 
modifying  the  basic  BIPOLITT  structure.  For  example,  we 
can  insert  a  thin  insulting  region  between  the  base  and 
collector,  or  employ  the  heterostructure  for  the  base- 
collector  junction,  too.  Many  works  [refi*eneces???]  has 
shown  those  two  techniques  can  improve  the  breakdown 
voltage  substantially.  If  somehow,  we  can  increase  the 
breakdown  voltage  to  about  4  volts,  which  leaves  2.6  volts 
for  the  drift  region  to  have  a  good  efficiency.  Of  course,  we 
also  can  use  other  materials  with  higher  breakdown  voltage 
instead  of  GaAs.  Nevertheless,  the  biasing  circuit  shown  in 
Fig.  5(b)  worthies  attentions. 


VI.  DISCUSSION  AND  CONCLUSION 

Our  BIPOLITT  diode  has  a  comparable  performance  to 
the  RTDs  at  millimeter  and  submillimeter  wave  fi*equency 
range.  Meanwhile,  because  BIPOLITT  does  not  have 
negative  resistance  at  lower  firequency,  it  will  not  have  the 
problem  of  spurious  oscillation  and  bias  circuit  design 
problem  [13,14]  as  found  in  RTD.  It  can  also  be  integrated 
in  series  without  the  dc  instability  problem  of  RTD  [15].  It 
therefore  has  the  potential  to  be  a  higher  power  device  than 
RTDs,  and  be  as  one  of  the  best  solid  state  sources  in  the 
submillimeter  and  deep  millimeter  wave  fi*equency  range.  A 
large  signal  analysis  shows  the  efficiency  of  the  BIPOLITT 


is  about  10%  [16],  which  is  comparable  to  that  of  the  RTD 
of  about  10%  in  this  frequency  range. 

As  mentioned  earlier,  the  major  difference  between  our 
BIPOLITT  diode  and  the  commercially  designed  BJTs  or 
HBTs  [4-6]  is  that  our  BIPOLITT  diode  has  a  relatively 
longer  base,  which  introduces  a  phase  delay  close  to  80°. 
In  fact,  the  shorter  base  in  the  commercially  designed  BJTS 
or  HBTs  is  one  of  the  major  reasons  why  the  expected 
negative  resistance  phenomena  are  not  observed  in  them. 
First  of  all,  due  to  die  shorter  bases,  the  commercial  BJTS 

or  HBTs  generally  have  base  phase  delays  <  45° ,  which  is 
not  enough  for  a  good  transit  time  delay.  Our  simulation 

[16]  found  out  that  at  about  30°  and  90°  base  phase 
delays,  the  efficiency  of  a  HBT  as  a  negative  resistance 
device  are  about  1%  and  5%,  respectively,  while  reaches 

about  10%  at  80°  base  phase  delay.  Second,  in  order  to 

have  the  optimum  360°  total  phase  delay,  the  shorter  the 
base,  the  longer  the  collector  (the  drift)  region  is.  When  the 
collector  becomes  longer,  the  diffusion  effect  will  degrade 
the  performance  more  seriously,  which  accounts  for  one  of 
the  reasons  why  our  BIPOLITT  diode  with  a  shorter  drift 
region  is  better  for  the  negative  resistance  device.  In  the 
meantime,  in  the  commercial  BJTs  or  HBTs  designed  to 
operate  at  microwave  fi^equency,  the  negative  resistance  in 
the  drift  region  is  usually  cancelled  out  by  the  parasite 
emitter,  base  and  contact  resistance.  Our  BIPOLITT  diode 
is  designed  for  millimeter  and  submillimeter  wave 
fi*equency  range,  where  the  parasite  resistance  are  greatly 
decreased  by  the  no  longer  negligible  shunt  capacitance; 
therefore,  the  BIPOLITT  diode  has  the  potential  to  be  a 
negative  resistance  device  at  millimeter  and  submillimeter 
wave  fi-equencies  if  carefully  designed. 

In  conclusion,  we  explores  the  potential  of  making  a  good 
transit  time  device  based  on  the  HBT  structure.  The 
proposed  BIPOLITT  has  the  great  potential  to  become  the 
solid  state  source  at  deep  millimeter  and  submillimeter 
wave  fi*equencies.  It  can  find  applications  in  millimeter- 
wave  oscillators,  amplifiers,  and  mixers. 

This  work  was  supported  by  the  UCLA  JSEP  program. 
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A  Post-Shrink  MOSFET  Design  Utilizing 
Schottky  Tunneling  Meehanism 

J.Y.  Liao  and  D.S.  Pan 


Abstract — We  investigate  a  very  small  transistor  design 
utilizing  Schottky  tunneling  mechanism.  The  structure  with  a 
metal/semiconductor  Schottky  junction  replacing  the  source  of 
a  conventional  MOSFET  is  extensively  analyzed  and 
simulated.  A  number  of  simple  analytical  formulae  to 
implement  the  design  are  derived.  Very  good  transistor 
characteristics  are  obtained 

Index  Terms — MOSFET,  tunneling 

I.  Introduction 

As  the  device  development  advances  toward  sub-lOOnm 
regime,  there  exist  a  number  of  significant  technical 
barriers  to  the  extension  of  current  technology.  Material 
properties  and  fundamental  operating  principles  of 
MOSFETs  will  soon  become  the  limiting  factors  of  device 
scaling  [1].  Furthermore,  non-scaling  parameters  such  as 
threshold  voltage  {Vj.)  and  subthreshold  current  ( 7^^,^ )  are 
still  challenging  the  sub-lOOnm  technology.  Moreover, 
fabrication  technology,  in  an  effort  to  keep  up  with  the 
extraordinary  advancement  of  integrated  circuits  (IC),  is 
facing  serious  problems  such  as  random  processing  defects 
and  interconnection.  To  overcome  these  difficulties, 
research  efforts  have  been  directed  toward  the  investigation 
and  development  of  alternative,  non-conventional 
MOSFET  structures  [2].  These  devices  should  offer  not 
only  high  transconductance  but  also  relatively  simple 
fabrication  procedures  to  minimize  parasitic  effects. 
Structures  such  as  Schottky  barrier  source  and  drain 
MOSFETs  have  been  proposed  [3-5]. 

The  voltage  loss  at  the  Schottky  barrier  between  the 
source  and  the  channel  has  been  a  major  concern  in  the 
development  of  Schottky  source/drain  devices  for  low 
voltage  applications.  As  device  engineers  are  confronted 
with  complications  associated  with  ultra-small  transistors, 
experiments  have  demonstrated  a  favorable  possibility  of 
velocity  overshoot  in  very  short  devices  under  low  voltage 
operation  [6].  For  devices  with  gate  length  below  0.1  pm, 
electron  speed  at  the  source  edge  has  become  the  most 
important  performance  factor  [6,7].  Electrons  penetrating 
the  Schottky  barrier  possess  significantly  higher  speed  than 
the  velocity  saturation  value,  since  they  are  oriented 
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directly  into  the  channel.  Also,  contact  resistance  is 
eliminated  because  electrons  enter  the  channel  directly 
from  the  metal.  Thirdly,  because  tunneling  phenomenon  is 
a  majority  carrier  effect,  issues  such  as  inversion  layer 
broadening  can  be  neglected,  and  the  turn-off  behavior  of 
transistors  is  characterized  by  thermionic  tunneling. 
Moreover,  no  gate  voltage  is  needed  to  change  the  surface 
potential  and  to  offset  the  depletion-layer  charge;  space 
charge  is  induced  by  the  comer  electric  field  at  the  source- 
channel  interface,  and  independent  of  channel  doping.  The 
combination  of  the  above-mentioned  offers  a  beneficial 
tradeoff  for  the  voltage  drop  loss  at  the  source  edge. 

Approaches  in  utilizing  tunneling  barrier  in  MOSFET 
operation  have  been  reported  [4,8].  In  these  MOSFETs 
where  tunneling  injectors  are  employed,  the  source  is 
typically  a  reverse-biased  Schottky  barrier.  The  Schottky 
barrier  depletion  width  at  the  reverse-biased  condition  is 
dependent  on  the  gate  bias;  this  dependence  results  in  the 
three-terminal  transistor  effect:  the  drain  current  being 
controlled  by  the  gate  bias.  In  other  words,  when  a  positive 
gate  voltage  is  applied,  the  tunneling  barrier  at  the  source- 
channel  interface  becomes  narrower,  resulting  in  the 
enhancement  of  tunneling  mechanism. 

In  research  work  where  a  Schottky-semiconductor 
interface  has  been  employed  to  both  the  drain  and  source  of 
the  MOSFET,  the  reported  drain  current  is  low  (in  the  order 
of  pA/pm),  and  the  transistor  does  not  leave  triode  mode 
[5].  Furthermore,  subthreshold  behavior  of  these  devices 
is  not  characterized.  Although  Schottky  barrier  transistors 
do  offer  a  fabrication  advantage  over  conventional 
MOSFETs,  researchers  have  to  address  the  enhancement  of 
transconductance  and  transistor  characteristics,  as  these  are 
the  major  concerns  for  the  development  of  tunneling 
transistors  for  IC  applications. 

II.  Device  Structure 

The  stmcture  proposed  in  this  research  work  is  similar 
to  a  MOSFET,  with  metallic  termination  at  the  source.  The 
significant  features  are  the  physical  dimensions:  the  gate 
length  varies  from  3 00 A  to  5 00 A  in  order  to  enhance 
Schottky  tunneling  mechanism,  and  the  thickness  of  the 
active  layer  is  generally  X  (100 A  in  the  stmcture  proposed 
in  Fig.  1)  of  the  gate  length  for  large  current  handling 
capability.  This  stmcture  consists  of  an  n-doped 
semiconductor  as  the  active  layer,  an  insulated  gate,  and 
intrinsic  substrate  which  is  necessary  to  provide  electrical 
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isolation  between  neighboring  electrodes.  To  characterize 
tunneling  mechanism  in  modeling  and  simulation,  silicon- 
germanium  parameters  are  employed.  For  better  device 
isolation  and  reduction  in  substrate  leakage  and  parasitic 
capacitances  [2],  a  backside  oxide  is  sandwiched  between 
the  active  thin  film  and  the  substrate,  similar  to  a  Silicon- 
on-Insulator  (SOI)  MOSFET. 


Oxide 


Figure  1.  Schematic  structure  of  a  Schottky  source 
MOSFET.  ^^^(=50A)  is  the  gate  oxide  thickness, 
/^-(-lOOA)  the  active  layer  thickness,  and  L  (=3  00 A) 
the  gate  length. 

III.  Simulation  Method 

ATLAS  2-D  device  simulator  program  by  Silvaco  is 
utilized  in  the  simulation  of  proposed  device  structures. 
The  simulator  does  not  provide  a  model  for  Schottky 
tunneling  but  for  interband  tunneling.  However,  because  of 
the  fact  that  both  tunneling  mechanisms  are  characterized 
as  functions  of  local  electric  fields  [9],  the  interband 
tunneling  method  can  be  modified  to  accommodate  the 
modeling  of  Schottky  tunneling  mechanism. 

Two-dimensional  numerical  simulation  is  performed 
on  a  Schottky  barrier  MOSFET  with  metallic  material 
replacing  both  the  conventional  source  and  drain.  Because 
electrons  generated  at  the  source  edge  have  to  overcome  the 
Schottky  barrier  at  the  drain  edge,  the  device  suffers  a 
delayed  tum-on.  A  voltage  applied  to  lower  the  drain 
barrier  is  required  before  large  amount  of  electrons  can  be 
collected  by  drain.  The  transistor  enters  saturation  mode  at 
a  higher  drain  voltage,  and  the  delay  in  turn-on  is 
unfavorable  toward  low  voltage  operation.  In  order  to 
achieve  a  better  transistor  characteristic  where  the  device  is 
saturated  at  a  low  drain  bias,  the  transistor  with  the  drain 
remaining  a  perfect  Ohmic  contact  is  extensively  simulated. 

IV.  Analysis  OF  Tunneling  Current 

Generally,  the  Schottky  tunneling  current  density 
discussed  in  previous  research  work  [4]  is  expressed  in  an 
integral  form,  and  the  dependence  of  current  on  transverse 
momentum  was  not  considered  in  the  derivation.  We 
attempt  to  include  the  measure  of  significant  range  of 
transverse  momentum  in  our  discussion,  and  provide  a 
simple  analytical  expression  for  the  Schottky  tunneling.  In 
this  section  we  will  discuss  the  derivation  of  tunneling 


currents,  current  dependence  on  electric  fields, 
approximation  of  electric  fields,  and  subthreshold  behavior. 

Derivation  of  Tunneling  Currents 

To  better  understand  and  analyze  the  tunneling 
currents,  we  categorize  the  total  current  into  four 
components  according  to  the  energy  range  within  which 
tunneling  occurs:  1.)  Schottky  tunneling  current  for 
tunneling  between  the  metal  Fermi  level  and  the  edge  of 
semiconductor  valance  band;  2.)  Band-to-band  tunneling 
current  for  tunneling  from  the  valance  band  edge  to 
infinity;  3.)  Thermionic  tunneling  current  for  tunneling 
between  metal  Fermi  level  to  the  edge  of  the  conduction 
band;  and  4.)  Over-barrier  thermionic  tunneling  current  for 
tunneling  from  the  conduction  band  edge  to  infinity.  For 
significant  interband  tunneling  to  occur,  a  large  drain 
voltage  must  be  applied  to  reverse  bias  the  junction  so  that 
the  conduction  band  is  lowered  below  Fermi  level. 
Moreover,  electrons  must  travel  a  distance  to  reach  valance 
band  edge  before  starting  to  tunnel,  and  the  electric  field  at 
the  location  for  interband  tunneling  to  occur  should  also  be 
sufficiently  high.  Because  of  these  constraints,  interband 
tunneling  is  not  discussed  in  detail  in  this  paper. 

Consider  a  metal- semiconductor  interface  as  illustrated 
in  Fig.  2.  Electrons  with  energy  above  metal  Fermi  level 
are  trying  to  tunnel  through  the  barrier  width  and  enter  the 
conduction  band.  When  a  small  bias  is  applied  to  the  drain, 
the  barrier  at  the  drain  side  is  lowered,  and  electrons  are 
collected.  When  the  drain  bias  reaches  approximately  the 
difference  between  the  metal  Fermi  level  and  the 
semiconductor  valance  band  edge,  the  integration  of  the 
number  of  electrons  constitutes  the  total  Schottky  tunneling 
current.  Additional  depletion  width  narrowing  beyond  this 
drain  voltage  becomes  negligible,  and  further  increase  in 
drain  bias  does  not  result  in  significant  current  increase. 
Therefore,  the  device  enters  saturation  mode  beyond  this 
point. 


Figure  2.  Schematic  diagram  of  Schottky  barrier  and 
the  definitions  used  in  the  current  derivation. 
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From  the  derivation  of  electrons  tunneling  through  a 
triangular  barrier,  and  the  theory  of  tunneling  with 
momentum  being  conserved  [2,9],  the  incident  current 
density  is  expressed  as 


The  integration  to  determine  Schottky  current  density 
is  taken  over  the  energy  range  from  the  metal  Fermi  level 
( is  defined  0  at  this  level)  to  the  edge  of  the 
semiconductor  valance  band  {E^- ).  A  close  form  of 
the  current  expression  can  be  closely  estimated  by  Taylor 
expansion.  Numerical  analysis  shows  that  the  simplified 
expression  is  only  slightly  deviated  from  the  original 
equation,  and  it  provides  a  good  approximation  for 
Schottky  tunneling  current.  For  large  bandgap  material 
(e.g.  Si)  where  tunneling  is  difficult,  a  large  electric  field  is 
required  to  achieve  a  large  current  density  for  useful 
applications  (Table.  1). 

For  Si-  or  SiGe-Al  interfaces,  Schottky  barrier  height  is 
roughly  %  of  the  semiconductor  bandgap.  We  can  further 
simplify  the  expression  of  as 


'tSCH 
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This  is  a  simple  formula  to  implement  the  transistor  design 
utilizing  Schottky  tunneling  mechanism  at  the  source- 
channel  edge. 


Max^o 

fe) 

Si 

SiGe 

2 

500 

1.22e6 

3 

2.7e4 

4.66e6 

4 

2.3e5 

9.86e6 

Table  1 .  Schottky  tunneling  current  density  for  Si- 
and  SiGe-Al  parameters  under  maximum  electric 
field. 


Current  Dependence  on  Electric  Fields 


From  simulations,  the  magnitude  of  the  parallel  electric 
field  is  varying  with  position,  and  can  be  approximated 
with  a  simple  exponentially  decaying  function  with  a  decay 
parameter  of  : 


^(x)  =  ^oexp 


In  order  words,  when  the  distance  of  tunneling  reaches 
beyond  oxide  thickness,  the  electric  field  becomes  too 
small,  and  tunneling  current  generated  under  this  electric 
field  is  negligible.  Changing  the  integration  variable  at 

X  =  -^ ,  we  can  re-write  the  current  density  as 
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Comparison  of  Schottky  tunneling  currents  with  and 
without  electric  fields  being  dependent  on  position  is 
illustrated  in  Fig.  3. 


Figure  3.  Calculated  maximum  Schottky  tunneling 
currents  with  SiGe  parameters.  Vj^  =  0.5V  .  =  50A. 

Dependence  of  electric  field  on  position  is  included  in 
dashed  line. 

Thermionic  Currents 


Electrons  with  all  energy  levels  are  also  tunneling 
through  Schottky  barrier.  These  electrons  contribute  to  the 
thermionic  current.  The  thermionic  current  consists  of  two 
components:  thermionic  tunneling  current  integrated  from 
the  metal  Fermi  level  to  the  edge  of  the  Schottky  barrier 
height,  and  over-barrier  thermionic  current  integrated  from 
the  edge  of  the  Schottky  barrier  height  to  infinity.  Over¬ 
barrier  thermionic  current  density,  Jovither  » given  by  [2] 

J  -  A* expf-  ^ 

oviTHER  —  ^  J 

where  A*  is  the  effective  Richardson  constant  for 
thermionic  emission.  For  device  operated  at  room 
temperature,  Jqvu'her  1^  negligible. 

The  thermionic  tunneling  current  component  is 
proportional  to  the  quantum  transmission  coefficient 
multiplied  by  the  occupation  probability  in  the 
semiconductor  and  the  unoccupied  probability  in  the  metal. 
The  tunneling  probability,  however,  is  significantly  smaller 
for  electrons  with  energy  below  Fermi-level,  although  they 
are  to  tunnel  through  a  narrower  barrier.  This  small 
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amount  of  current  is  also  dependent  on  the  electric  field. 
The  current  density  expression  for  the  thermionic 
component  is  derived  with  a  similar  approach  to  that  of 
Schottky  tunneling  current,  with  modification  of  the 
integration  boundaries  and  tunneling  constants. 


where 


A^Jlrn 

3qfi 
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and  F{E^)  is  the  occupancy  factor  given  by: 
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The  integration  is  taken  (with  positive  upwardly) 
from  Fermi  level  (defined  0)  to  the  edge  of  the  conduction 
band.  is  dominant  under  small  gate  bias,  and 

contributes  to  the  total  current  at  high  gate  voltages. 


Electric  Field  (MV/cm) 

Figure  4.  Total  current  (  =  Ifscfi  +  hwER  )  versus 

electric  field  with  SiGe-Al  parameters.  is  lOOA  in 
IsttlOO,  and  50A  in  IsttSO. 

Threshold  Voltage 

For  conventional  MOSFETs,  the  threshold  voltage, 
Vj,  is  defined  as  the  gate  bias  when  diffusion  current  is 
equal  to  the  conduction  current.  In  Schottky  source/drain 
MOSFETs,  these  are  the  thermionic  current  and  the 
Schottky  tunneling  current,  respectively.  The  ON  current, 
7^^ ,  can  be  considered  as  the  current  flowing  when 
Schottky  tunneling  current  is  equal  to  thermionic  current, 
or  when  current  is  increased  a  decade  higher  than  the 
thermionic  current  under  no  gate  bias.  With  numerical 
approximation,  a  close  form  of  expressed  as 
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is  integrated 


Consider  the  first  definition.  While  7^ 
up  to  the  distance  roughly  equal  to  t^^  (obtained  from 
simulation  and  two-dimensional  analysis),  is 

integrated  along  the  entire  active  layer  thickness.  Solving 
for  ^0  when  7,7^^^  =  fscn  ^  ^he  following 

expression: 
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Numerical  solution  with  geometry  parameters  yields 

t  ^SCH 

4 

If  we  consider  as  an  order  of  magnitude  higher 
than  at  =  OF  ,  equating  the  two  currents  yields 
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The  solution  with  proper  geometric  parameters  gives 
similar  results  as  in  first  case.  This  is  the  threshold  electric 
field,  and  the  corresponding  threshold  voltage  is  given  by 
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^ox  ^SCH  _|_  Y 


FB 


I^CltTHER) 

log(ItSCH) 


E-field  (MV/cm) 

Figure  5.  Calculated  currents  versus  electric  field. 
C  =50A. 


Electric  Fields 

To  understand  the  electric  field  distribution  is  the  most 
important  step  in  analyzing  Schottky  tunneling  mechanism, 
since  it  is  characterized  as  a  function  of  local  electric  field. 
The  electric  field  (4)  the  comer  of  the  source-channel 
interface  below  gate  oxide  strongly  affects  the  current 
density.  For  simplicity,  the  electric  field  is  expressed  as: 


4-^ 
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A  correction  factor,  k  ,  necessary  to  modify  the  field 
expression,  is  dependent  on  the  geometries  of  and  . 
Parallel  electric  field  distribution  against  distance  along  the 
channel  under  different  gate  biases  is  plotted  in  Fig.  6(a). 
For  all  figures  plotted  against  the  distance  along  the 
channel,  the  device  dimensions  are  described  in  Fig.  1 
unless  otherwise  noted.  Cross-section  is  obtained  lA  into 
the  active  thin  film  from  the  Si/Si02  interface. 


Distance  along  Channel  (um) 


(a)  =0,  steps  from  0  to  2V. 


(b)  =0,  Vj^  steps  from  0  to  2V. 

Figure  6.  Simulated  parallel  electric  fields  along  channel. 

The  varying  voltage  increases  evenly  with  0.5V  per  step. 

The  electric  field  is  not  only  modulated  by  but 
affected  by  as  well.  Figure  6(b)  illustrates  the  drain 
influence  on  the  parallel  electric  fields  at  =2V.  The 

phenomenon  is  possibly  resulting  from  the  two- 
dimensional  effect.  To  understand  the  dependence  of 
drain  voltage  on  the  maximum  comer  electric  fields,  we 
employ  the  two-dimensional  analytical  model  for  fully 
depleted  SOI  MOSFETs  proposed  by  Young  [16],  and 
modify  the  boundary  conditions  to  accommodate  the 
Schottky  barrier  source.  It  should  be  noted  that  when 
drain  is  biased  beyond  depletion,  space  charge  effects 
should  also  be  considered  and  included  in  the  calculation. 
The  variation  of  electric  field  due  to  generated  carrier 
space  charge  is  more  significant  at  smaller  drain  biases 
(Fig.  7). 


The  derivative  of  potential  distribution  with  respect  to 
the  position  in  the  tunneling  direction,  yields  the  parallel 
electric  field  expression.  For  the  device  stmcture  in  Fig.  1 
with  a  backside  oxide  of  thickness  of  50A 

sandwiched  between  the  active  silicon  layer  and  the 
substrate,  the  maximum  electric  field  can  be  approximated 
by 

where  =  0.2*^„  from  simulations. 

If  is  30A, 

^<1/  +  0.054F,;)+ 

where  = 


(a)  K^,=0.2V. 
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(b)  ^as  =2V 


Figure  7.  Calculated  potential  distribution  (V)  versus 
distance  along  channel  under  Si/Si02  interface  from 
two-dimensional  model.  Space  charge  effect  is  included 
in  .  Gate  is  unbiased. 
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1 

1.87 

1.86 

2.15 

2.17 

1.5 

2.36 

2.41 

2.68 

2.83 

2 

2.79 

2.95 

3.15 

3.49 

Table  3.  Simulated  and  calculated  maximum  corner 
electric  fields  for  different  ^,,^of  an  SSOD  device. 

=  0.5V .  A  deviation  of  the  calculated  values  from 
the  simulated  is  observed  at  high  gate  biases.  This  is 
because  of  electric  fields  being  less  dependent  on  V^^. 
beyond  saturation. 


Figure  8.  Simulated  (dots)  and  calculated  (rhombuses) 
maximum  comer  electric  fields  of  an  SSSD  device. 
/„,=50A. 
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Figure  9.  Calculated  dependence  of  on  and 

ta.-  f^.-5000A. 

Figure  8  illustrates  the  simulated  and  calculated 
versus  V^^.  under  different  V^^. .  A  significant  difference 
between  the  simulated  and  calculated  4  is  observed  at 
small  V^^. .  Because  at  small  ,  electrons  penetrating  the 
source  barrier  are  blocked  by  the  Schottky  barrier  at  drain, 
the  accumulation  of  electrons  confined  by  the  barrier  well 
(Fig.  11(a))  lowers  the  electric  field  at  the  source-channel 
interface.  This  phenomenon  is  not  included  in  the  two- 
dimensional  derivation.  However,  the  equations  still 
provide  a  close  estimate  of  when  the  device  saturates. 
In  fact,  the  lowered  electric  field  at  small  results  in 
small  current  flow,  and  therefore  superior  turn-off 
behavior. 

Intuitively,  would  be  inversely  proportional  to  ; 
i.e.,  when  is  reduced  from  50A  to  30A,  should 
increase  accordingly.  However,  simulation  and  calculation 
do  not  agree  with  the  ideal  expectation.  This  is  because  of 
the  fact  that  in  additional  to  ,  other  geometric  parameters 
such  as  channel  thickness  ( )  and  backside  oxide  thickness 
( )  are  also  affecting  the  magnitude  of  .  These  terms 
are  expressed  in  ,  and  grouped  as  the  correction  factor 
K  .  The  following  figures  show  the  geometric  dependence 
of  a^.  r^,  a  parameter  which  characterizes  two- 

dimensional  effect,  is  the  product  of  a^^^^and  the  channel 
length.  The  larger  ,  the  lesser  dependence  of  potential 
and  electric  fields  on  two-dimensional  geometry. 
Numerical  calculation  shows  that  with  given  channel  length, 
short  channel  effect  is  significant  when  or  is  reduced 
(Fig.9). 

Integration  Distance 

Gate  oxide  thickness  is  not  only  an  important 
parameter  for  approximating  maximum  electric  field  at  the 
comer  of  the  source-channel  interface,  but  also  assumed  as 
the  integration  distance  for  current  estimation  along  the 
direction  perpendicular  to  the  Si/Si02  interface. 
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The  tunneling  current  is  obtained  by  integrating 
the  current  density  along  the  distance  perpendicular  to  the 
Si/Si02  interface.  The  determination  of  depends 
strongly  on  the  decay  of  .  As  decays  to  a  certain 
value  when  tunneling  under  this  electric  field  becomes  less 
significant,  this  will  be  the  upper  integration  limit.  There 
are  a  few  possible  estimates  of  ty.  It  is  simplicity  to 
assume  that  the  electric  field  decays  with  the  decay  factor 
equivalent  to  the  gate  oxide  thickness;  therefore,  ty  equals 
to  .  This  assumption  provide  an  acceptable 
approximation  of  drain  current  for  =50 A.  When  is 
reduced,  we  expect  much  larger  drain  current  due  to  the 
increase  in  electric  field.  However,  as  ty  is  limited  by  , 
the  total  current  may  not  be  much  higher  than  that  for 
thicker  gate  oxide. 

V.  Simulation  Results 

The  analysis  is  compared  with  the  simulation  results 
for  a  better  understanding  of  Schottky  tunneling 
mechanism.  Structures  with  Schottky  source  and  Schottky 
or  Ohmic  drain  are  simulated.  Effects  of  geometric 
dependence  are  also  included.  I-V  characteristics, 
threshold  voltage,  transconductance  and  cutoff  frequencies 
are  discussed  in  detail. 

Transistor  Characteristics  of  Schottky  barrier  Source  and 
Drain  (SSSD)  MOSFETs 

The  structure  with  Schottky  barrier  source  and  drain 
(SSSD)  was  simulated.  The  drain  current  is  plotted  against 
a  sweeping  drain  voltage  in  Fig.  10,  with  gate  bias  stepping 
from  0  to  2V.  A  delayed  tum-on  and  the  near-saturated  I-V 
curves  are  the  major  characteristics  observed  in  this 
simulation.  The  delay  in  tum-on  is  due  to  the  voltage  loss 
at  the  Schottky  barrier  between  the  drain  and  the  channel, 
when  the  high  speed  electrons  generated  by  tunneling 
mechanism  at  the  source-channel  interface  have  to 
overcome  the  barrier  before  entering  drain  (Fig.  11(a)). 
Although  tunneling  width  narrows  as  increases, 
electrons  penetrating  the  source  barrier  are  blocked  by  the 
drain  barrier  when  is  small.  The  drain  barrier  has  to  be 
lowered  for  a  large  amount  of  electrons  to  be  collected  by 
drain  (Fig.  11(b)).  A  slope  in  transistor  characteristics 
beyond  saturation  is  observed  in  the  simulation.  This  is 
attributed  to  interband  tunneling,  when  electrons  begin  to 


tunnel  from  valance  band  directly  into  conduction  band  at 
sufficiently  high  drain  voltage. 

Figure  10.  Simulated  (solid  lines)  and  calculated  (dotted 
lines)  transistor  characteristics  of  an  SSSD  MOSFET. 
steps  evenly  from  0  to  2V. 

The  resulting  electric  field  (Fig.  12(b))  at  F^=0  is 
significantly  lowered  due  to  the  confinement  of  electrons 
between  the  two  Schottky  barriers.  The  current  flowing  is 
the  thermionic  tunneling  current  under  high  electric  fields. 
It  is  also  observed  that,  although  the  drain  current  exhibits 
an  exponential  increase  after  tum-on,  the  device  will  not 
reach  saturation  until  sufficiently  high  drain  voltages  (IV  is 
necessary  in  this  simulation.)  For  it  is  necessary  to  lower 
the  conduction  band  edge  of  the  drain  edge  to  the  level  of 
the  metal  Fermi  level  before  electrons  can  be  collected. 

^ds,sA7'  approximately  equal  to  ;  i.e.,  after  the 
conduction  band  edge  at  drain  is  lowered  to  the  level  as  the 
valance  band  edge  at  source,  further  increase  in  drain 
voltage  will  not  result  in  significant  increase  in  current. 

(a) 

(b)  F,,.=0.5V. 

Figure  11.  Simulated  conduction  band  bendings  of  an 
SSSD  MOSFET.  increases  evenly  from  0  to  2V 
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(a)  =0,  steps  from  0  to  2V, 
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(b)  =0,  F^^.  steps  from  0  to  2V. 

Figure  12.  Simulated  parallel  electric  fields  of  an  SSSD 
MOSFET. 


Transistor  Characteristics  of  Schottky  Source  Ohmic 
Drain  (SSOD)  MOSFETs 
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(d)  F^^.  =2V,  F^^.  steps  from  0  to  2V. 

Figure  13.  Simulated  conduction  band  bendings  of  an 
SSOD  MOSFET  along  channel  under  different  bias 
conditions.  F^^,  is  included  in  (a)  and  (b).  The 
varying  voltage  increases  evenly  with  0.5V  per  step. 

To  improve  saturation  characteristics  for  low-power  IC 
applications,  simulations  were  performed  on  a  modified 
structure  with  drain  remaining  a  perfect  Ohmic  contact. 
The  profiles  of  conduction  band  bendings  along  channel 
are  illustrated  in  Fig.  13(a)  with  drain  unbiased.  As  gate 
bias  is  equal  to  the  flatband  voltage,  the  potential  barrier  is 
tall  enough  to  block  electrons  flowing  from  the  source  to 
the  drain.  When  gate  bias  is  more  positively  increased, 
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electrons  at  the  Si/Si02  interface  are  induced,  and  a  much 
thinner  tunneling  barrier  is  formed  between  the  source  and 
the  channel,  therefore  resulting  in  a  significantly  enhanced 
channel  conductance.  Under  this  condition,  current  starts 
to  flow  when  a  small  drain  bias  is  applied  (Fig.  13(b)). 

A  small  barrier  at  the  drain  which  results  from  the 
built-in  potential  between  two  differently  doped 
semiconductor  regions  has  to  be  overcome  before  current 
enters  drain  (Fig.  13(a)).  When  a  moderate  drain  bias  is 
applied,  current  starts  to  flow,  and  increase  linearly  with 

.  When  reaches  beyond  “  A/? )  further  increase 
in  F^will  not  result  in  linear  variation  in  current.  The 
device  enters  saturation  mode  at  this  The  device 

enters  saturation  when  the  drain  bias  reaches  approximately 
0.5V  (Fig.  15). 

It  is  also  observed  that  when  drain  is  biased  beyond 
certain  range,  its  effect  on  the  depletion  width  narrowing 
becomes  insignificant  (Fig.  13(d)).  This  is  because  of  the 
fact  that  the  dependence  of  electric  field  is  on  drain  voltage 
is  modulated  by  the  channel  length,  which  is  much  larger 
than  the  gate  oxide  thickness. 

The  effect  of  different  channel  thicknesses  is  also 
simulated  (Fig.  14).  However,  the  difference  in  drain 
current  cannot  be  attributed  only  to  the  same  thermionic 
tunneling  current  density  integrated  up  to  •  Fur  in 
reducing  ,  the  interdependence  of  geometric  parameters 
change,  and  the  potential  and  electric  field  distributions 
also  vary,  according  to  the  two-dimensional  analytical 
modeling.  Numerical  calculation  shows  that  the  comer 
electric  field  and  tunneling  current  are  enhanced  at  = 
50A. 


Figure  14.  Simulated  transistor  characteristics  of  a 
SSOD  MOSFET  with  germanium  as  the  active  layer, 
increases  evenly  from  0  to  2V  with  0.5V  per 

step.  The  only  varying  parameter  is  :  lOOA  (dotted 
lines)  and  50A  (solid  lines). 

It  is  also  noticed  that  a  large  amount  of  current  exists 
with  unbiased  gate.  This  could  be  explained  by  the  fact 
that  the  built-in  potential  at  the  drain-channel  interface 


induces  an  electric  field  which  might  be  high  enough  for 
interband  tunneling  to  occur  before  gate  bias  is  applied. 
This  provides  a  possible  explanation  of  drain  current  when 
F^,=2Vand  F^,=0V  (Fig.  13(c)). 

The  transistor  characteristics  are  illustrated  in  Fig.  15. 
The  solid  lines  depict  the  simulated  family  curves  and  the 
dots  are  the  calculated  drain  currents  with  maximum 
electric  fields  extracted  from  simulation.  With  these 
modifications,  it  is  observed  that  the  values  of  total  drain 
currents  calculated  with  the  theoretical  expression 
correspond  to  the  saturation  currents  in  Fig.  15.  It  should 
also  be  noted  that  when  the  electric  field  drops 

significantly  fast,  reaching  e~^  of  its  maximum  value  at  the 
location  approximately  )4  of  beneath  the  Si/Si02 
interface,  where  is  the  decay  factor.  Therefore,  the  total 
tunneling  current  under  this  condition  is  integrated  up  to  X 
of  .  These  calculated  values  are  included.  When  is 
smaller  than  drain  current  is  controlled  by  both 

and  .  As  the  device  enters  saturation  mode,  the 
dependence  of  the  drain  current  on  F^,.  becomes  less 
significant,  and  F^,  is  modulating  the  ON-OFF  currents  of 
the  transistor.  These  are  similar  to  the  typical  /-F 
characteristics  of  conventional  MOSFETs,  and  ideal  for 
low-power  IC  applications. 


Figure  15.  Simulated  transistor  characteristics  of  an 
SSOD  MOSFET.  Simulated  (solid  lines)  and  calculated 
(dotted  lines)  I-V  characteristics  with  F^^.  stepping 
evenly  from  0  to  2V. 

Subthreshold  Behavior 

By  defining  the  ON  state  at  the  condition  where 
Schottky  tunneling  current  equals  thermionic  current,  we 
derive  a  simple  approximation  for  threshold  voltage.  It  is 
expressed  in  terms  of  Schottky  barrier  height,  effective 
mass,  front  oxide  thickness,  and  flatband  voltage.  With 
SiGe-Al  parameters,  Fy^  is  roughly  estimated  to  be  0.5V, 
which  corresponds  to  the  ON  position  (Fig.  16).  A  F^^.  of 
approximately  IV  is  necessary  to  raise  the  drain  current  by 
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3  decades  when  F^^.=0.5V.  Figure  17  shows  the 
subthreshold  swing  versus  gate  bias  under  different  drain 
voltages.  At  higher  ,  it  is  more  difficult  to  turn  off  the 
device,  as  a  larger  .S'  is  required. 


Figure  16.  Calculated  subthreshold  characteristics  for  an 
SSSD  device  with  Vj^  stepping  from  0  to  2V. 


Figure  17.  Subthreshold  swing  versus  gate  bias  with 
stepping  from  0  to  2V. 

Optimization  of  Gate  Oxide  Thickness 

The  derivation  and  simulation  of  total  drain  current 
from  the  previous  sections  reveal  that  both  the  magnitude 
of  the  parallel  electric  fields  and  the  integration  of 
tunneling  current  in  the  direction  perpendicular  to  the 
Si/Si02  interface  are  strongly  dependent  on  It  is 

therefore  important  to  characterize  in  order  to  achieve 
optimum  drain  current  for  the  device.  Figure  18  provides  a 
comparison  of  maximum  electric  fields  and  tunneling 
currents  at  different  for  the  SSSD  devices,  and  the 
transistor  characteristics  for  an  SSSD  device  with  =30A 
are  illustrated  in  Fig.  19.  Because  of  the  Schottky  barrier 
at  drain,  SSSD  devices  exhibit  sufficiently  small  , 
until  the  drain  barrier  is  overcome  by  biasing.  Although 
higher  electric  fields  result  in  higher  current  densities  at 
smaller  at  ON  state,  it  becomes  more  difficult  to  turn 


off  the  SSOD  devices,  because  a  very  small  barrier  is 
easily  overcome  by  the  high  speed  electrons  entering  the 
channel.  There  exists  at  0  a  noticeable  amount  of 
current  increase  when  is  decreased  to  30A. 


Figure  18.  Calculated  transistor  characteristics  of  an 
SSSD  MOSFET  with  SiGe  parameters  for  ^^^=50A 
(dotted  lines)  and  30A  (solid  lines).  ^  id  hx  ■ 


Vds  (V) 

Figure  19.  Simulated  (solid  lines)  and  calculated  (dotted 
lines)  transistor  characteristics  of  an  SSSD  MOSFET 
with  stepping  evenly  from  0  to  2V.  =30A. 

VI.  Conclusion 

We  have  presented  a  small  transistor  design  utilizing 
Schottky  tunneling  mechanism.  The  simple  analytical 
current  expression  is  derived  from  the  tunneling  theory. 
Simulations  and  discussion  of  transistor  behavior  are 
performed  to  verify  and  modify  the  current  expression. 

The  drain  current  due  to  tunneling  at  the  Schottky 
barrier  source-channel  interface  is  estimated  by 
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^QsCH^g 


^0  exp 


^SCH 

^0 


where 


PscH  and  Q, 


SCH 


are  tunneling  constants,  ^qXIiq 


maximum  comer  electric  field,  and  tjj  the  integration 
distance  along  the  direction  perpendicular  to  the  Si/Si02 
interface. 

The  drain  current  is  in  the  0.2mA/|im  range  when  the 
transistor  enters  saturation  mode.  Gate  capacitance  is  also 
reduced  comparing  to  a  conventional  MOSFET;  therefore, 
the  cutoff  frequency  can  reach  approximately  150GHz. 
With  Ohmic  drain  (SSOD  devices),  the  saturation  voltage 
^ds,sAT  of  the  device  is  approximately  equal  to  the 


difference  between  the  metal  Fermi  level  and  the  valance 
band  edge  of  the  semiconductor  at  the  metal-semiconductor 
interface  (~0.5V).  However,  with  Schottky  drain  (SSSD 
devices),  the  Schottky  barrier  height  also  contributes  as  an 
additional  term  in  (-1 V).  The  tunneling  parameters 


in  this  analysis  are  extracted  from  Al-SiGe.  Current 
fabrication  technology  provides  possible  tuning  of  Schottky 
barrier  height  to  semiconductor  by  employment  of  silicided 
material.  By  lowering  the  barrier  height,  Schottky 
tunneling  current  can  be  further  enhanced. 

The  characterization  of  is  also  advantageous  in  the 
tunneling  device.  As  subthreshold  current  is  dominated  by 
thermionic  tunneling,  is  determined  by  Schottky  barrier 
height  and  flatband  voltage,  and  independent  of  device 
doping.  Since  Cx  is  required  for  high  electric  field  only  at 
the  comer,  its  thickness  over  the  entire  channel  is  not 
critical  for  high  transconductance.  Current  fabrication 
technology  suggests  a  step-oxide  possibility.  The  critical 
thickness  is  only  necessary  for  a  small  area  covering  the 
source  comer  (Fig.  20).  Because  lower  electric  fields  result 
from  larger  at  the  drain  edge,  simulation  shows  that  the 
device  exhibits  a  better  saturation  characteristic:  flatter  IV 
curves  at  an  earlier  .  This  offers  another  superb 
advantage  in  utilizing  Schottky  tunneling  mechanism  in 
this  structure. 


Figure  20.  Schematic  structure  of  a  Schottky  source 
MOSFET  similar  to  Fig.  1,  except  gate  oxide  thickness 
is  thinner  near  the  corner  of  source-channel-  front 
Si/Si02  interface. 


VII.  Appendix 


Schottky  Tunneling  Current  Density 
The  tunneling  probability  is  of  the  form 
"  -4V2; 


T,  «  exp 


tm " 


3^qn^ 


2E^ 

expl-y- 


V  J 

where  m*  is  the  effective  mass,  the  bandgap  of  the 

semiconductor  and  ^  the  electric  field.  E  ,  the  measure  of 
the  significant  range  of  transverse  momentum,  is  given  by 


E  =■ 


3mn 


p\/2 


The  induced  current  per  unit  area  in  the  energy  range 
dE^dE^  is  given  by 


dJ^  = 


-^dE^dE, 


where  E^  is  the  energy  associated  with  momentum  in  the 
tunneling  direction.  The  tunneling  current  per  unit  area  can 
be  obtained  with  the  above  equations. 


qm* 


00 

^T,(FAE)-FAE))dEdE, 


where  F^.(£)and  F^(£)  are  the  occupancy  factors.  Since 
the  device  is  reverse  biased,  F^{E)^  0,  and  F^{E)=  1. 
The  integral  becomes 

Eg-q<i>,sb 

K 


^tSCH  ~  2 


^-A  ^ 

^SCH 


where  tunneling  parameters  ,  and  K  are  given 

by 

\3/2 


^SCH  ~ 


^SCH  - 


4^^2m"(q^,,+E,f 
3qh 

4^|2qn 

37tm’''^{q^,^+E^) 


1/2  > 


K  = 


qm^ 


Material  parameters  in  these  tunneling  constants  are 
separated  from  the  varying  term  E^ ,  grouped  and  defined 
as 


'  >  Qsch  ~  ' 


A^lqh 


3^^  37mi  {q(j),f,) 

The  incident  current  density  is  re-written  as 
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Qsch^ 


where  is  the  Schottky  barrier  height.  A  close  form 
expression  is  provided  as 


r  ^  I  I  ^  SCH 

-  2  ^ 


-p. 


3/2  ^ 


.  ^  N./2 

1  + 


<1^.4 


dE^ 


f  p  \ 

^SCH 

t: 

exp 

V  ^0  / 

V  V 

(  ^sb  ^SCH  ■•■  ^  ^SCH  ( ^ ^0  )  ^ ^^C//  ^0  ) 


For  Si-  or  SiGe-Al  interfaces,  Schottky  barrier  height 
is  roughly  of  the  semiconductor  bandgap.  We  can 
further  simplify  the  expression  of  as 


^iSCH  ~ 


^Qsch  j:2 


4  exp 


^SCH 


4  J 


c. 


where  the  correction  factor  represents  barrier  heights 
and  bandgap  parameters: 


C  = - ^ — 

corr  2 

SCH 


9P. 


^  f -3P  ^ 

exp 

V  V 


4^0 


■"corr  Q  p  Z 

^^SCH 


44 


_44 
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Since  {^Psch  » ■^).  C^on-  is  further  approximated,  and  the 
current  density  expression  is 


4< 


f^QscnP'g  „2  I 

-6P—^  “Pi 


SCH 


p 

pSCH 


Derivation  of  Parallel  Electric  Field 

The  derivative  of  potential  distribution  with  respect  to 
the  position  in  the  tunneling  direction,  yields  the  parallel 
electric  field  expression: 

= - 7 - \  {(‘^/  ~  )  exp(-  ))  exp(-  a/^x) 

1  —  expl  —  2r y  I 

-((<r/  +4')-o'/  exp(-r^))exp(-«^'"(Z.-x))| 


where 


^  ft  f  t 

2  ^  ^  ox  *  xi  ^  ox*'  box 

V  ^si^ox  ^ox^ox  J 


1  +  2 


T^=L-a/'^ 


At  the  metal-semiconductor  interface  (x  =  0 ), 


e  1/ 

4/  “  «/ 


1  -  expl 

(-r/] 

1  .  exp( 

07/  \ 

-r/] 

1  1 
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Theory  of  Generalized  Drift-Diffusion  Model 
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Abstract 

By  taking  the  appropriate  moments  of  the  Boltzmann  transport  equation,  we  show  that  the 
mobility  and  diffusion  constants  can  be  rigorously  defined  for  the  most  general  nonstationary 
transport  condition,  including  the  velocity  overshoot  and  the  ballistic  case.  The  moments  are 
obtained  by  using  a  special  parameter  called  Asymmetric  Relaxation  Time.  The  theory  can 
provide  a  fundamental  basis  for  constructing  compact  circuit  models  for  all  deep  submicron 
devices  and  sub  lOOnm  devices.  It  can  also  be  used  to  construct  efficient  models  for  device 
TCAD  simulators. 


I.  Introduction 


1 


The  Boltzmann  Transport  Equation  (BTE)  provides  the  foundation  of  modeling  carrier 
transport  in  modem  devices  and  the  devices  of  near  future.  The  conventional  drifl-difiusion 
(DD)model  for  carrier  transport  is  the  first  moment  of  the  BTE.  [1]  As  the  minimum  feature  sizes 
of  the  devices  came  to  the  submicron  regime,  it  became  necessary  to  include  higher  order  of 
moments  of  the  BTE  to  account  for  the  nonstationaiy  effects.  The  resultant  set  of  equations  are 
usually  referred  to  as  the  hydrodynamic  equations  (HD)M  [2].  Today,  most  of  the  existing  device 
simulators  incorporate  some  extent  of  a  HD  model.  The  differences  among  the  various  HD 
models  are  due  to  different  ways  of  modeling  the  carrier  distribution  function,  the  carrier  energy 
band  stracture,  and  the  microscopic  relaxation  times  of  the  collision  term  [2].  Of  course,  the 
Monte  Carlo  (MC)  method  is  the  most  accurate  solution  of  the  BTE  [3].  It  is  still  considered  to 
be  too  uneconomical  to  be  directly  used  for  semiconductor  device  design. 

Even  though  the  HD  model  is  more  efficient  than  the  MC  method,  it  is  still  much  more 
time  consuming  than  the  DD  model.  As  a  engineering  design  tool,  both  the  accuracy  and 
economy  of  the  model  are  extremely  important.  To  meet  the  trend  of  a  shortened  development 
cycle  and  more  stringent  performance  demand,  the  semiconductor  industry  has  an  immediate 
need  in  accurate  physics  based  device  simulators  with  fast  computation  tum-around  time.  There 
are  many  reasons  that  a  new  approach  should  be  sought.  We  have  endeavored  an  investigation  in 
this  direction  and  found  that  a  generalized  drift-diffusion  (ODD)  model  can  be  developed 
rigorously  based  on  BTE.  The  ODD  approach  can  be  very  accurate  since  the  two  basic 
parameters,  generalized  mobility  and  diffusivity,  are  precisely  defined.  We  then  need  to  model 
these  parameters.  Once  they  are  properly  modeled,  the  simulation  time  will  be  similar  to  the 
conventional  DD  model  which  is  known  to  be  very  efficient. 


II.  Theory 

The  Boltzmann  Transport  Equation  is  given  by 


dt 


n  dt 


(1) 


2 


^  -f 

For  the  analysis  of  very  small  devices,  (-^)  s  0 . 

ot 


The  distribution  ftinction  is  determined  by  three 


terms.  The  term  with  corresponds  to  diffiision  motion,  the  term  with  corresponds  to  drift 


motion  and  the  last  term  represents  the  effect  of  collisions.  For  non-degenerate  carriers. 


V  ^  /coll 


k  ,r,t)-S(k  ,k) 


f(k,r,t)-S(k,k') 


(2) 


To  make  the  most  general  interpretation  of  drift  and  diffusion,  we  will  define  a  relaxation 
time  called  Asymmetric  Relaxation  Time  (ART).  The  ART  will  be  the  momentum  relaxation 
time  (MRT)  when  the  latter  exist.  Note  that  the  existence  of  momentum  relaxation  time  gives  a 
most  simple  but  rigorous  description  of  the  drift-diffusion  model.  When  scattering  mechanisms 
are  velocity  randomizing  or  quasi-elastic,  the  momentum  relaxation  time  exists  [4].  In  silicon 
these  conditions  are  satisfied  in  an  approximate  sense.  But  for  scattering  mechanisms  such  as  the 
polar  optical  phonon  scattering  in  GaAs,  the  conditions  are  not  satisfied.  The  ART  can  play  the 
same  role  as  MRT. 

From  the  BTE,  it  is  natural  to  divide  the  current  density  into  two  parts  even  in  the  most 
general  case.  As  shown,  the  time  independent  BTE  has  three  terms:  diffusion,  drift  and  collision. 
If  we  can  take  an  appropriate  moment  of  the  collision  term  to  make  it  a  current  density,  we  then 
can  obtain  a  general  drift  and  diffusion  model.  Of  course,  in  order  to  be  useful,  the  chosen  form 
of  the  moment  should  be  simple  and  have  clear  physical  meaning.  It  turns  out  such  a  form  is 
indeed  possible.  We  find  that  we  can  divide  the  distribution  function  into  symmetric  and 
asymmetric  parts.  Due  to  time  reversal  symmetry,  the  collision  matrix  elements  have  some  very 
general  properties.  It  can  be  shown  that  these  properties  will  maintain  the  symmetry  and 
asymmetry  of  ftie  distribution  function.  Only  the  asymmetric  part  of  the  distribution  provides  the 
current.  We  can  use  the  relaxation  time  ART  for  the  asymmetric  distribution  to  define  the 
required  form  of  the  moment.  As  will  be  shown  in  the  following,  the  symmetric  part  of  the 
distribution  makes  zero  contribution  to  the  moment  due  to  time  reversal  symmetry.  The  collision 
term  then  produces  the  current  density  and  the  other  two  terms  generate  the  precise  definitions  of 
the  generalized  drift  and  diffusion  current.  The  ART  is  similar  to  the  MRT  in  physical  meaning 
but  is  much  more  general.  It  encompasses  all  the  transport  conditions  that  BTE  can  describe. 
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including  the  nonstationary  transport  and  nearly  ballistic  case.  A  more  detailed  account  is 
provided  in  the  following. 

To  define  ART,  separate /into  symmetric  and  asymmetric  parts 


f{r,k,t)  ■■ 

=  foirJd)  +  Mr,k,t) 

(3a) 

f{r,k,t)  +  f(r-k,t) 

2 

(3b) 

fxir,k,t) 

f(f,k,t)-f(r-k,t) 

2 

(3c) 

Therefore, 

II 

1 

(%)  +  (Ml) 

V  0^  Jcoll  ^  ^  Qj^ 

(4a) 

Define  (^ 
dt 

s  _  -fxir,k,t) 

T(r,k,t)  ’ 

(4b) 

where  the  relaxation  time  T(f,k,t)  is  called  ART  because  it  describes  the  collision  effects  on  the 

asymmetric  part  of  the  distribution  function.  It  can  be  shown  that  this  relaxation  time  alone 
determines  the  local  current  density.  Take  the  momentum  of 

Q(r,k,t)  =  -  ev(k)r{r,k,t)  .  (5) 

Then,  we  have 


-  «v-V,/)  -  =  0 


h  *"  “'a 

dk 


dt 


=  J{r,t) 


{Qhf-Xou)  =  -  e 

ot 


It  can  be  shown  in  Appendix  A  that 


dk 


dt 


'it 


(6a) 

(6b) 

(6c) 
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T{r,k,t)  =  T{r-k,t), 

o(^.^’0!L  _  1^0  (^-^>0 


CO/? 


4 


rv  n  ^ 

From  Eq.(7a)  and  (7b),  it  can  be  shown  that 

=  0. 

ot 


(7a) 

(7b) 


(8) 


We  can  neglect  the  r  dependence  of  T(r,k,t)  when  compared  with  f(r,k,t)  (which  changes 
much  faster  than  T).Then,  the  diffusion  term  can  be  written  as 

dk  „  _ r  ..  dk 


-e7{kyt{rjMv-^Jir,k,ty)^=  -  e  V-2^/(r,^,/) 


Att 


(9) 


We  assume  that  /  changes  in  the  direction  of  the  electric  field,  the  diffusion  current  deisity  can  be 
defined  as  (if  not,  it  can  be  easily  generalized  as  a  tensor) 


'^diffusion  e  Vj:D(r,t)-n(f,t) 


(10a) 


Then,  the  generalized  diffusivity  is  given  by, 


k,  tyifir,  k,  t)  ^ypfo  (r,  k,  t) 


D(r,t)  s 


4;r  _ 


dk 

At^ 


-  7  dk 


^r,k,t) 


^(r,k,t) 


dk 


(10b) 


4;r' 

Similarly,  drift  term  can  be  used  to  define  a  generalized  mobility 

-  c  7(^)r(r,^,/)^-V.-/~(F,^,/)|^s  +  en^(r,E,t)E(r,t) 


(11a) 


with 


ju(f,t)  = 


dk 

An^ 


dk 

4;r^ 


(11b) 


The  generalized  diffusivity  and  mobility  defined  in  Eqs.(lOb)  and  (lib)  applies  to  any  situation 
of  nonstationary  transport,  including  the  limiting  case  of  nearly  ballistic  transport. 
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III.  Simple  Cases 

Some  simplified  expression  can  be  derived  if  the  distribution  function  can  be  simplified.  In  the 
following,  it  is  shown  for  the  case  that  a  single  hot  electron  temperature  can  be  used  to 
characterize  the  distribution  function. 

When  / (r,k,t)  =  ,  it  can  be  shown  that 


juir,t)  =  -  =  —D(r,t), 

—  dk  kZ 


Z 


(12) 


4;r^ 


where  7^  is  the  hot  electron  temperature. 

For  a  more  general  case  that  a  single  hot  electron  temperature  can  not  be  used  to  characterize  the 
distribution  but  the  asymmetric  part  of  the  distribution  can  be  neglected  in  the  generalized 
mobility,  then 

dk 


Krd)  = 


k,t)E-  V^-/o  {r,  k,  t) 


Att 


dk 


^  E-Vj-fo(r,k,t)  = 

uJL 


(13) 


fo(r,  k, 


M(rd)  = 


dE 


dk 


The  generalized  Einstein  relationship  becomes, 

=  —D 

kZ 

The  carrier  temperature  for  Einstein  relation  is  given  by. 


kZ 


...  ....  -  -'4^^ 


dk 

Ur,k,t)^ 

An 


(14) 


(15) 


(16) 
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If  the  band  is  isotropic, 

TVp  =  —  Tv^  and 
£  3 


^  ~  '^f,D{E)dE 


(17) 


In  Eq.(17),  the  hot  carrier  temperature  for  generalized  Einstein  relation  is  given  by  the  local 

d  y 

energy  temperature - In weighted  by  the  product  t(E)v(E)  fQ(E)D(E). 


IV.  Discussion 

There  are  many  good  reasons  why  this  approach  should  be  very  useful. 

1  The  TCAD  field  is  becoming  a  mature  industry.  The  time  to  industry  availability  of  new 
physical  models  is  shrinking  from  years  to  days.  This  important  trend  encourages  the 
development  of  modules  of  device-specific  models.  An  outstanding  example  is  the  ultra-fast 
device  simulation  with  MC  tuned  transport  models  in  FastBlaze.  The  speed  is  greatly 
enhanced  while  accuracy  is  the  same  or  even  improved.  With  the  GDD  model,  the  speed  will 
be  further  improved  when  compared  with  the  energy  balance  simulation.  The  accuracy  can 
also  be  improved  because  the  GDD  parameters  are  precisely  defined  while  the  energy 
balance  parameters  are  approximately  defined.  Of  course,  the  GDD  parameters  can  be  MC 
tuned. 

1.  As  the  minimum  feature  sizes  of  the  devices  are  moving  into  the  sub-lOOnm  regime,  nearly 
ballistic  types  of  motion  may  appear  in  a  substantial  part  of  the  active  device  region. 
Therefore,  more  and  more  higher  order  moments  of  BTE  are  required  to  ensure  accuracy. 
This  can  be  easily  expected  since  a  nearly  ballistic  distribution  function  requires  many 
moments  to  represent  it.  The  required  simulation  time  will  then  increase  greatly.  As  will  be 
shown  later,  the  GDD  parameters  are  precisely  defined  in  terms  of  the  local  distribution 
function.  It  takes  care  of  nearly  ballistic  electrons  as  easy  (  in  some  sense  easier)  as  the 
thermalized  hot  electrons. 
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2.  As  more  hetorojunctions  are  used  in  special  devices,  the  tunneling  current  may  affect  the 
mobility  and  diffiisivity  in  the  future  devices.  The  HD  model  will  be  again  very  tedious  in 
dealing  with  these  effects.  For  the  GDD  model,  it  is  natural  to  treat  these  effects  as  in  the 
nearly  ballistic  type  of  transport. 

4  As  the  Semiconductors  Roadmap  [6]  has  emphasized,  timely  access  to  accurate  computer 
models  of  each  technology  area  will  be  necessary  to  managing  the  growing  complexity  of  the 
semiconductor  technology.  Different  levels  of  abstraction  for  models  are  essential  to 
effectively  shorten  time  scales,  lower  cost,  and  increase  quality.  In  device  TCAD  area,  how  to 
build  compact  circuit  models  for  device  is  as  important  as  the  device  simulation.  In  doing 
this,  effective  parameter  extraction  for  device  simulation  results  are  of  increasing  importance. 
Since  most  analytical  works  for  device  modeling  have  been  done  in  terms  of  drift-diffusion 
model,  the  usage  of  GDD  will  naturally  make  the  task  of  compact  circuit  modeling  much 
easier  than  the  HD  or  MC  approach.  For  example,  there  is  up  until  now  no  rigorous  Gummel- 
Poon  type  model  for  HBT  since  the  nonstationary  transport  involved  here  cannot  be  put  into 
simple  analytical  forms.  Using  GDD  approach,  it  will  be  straightforward  to  set  up  such 
compact  circuit  models. 

This  work  is  supported  by  UCLA  JSEP  program. 


Appendix  A:  Proof 


Note:  in  the  following ,  we  present  the  proof  of 


(1)  T{f,k,t)  =  r{f-k,t) 


To  prove  (1),  let 


1 

T(r,kJ) 
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=  l^sikS) - ^-5 - — 

-  ‘ 

T(r-k,t) 

By  time  reversal  symmetry  [5] 

S(k,k')  =  S(-k-k') 

Si-k,k')  =  S(k-k') 


Z^n[(r,f,t)-/(n-k',OS(k',k) 

^oTT 

f(F,-k,0-f(r,k,t) 


- ^  =  "^Sik.-k)  -  -  ‘  ^  '^Xi(f,k\t)-f(.rA\t)S(.-k\k) 

i(r,-k,t)  ‘Sir’  f(r,-k,t)-f(r,k,t)  ‘Sk' 

=  ’^jSikX) - -  ^  ''i^U(f-k\t)-f(r,k\t)Sii\k)=-^^ 

^7^  f{r,k,t)-f{r-k,t)  r{r,k,t) 


This  is  reasonable,  because  only  when  r{r,k,t)  =  r{f-k,t) ,  the  asymmetry  of  f^(r,k,t)  can  be 
maintained. 

To  prove  (2) 


na  Tv 

cou 


(r,k',t)-S(k',k)  -  fQ(r,k,t)-S(k,k') 


ir,k\t)-Sik'-k)  -  f,(r-k,t)-Si-k,k') 


S(k' -k)  =  S(-t,k)  (time reversal) 


(r,k',t)-Si-k',k)  -  f,(r-k,t)-Sik-£) 
(r,k",t)-S(k\k)  -  f,ir,k,t)-Sik,£') 


W/ 
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150GsampIe/s  wavelength  division  sampler 
v\^h  time-stretched  output 

A.S.  Bhushan,  F.  Coppinger,  B.  Jalali.  S.  Waag  and 

H.F,  FetteiTnan 


A  new  teduiiquc  for  optical  sampling  and  time  stretdiiiig  of  an 
arifJogiie  signal  is  planted.  The  signal  is  ‘wavelength  stamped’ 
by  a  chirped  optical  pulse.  An  annyed  waveguide  grating  provides 
parallel  thne-iT^terleax'ed  samples  us  well  ns  a  tiinc-strctchcd  serial 
sanipkxl  waveform. 

High-speed  a-naloguo-to-digitid  conv'crsion  systems  require  sam¬ 
pling  techniques  that  have  high  temporal  resolution  and  fast  sam¬ 
pling  rates.  Also  of  great  impoitance  is  a  mechanism  for 
demultiplexing  or  slowing  the  .samples  in  tijne  as  they  approach 
the  electronic  quantiser.  If  time  could  be  mapped  onto  wave¬ 
length,  then  the  problem  of  demullipkxing  in  time  could  be  trans¬ 
lated  into  a  problan  of  demultiplexing  in  wavelength.  The  latter 
can  be  resolved  easily  by  using  just  passive  optical  filters.  Tinic-to- 
wavelength  mapping  by  mixing  a  chirped  pulse  with  a  digital  data 
stream  has  been  used  to  demullipk.x  digital  serial  data  [1].  Chirp¬ 
ing  a  broadband  transform-limited  pulse  separates  different  spec¬ 
tral  components  of  the  pulse  in  time.  When  the  chirped  pulse 
enters  the  mixer,  diirerent  wavelength  components  arrive  sequen- 
licdly.  Thus,  time  Is  mapj^  into  the  optical  wavelength.  Tills  tech¬ 
nique  has  also  been  applied  to  transmission  of  WDM  data  [2]  and 
to  sample  analogue  signals  [3,  4],  In  this  Letter,  we  demonstrate  a 
novel  sampling  ardiitecture  ba^  on  two  comf^imentary’  con- 
cepb’.  The  first  is  the  diirped-pulse  timc-to-w'avelength  transfor¬ 
mation.  The  seamd  is  the  use  of  discrete  dispersion,  obtainable 
with  an  arrayed  waveguide  grating  (AWG)  [5],  to  slow  down  the 
j^ampled  waveform. 


Kig.  1  Wavelength  division  .sampiing 

Different  shades  of  grey  represent  difTereiit  wavelengths 


Tb&  cliirpcd  puLse  sampling  is  shown  in  Fig.  1.  To  generate 
chirped  pulses,  vve  use  a  passively  modelocked  erbium-doped  fihm 


persioii  m  a  singlcmode  fibre  (SMF,  D  -  1 7pN^mn.km).  Sampling 
is  performed  by  a  silica  .AW<3  with  a  channel  spacing  of  0,8  nm 
and  a  pa.ssband  of  0.3 nm.  Ihc  clcctriail  signal  and  the  chirped 
pulse  interact  in  the  optical  modulator  to  produce  a  'wavelength 
stamped*  intensity-modulated  optical  signal.  Tlic  AWG  samples 
the  optical  spectrum  and  thus  samples  the  dmc  evolution  of  tlie 
electrical  trajxsienl.  In  Fig.  2,  w^e  report  the  timing  measuremenls 
for  a  chirp  of  I20ps/hm  (7km  of  SMF).  Outputs  from  individual 
AW'G  ports  are  phol.odctectcd  and  displayed  using  a  digital-sam¬ 
pling  oscilloscope.  The  puLse-lo-pulse  spacing  of  l.OOps  in  Fig.  2 
corresponds  well  to  the  channel  spacing  of  O.Snm  for  the  AWG. 


Haeh  trace  is  a  dillerent  outpiit  of  the  .AWG 


FSL  denotes  femtosecond  laser 


The  setup  for  the  1  oOCjsamplc's  experiment  is  showr  in  Fig.  3. 
Vhe  transform-limited  pulse  jQrom  the  EDFL  was  dispersed 
through  a  Coming  SMF*  lo  produce  a  chirped  pulse.  A  portion  of 
die  temto-second  laser  output  is  fed  to  a  photodelecior  to  provide 
the  syttchronising  clock  for  the  entire  sysicuL  The  diirped  pulse  is 
modulated  and  then  fed  to  an  8-channcl  silica  AWG.  A  Smn 
bandi^ss  filter  piccedes  the  AWG  so  that  dK  serialised  output 
contains  only  one  free  spectral  range  of  the  A  WG.  The  outputs  of 
ihe  AWG  arc  fed  back  to  the  ^responding  inputs  witli  inci*emen- 
tal  delays  of  around  L3ns.  Tliis  feedback  architecture  provides 
discrete  dispersion  and  separates  samples  in  time.  In  an  analoguc- 
to-digilal  (A/D)  converter  sy'stem.  a  sini^c  slow  electronic  quan¬ 
tiser  can  ihcn  be  used  to  digitise  a  fast  dcctrical  tiaasient.  The 
feedback  paths  include  variable  attenuators  to  provide  spectral 
equalisation.  Spectral  equalisation  con  also  be  attained  by  using  a 
dispersion-decreasing  super-a>nlinuum  fibre  [6]  or  by  a  fe«i  for¬ 
ward  archilecturc  [2J.  The  sampling  rate  is  given  by  £  =  l,^.  D, 
where  A/c  is  the  AWG  ehannel  spacing  and  O  is  the  total  disper¬ 
sion.  Wc  use  a  SOOm  spool  of  fibre  for  disjiersion  corresponding 
to  a  chirp  rate  of  8.5ps/hui  at  the  modulator  input.  The  pulse  sep¬ 
aration  (in  lime)  at  the  output  of  each  A.WG  channel  is  6.8ps,  cor- 
re.sponding  to  a  sampling  rale  of  147GsampIe^s.  The  response  time 
of  the  photodetector  {30 ps)  and  the  jitter  noise  in  the  measuring 
electronics  (Ups)  make  measuring  the  actual  pulse-to-puLse  sepa¬ 
ration  difficult  However,  the  separation  between  the  first  luid  the 
last  channel  oulpuls  (before  feedback)  was  mcasui^  to  be 
wliich  corresponds  to  an  average  pulse  separation  of  6.4ps.  Fig.  4 
shows  Ihc  time-stretched  serialised  output.  The  output  rate  is 
reduced  from  147GH7  to  770MFlz  by  i.3ns  incremeniai  delays  in 
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sample  stream  is  l4dB.  In  general,  the  SNR  has  a  contribution 
froTTi  the  spectral  nonuniformity  (from  laser  and  AWG)  as  well  as 
crosstalk  within  the  AWG.  Since  the  measured  AWG  opLiciii 
crosstalk  is  27dB,  we  conclude  that  the  measured  SNR  is  limited 
by  the  incomplete  spectral  equalisation,  w^hich  can  be  improved 
with  higher  quality  attenuators. 


0  1  2  3  4  3  3 

time, ns 


Fig.  4  Time  stretched  and  serialised  output  for  I47Gsample/s  sampiirtg 

Modulator  voltage  -  IV  (top  trace)  and  8V  (bottom  trace) 

^"rt=10V 

(i)  X  -  I565.4nm 

(ii)  X  “  1 566.2 nm 

(iii)  X  =  1567wn 

(iv)  X  =  1 567.8 mn 

(v)  X  =  1568.6iim 

Tn  suntmary,  we  have  deanonslrated  a  novel  optical  sampling 
method.  The  chirped  pulse  technique  provides  time-interleaved 
samples  compatible  with  interleaved  A/D  architecture.  Using  an 
AWG  in  a  recirculating  coniiguration,  a  stream  of  samples  can  be 
stretchal  in  time  and  serialised,  allowing  an  ultra-fast  cleclxical 
transient  to  be  digitised  with  a  single  slow  AD. 

AcknowlediiimenLs:  This  work  is  supported  by  JSEP  (Joint  Service 
Electronic  Program)  and  the  CNR  (Office  of  Naval  Reseamh) 
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S.  Yamamoto  and  S.  Akiba 

i60Gbit/s  (8  X  20  GbiVs)  sollton  WDM  signals  were  successfully 
tiansmittcd  over  4000kiii  with  a  BF.R  of  <  10^  by  using 
periodically  dispersion-compensated  dispcisioa-flattencd  tibre 
'  with  an  average  dispersion  slope  of  0.0005  piwOoii/iiiir. 

Soliton-WDM  transmissioii  is  attractive  for  large-capacity  long- 
haul  optical  transmission  systems  because  soliton-based  systems 
liavc  the  potential  to  carry  higher  bit  rate  signals  than  NRZ  .sys¬ 
tems.  The  key  technological  issues  in  soliton  WDM  transmission 
arc  the  reduction  of  collision-induced  timing  jitter,  four-wave  mix¬ 
ing  [1]  and  the  dispersion  slope  of  the  transmission  fibre  [2].  The 
first  tw'o  requirements  cm  be  achieved  using  dispersion  tapeied 
fibre  spans  with  soliton  control  techniques  [3,  4].  To  overcome  the 
impact  of  the  dispersion  slope,  the  use  of  dispersion-flattened  fibre 
seems  quite  attractive  and  practical  compared  with  a  cliannel-by- 
channel  dispersion  compemjation  scheme  [2,  4].  So  fai\  we  have 
made  a  preliminary  report  on  the  effectiveness  of  the  soliton 
WDM  systems  using  dispersion-flattened  transmission  fibre  and 
80Gbil/s  (4  X  20Gbil/s)  signals  [5].  Tn  thLs  Letter,  we  show  the 
experimental  results  on  160Gbil/s  (8  x  20  Gbit/s)  soliton  WDM 
traiismission  over  4000kin  using  dispersion-flattened  fibre  with 
periodic  dispersion  compensation  with  a  w^avelenglli  bandwidtlt  of 
11.2iun. 


Fig.  1  Jsxperimenral  setup  of  I60GhiL<s  {8  x  20  Gbit/s)  soiiton  WDXf 
transmission 
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Fig.  1  show^  a  schematic  diagram  of  the  experimental  setup  for 
dghl-channel  20GbiL/s  soliton  WDM  transmission.  We  used  four 
2WDMx20  Gbitf's  transmitters  as  a  IbOGbit^  transmitter.  The 
w'avdengths  of  eight  diannels,  ranging  from  1546  to  1558 nm,  were 
equally  spaced  by  1.6nm.  Each  2WDM  transmitter  consists  of  two 
EA-modulalor-based  soliton  pulse  generators  with  dilTerent  wave¬ 
lengths.  A  2DGbiys  optical  soliton  data  stream  was  produced  by 
using  two  LiNbOj  intensity'  modulators  operated  at  lOGbit^s  by  a 
2’^  -  I  pseudorandom  binary  sequence  and  optical  time  division 
multi{^cxing.  The  pulse  tvidth  obtained  was  from  10  to  15ps,  In 
addition,  each  20Gbit/s  signal  was  modulated  at  20GHz  by  an 
IJNbOj  phase  modulator  to  improve  the  ti'ansmissioii  perform¬ 
ance  [6,  7].  Afl^*  combming  ei^t  20  Gbit/s  signals  so  that  the 
adjacent  channels  wei'e  orthogonally  polarised,  a  polarisation 
scrambler  w'as  used  to  suppress  the  PHB  of  the  EDFA  repeaters. 
The  transmission  line  comprises  three  spans  of  ^27iaii  long  disper- 
sion-flaiiened  fibre  (DFI^  and  dispersion-compensation  fibre 
(DCF).  The  loss  of  the  DFF  spans,  including  the  splicing  loss,  was 
still  large  and  the  average  span  lass  w'as  ^dB.  The  average  disper¬ 
sion  slope  of  the  DFF  was  0.0005ps/km/nm^  and  the  average  chro¬ 
matic  disperaon  of  the  DFF  at  1555nm  was  l.lps/km/nm.  The 
accumulated  cliromatic  dispersion  of  the  k>op  was  compensated 
for  by  the  DCF  with  82p^nm  and  the  residual  dispenion  in  the 
loop  was  0.02ps/km/nin.  As  a  frequency  guiding  filter,  an  FP- 
etalon  filter  {FSR  -  1.6nm)  was  placed  in  the  loop  to  stabilise  tlie 
transmission  pcrfonnancc.  In  the  receiver,  tfie  desired  cliannel  was 
selected  by  optical  bandpass  filters  and  the  transmitted  20 Gbit/s 
signals  were  optically  time-dhasion-demuluplcxcd  to  a  lOGbib's 
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chirped  optical  pulses 


gauon  mrougii  me  secona  aispersive  cicmcni.  me  iiign;u  ij> 
slreiched  in  time  by  a  factor  of  3.25  with  high  fidelity. 


F.  Coppinger,  A.S.  Bhushaa  and  B.  Jalali 

Ihc  authors  report  a  method  for  stretching  electrical  signals  in 
time.  A  liigh  chirp  laie  is  imposed  on  tlie  electrical  signal  by 
mixing  it  with  a  dispersed  u!ti*a-short  optical  pulse  in  an  electro* 
optic  intensity  modulator.  This  is  followed  by  a  passive  optical 
dispersion  element  to  produce  a  tiinc^nagnified  copy  of  the  input 
cicctricai  signal 


Time  manipulation  of  signals  has  been  proposed  as  a  possible 
technique  ro  match  the  data  rates  of  signals  to  the  receiver  or  to 
observed  very'  fast  phenomena.  Considering  the  limc'Spacc  equiva¬ 
lence  of  dispersion  and  diffraction,  and  in  direct  analogy  to  a  spa¬ 
tial  lens,  a  signal  nuiy  be  stretched  in  time  by  subjecting  il  to  a 
dispersion-quadratic  phase  modulation-dispersion  sequence  [1,  2], 
While  this  concept  has  been  known  for  some  lime,  few  sucocsslul 
experiments  have  been  carried  out  due  to  tlie  difficulty  of  obtain¬ 
ing  high  quadratic  pliase  modulation  rate  and'or  large  bandwidth 
dispersive  elements.  In  his  pioneering  work,  C'aputi  demonstrated 
an  all-elccirical  time  simtch  system  [1].  In  the  electrical  domain, 
liighly  dispersive  elements  are  available,  however  the  small  band¬ 
width  of  these  elements  limits  the  capability  and  performance  of 
an  all-elcctrical  system.  Large  bandwidth  dispersive  elements  are 
easily  obtained  optically,  but  the  dispersion  remains  low  requiring 
a  high  quadratic  phase  modulation  rale  (or  linear  frequency'  chirp¬ 
ing).  Resonant  optical  phase  modulators  were  pi'oposcd  as  a  way 
to  attain  a  high  modulation  rale  [3].  While  high  chirp  rates  are 
iichicvablc  with  this  device,  the  duration  of  the  chirp  is  limited. 
Pulses  from  a  Nd:YAG  laser  were  also  used  to  create  Llic  chiip  in 
an  all-optical  time  lens  12].  The  time  aperture  of  the  lens  was  then 
limited  by  the  bandwidth  of  the  pulse  (620GHz).  In  this  Letter,  wc 
demonstrate  an  optoelectronic  time  magnification  system  with  an 
electrical  input/dulput  and  with  the  linear  chirp  provided  by  a  dis¬ 
persed  optical  puLse  from  a  modclockcd  (^bium-doped  fibre  laser 
(EDFL).  In  addition  to  the  ultra-high  bandwidth  ('*-7.5THz)  tmd 
chirp  rate,  this  system  differs  from  previously  reported  time  mag¬ 
nification  systems  in  that  the  bandwidth  of  ihe  input  signal  is  neg¬ 
ligible  compar*cd  to  the  chirp  bandwidth.  With  this  property,  an 
approximation  to  an  ideal  time  lens  can  be  obtained  without  the 
need  to  disperse  the  input  signal  This  results  in  a  more  relaxed 
and  simpler  design  of  a  time  stretching  system  that  can  readily  be 
implemented  with  commercially  available  components.  Tliis  tedi- 
nique  is  promising  for  analogue-to-digital  conversion  of  ultra-fast 
electrical  signals. 


EDFL 


input  signal 

J 

photodetector 


Fig.  1  Experimental  for  lime  magnification  using  chirped  optical 
pulse 


Fig.  1  describes  our  implementation  of  the  optoelectronic  time 
magnification  system.  The  160  fs  pulse  generated  by  the  EDFL  is 
dispersed  in  the  input  fibre.  When  the  pulse  reaches  the  modula¬ 
tor,  different  frequency  components  arrivic  at  difT^rent  times.  The 
electrical  signal  to  be  stretch^  is  efficiently  mixed  with  the  diirped 
optical  signal  by  the  Mach-Zehnder  ekctro-oplic  modulator 
(12GH2  bandwidth).  The  efficiency-  of  the  electro-optic  mixing  is 
significantly  higher  than  opdcal-optical  mixing  [41  in  a  nonfinear 
crystal  Tlie  resulting  intensity'-modulalcd  chir^  signal  is  dis¬ 
persed  in  a  second  spool  of  fibre  and  its  envelope  is  detected  by  a 
fast  pholodelector  with  a  30ps  response  time. 

In  the  experiment,  wc  use  a  2km  spool  of  fibre  to  generate  the 
fixx|ucncy  ramp.  Wlh  a  dispersion  of  ITps/nm^km  in  ic  fibre,  the 
diup  rate  at  tlic  modulator  is  a.SGHz^'ps.  To  generate  an  arbitrary 
transient,  vto  bias  ihe  modulator  at  and  apply  an  electrical 
pulse  to  il.  Fig.  2a  shows  the  resulting  waveform  measured  at  the 
modulator  output.  Tt  represents  tlie  applietl  electrical  pulse,  its  sec¬ 
ond  harmonic  originating  from  biasing  the  modulator  at  The 
envelope  is  also  somewhat  shaped  by  the  spectral  envelope  of  the 
chirped  pulse.  The  second  dispersion  consists  of  a  5.5  km  spool  of 


Fig,  1  Signal  envelope  after  modularion  and  niagnifed  signal  afier  sec¬ 
ond  dispersion 


a  Signal  envelope  after  modulation 
b  Magnified  sigmil  after  second  dispersion 


The  magniilctitioQ  factor  in  our  sy'stan  depends  on  both  the 
iunount  of  dispersion  in  Ihe  second  dispersive  stage  D2  and  the 
chirp  rate.  The  magnification  factor  can  he  increased  by  increasing 
D2  or  by  raising  the  chirp  rate  (decreasing  A)  However,  increas¬ 
ing  the  diirp  rate  also  reduces  the  time  aperture  of  the  sy'Slcm, 

The  fundamental  resolution  of  a  conventional  fime  lens  in 
which  the  image  is  at  the  focal  time  of  the  lens  is  limited  by  its 
bandwidth  and  is  estimated  to  be  where  B  is  ihe  chirp  band¬ 
width  [5].  If  the  image  is  not  at  the  focal  time  (the  case  in  our 
experiment),  the  image  will  be  distorted.  However,  the  resulting 
distortion  is  negligible  as  the  electrical  signal  bandwidth  ('-12GHz) 
is  negligible  corapaicd  to  the  bandwidth  of  ihc  chiip  (--T.STHz). 
Further,  the  distortion  is  averaged  out  by  the  finite  response  dine 
of  the  photodetector  as  long  as  a/;A  <  V» 
width  of  the  decirical  signal  and  is  the  response  time  of  the 
pholodelector  (30  ps).  The  resolution  and  linearity  of  the  lens  are 
also  limited  by  the  non-constant  dispersion  inside  the  fibre  over 
tlic  spectrum  of  the  EDFL.  DLspersion-llaitened  fibre  would  miti¬ 
gate  this  problem.  The  influence  of  the  non-flat  optical  spectrum 
on  the  input  signal  can  be  minimised  by  spcctnim  equalisation 
using  a  fibre  Bragg  grating  filter,  disperskm  deaeasing  supercou- 
tinuum  fibre  [6],  or  a  feedtbnv'ard  architecture  [7], 

In  summary,  w-e  have  used  the  high  chirp  rale  offered  by  a  dis¬ 
persed  pulse  from  an  EDFL  to  stretch  an  cleclrical  signal  in  time. 
The  frequency  ramp  is  inleiisity  modulated  by  the  electrical  signal 
in  a  Mach-Zehntte  doclro-optic  nuxlulator.  Tlie  chirped  signal  is 
that  further  dispersed  leading  to  time  mt^ification.  This  tech¬ 
nique  has  potential  for  the  digilisaiion  of  high  frequency  electrical 
signals. 
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uimeo  on  or  on.  i  no  arrays  were  laDiieatea  on  a  iiommaiiy  ^-ivpc 
Ga.4s  epitaxial  layer,  with  a  carrier  cx)nceiitration  of  <  IO^^ctyj-.  i 
Dielectric  isolation,  layers  were  omitted,  but  the  dark  current  was  | 
<  lOOnA  at  an  operating  bias  (0-10  V).  Tlie  interdigitated  elec-  | 
nodes  have  finger  spacing  and  a  width  of  2pitL.  After  dicing,  the  ‘ 
chips  were  mounted  on  dielectric  carriers  with  RF,  coplimur-  i 
waveguide  aixess  lines.  The  bias  pad  of  each  MSM  was  wire  • 
bonded  to  an  AC  ground  via  a  bia.s  decoupling  capacitor.  The 
common  bus  was  connected  to  the  coplanar  lines  on  die  chip  car-  I 
Her,  via  multiple  wire  bonds  at  eadi  end.  | 


photodeiector  1  photodetector  2 


Fig.  1  Simple  equivalent  circuit  for  curay  described  in  text 
VACh  cud  of  array  is  terminated  in  system  impedance 


Packaged  array  of  eight  MSM 
photodetectors  with  uniform  12GHz 
bandwidth 


R.G.  DeCorby,  R.I.  MacDonald,  A.J.P.  Hnaliw, 

D.  Boertjes,  J.N.  McMullin,  F.  Gouin  and  J.  Noad 

Tran.smi<;sio[i  liiK  concepts  wiere  used  in  the  design  of  a 
broad  baiKl,  intcgi‘atcd  pliotodctector  array.  The  array  consists  of 
eighl  meUil  semiconductur-metiU  (MSM)  photodetectors,  with 
aciive  regions  of  (26pm)^,  disliibuled  along  a  common  bus.  A 
uniform  3dB  bandwidth  of  12GH/.  is  clernonsiTHie^i.  Tlnese  arrays 
are  central  to  opiockctronic  switching  tedmiques,  and  are 
applicable  to  luicrow'ave  llbrt:  optic  sysiems. 


Introduction:  Arrays  in  which  several  phoUxIcLcctors  drive  a  single 
broadband  load  are  becoming  increasingiy'  important.  They  were 
originally  conceived  as  part  of  the  opti^dectronic  switching 
method  [1],  where  they  are  used  as  column  vectors  in  the  construe- 
lion  of  a  cross-point  switch  or  signal  processor  [2].  Metal-semicoii- 
ductor-metal  (MS\f;  phoUHktcctors  arc  particularly  well  suited  to 
this  application  [3],  due  to  their  large  bandwidth,  large  active  airea, 
simple  i^anar  structure,  and  low  capacitance.  Tn  addition,  their 
symmetiv'  allows  bipc^lar  opt^ation  and  implies  good  on-off  isola¬ 
tion  [2].  .A  scco.nd  emerging  application  of  photodetector  arrays  k 
in  the  microwave  Obre  optics  field,  wh«:e  the  link  insertion  loss 
and  si^al-to-noise  ratio  may  be  limited  by  the  power-handling 
capabilities  of  photodetectors  in  the  receiver  [4,  5], 

We  report  a  bandwidth  of  12GHz  for  an  integrated  array  of 
eiglii  lop-iilinninatcd  MSM  detectors,  packaged  on  a  dielecliic 
carrier  with  coplanar  w^aveguide  aaess  lines.  To  our  knowledge, 
this  is  the  hipest  speed  ratted  for  an  array  of  detectors  with 
individually  accessible  bias  lines,  as  required  for  optoelectronic 
switching.  Previously,  Lhi  el  al  [6]  reported  5  GHz  operation  for 
an  airay  of  4  PIN  drtcctors. 


Fig.  1  Photop'oph  of  packa:i^ed  detector  array 


Fabrication:  Fig.  1  show's  a  photograph  of  the  array  r^rled  here. 
Fight  incegrdted  MSM  photodetectors,  with  (26iJmy  active 
regions,  are  attached  to  a  oonmion  bus,  with  an  inler-detcctor 
spacing  of  250jjm.  The  opposite  contact  of  each  MSM  provides 


7'heory:  The  basic  theory  of  travcllmg-w-avc  detector  anays  is 
available  elsewhere  [5].  Essentially,  the  capacitance  of  a  single  long 
detector  or  several  parallel  photodetectors  can  be  compensated  by 
the  inductance  of  an  appropriate  intei'connect  netw-ork,  such  lhat 
an  artificial  transmission  line  matched  to  the  system  impedance  is 
formed.  This  is  illustrated  schematically  in  the  lumped-element 
model  of  Fig  2.  The  capacitance  derives  primarily  from  the  pho- 
todetcctors,  the  resistive  loss  terms  are  rnainly  due  to  finite  con¬ 
ductivity  and  charge  storage  in  the  semiconduaor  layers,  and  the 
inductance  is  provided  by  an  appropriate  intcixx)mioct  network,  as 
mentioned. 

At  microwave  IfequerK^ies,  it  is  more  accurate  to  treat  ihe  inter- 
connecting  lines  with  a  distributed  model.  For  the  arrays  discussed 
here,  the  interconnect  network  is  provided  by  an  on-wafer  micros- 
tr^  line,  which  forms  the  common  bus.  The  bus  is  100|im  wide 
and  SOOpm  from  the  ground  plane.  The  impedance  of  this  line  is 
approximately  80Q,  and  the  effective  propagation  index  is  approx¬ 
imately  2.78.  T  he  photodctcctor  array  closely  resembles  a  periodi¬ 
cally  loaded  transmission  line.  The  matdiing  oriteria  for  a 
terminated  periodic  structure  are  weU-knovvn,  and  iciid  to  the  con¬ 
dition  [5]: 


i 


I 


where  dr,  A.,  aJid  ni  arc  the  intcr-dctcctor  length,  characteristic  : 
impedance,  and  effective  propagation  index  of  the  high-impedance 
inteTConnect,  respectively,  Q  is  the  capacitance  associated  vrilli  i 
each  photodctcctor,  and  Zo  is  the  impedance  of  the  external  sys-  r 
tern.  Eqn.  1  i$  valid  for  w'avclcngths  greater  than  the  periai  of  the  i 
loaded-line  structure.  ; 

The  capadtance  of  discrete  MSM  photodetectors  is  often  esti-  ^ 
mated  using  amformal  transformation  [7],  Those  methods  predict 
the  capacitance  assTxiatcd  with  tlic  interdigitated  tmgers  only, 
asSLoning  no  excess  charge  storage  in  the  semiconductor  layers.  In  I 
real  MSM  detectors,  two  other  sources  of  parasitic  capacitance  ; 
must  be  considered  These  are  the  capacitance  associated  with  the  1 
tx^nd  pads  and  interconnect  metal  [7],  and  the  so-called  depletion  | 
capacitance  [8]  arising  from  charge  storage  deep  in  the  semicon- 
ductyr  bulk  or  at  hcterojunctioiofi.  By  fonnmg  these  aiTiiys  on  a  | 
rcsistiw  homojunction  layer,  grown  on  a  semi-insulating  GaAs  I 
substrate,  the  depletion-capadtancc  is  minimised  [7].  Tliis  simpli¬ 
fies  the  design  of  an  impeclancc-matchcd  array,  since  both  the  fin-  ■ 
ger  and  bond-pad  capacitance  contributions  may  be  predicted  v 
from  layout  geometry.  * 


Experimented  rextdis:  T)tc  capadtance  of  detectors  (in  the  anny)  \ 
w^as  detennined  using  an  on-wafer  probe,  and  by  fitting  reOec-  ‘ 
tioii  data  to  a  simple  equivalent  dreuif  [7].  The  capacitance  was 
approximately  40-^fF,  with  little  bias  depemienoc.  From  cqn.  1,  , 
and  using  the  high-impedance  line  characteristics  mentioned 
above,  the  lequired  detector  capadtance  for  matcliing  the  array  to  ; 
a  50Q  system  is  45 fF.  Assuming  perfect  matching  (Q  =  45fF),  the 
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ABSTRACT 

A  two-element  Yagi-Uda  array  with  two  physically  equal-length  slots  tuned  by  DC 
bias  is  presented.  By  changing  the  DC  bias  of  two  reactive  FET  circuits  imposed  upon 
both  ends  of  the  slot,  one  of  the  two  slots  can  be  made  shorter  electromagnetically  and 
serves  as  a  director  so  that  a  two-element  Yagi-Uda  antenna  array  is  achieved. 


L  INTRODUCTION 

Printed  slot  antennas  fed  by 
microstrip  lines  have  been  extensively 
used  in  radar  and  communication 
systems,  mainly  due  to  their  wider 
bandwidth  in  comparison  with  patch 
antennas,  low  cross-polarization,  as  well 
as  good  isolation  with  the  active 
circuitry.  Recently,  slot  antennas 
integrated  with  active  components  have 
expanded  their  application  even  more 

[1] .  In  general,  the  integration  of  active 
devices  provides  a  control  mechanism 
for  antenna  performance  that  is  not 
possible  with  passive  antennas. 

Traditionally,  a  two-element 
Yagi-Uda  antenna  array  is  made  of  two 
dipoles  or  slots  with  different  physical 
lengths,  in  which  one  of  the  two  slots  is 
shorter  or  longer  than  the  other  to  serve 
as  a  director  or  a  reflector,  respectively 

[2] .  Such  configuration,  however,  has 
the  direction  of  the  mainbeam  fixed  as 
the  geometrical  structure  is  decided.  In 

[3] ,  a  slot  antenna  with  electronically 
tunable  length  was  reported.  In  this 
paper,  a  novel  idea  of  utilizing  two 
physically  equal-length  slot  anteimas  in 
which  either  one  slot  is  tunable  in 


electrical  length  to  achieve  a  Yagi-Uda 
antenna  array  is  presented. 

n.  DESIGN 

Figure  1.  shows  the 
configuration  for  a  two-element  Yagi- 
Uda  antenna  array.  The  FET  is  used  to 
provide  a  reactive  output  as  a  one-port 
circuit.  The  varying  reactance  provided 
through  the  change  in  the  Schottky 
barrier  as  the  gate  bias  is  changed  can  be 
parallelly  inserted  to  the  slot  antenna 
through  electromagnetic  coupling.  With 
this  additional  reactance,  an 
electronically  tunable  slot  can  be 
obtained.  The  input  impedance  of  the 
reactive  FET  circuit  is  designed  to  be  (0- 
j20)  Q  around  10  GHz,  which 
corresponds  to  0.8  pF  at  10  GHz.  This 
input  impedance  is  to  be  placed  at  the 
ends  of  the  slot,  where  the  impedance  is 
low,  so  that  it  does  not  significantly 
affect  the  field  distribution  in  the  slot 
antenna.  The  impedance  at  the  very  end 
of  the  slot,  however,  is  0  £2.  Therefore, 
an  offset  of  0.067X  between  the  very  end 
of  the  slot  and  the  center  of  the 
microstrip  line  is  implemented. 


In  order  to  provide  a  RF  short  at 
the  crossing  point  of  the  slot  and 
microstrip  line,  a  capacitor  is  seriesly 
connected  to  the  end  of  the  microstrip 
line  before  it  is  grounded.  Two  iXjw  long 
slots  are  presented  here.  The  spacing 
between  two  slot  lines  is  0.25  Xo. 

A  two-element  Yagi-Uda 
antenna  array  is  built,  as  shown  in 
Figures  2  and  3.  As  the  DC  biases 
imposed  upon  the  gate  and  drain  of  the 
FET  are  changed,  the  slot  #2  is  increased 
in  its  resonant  frequency  while 
decreased  in  its  electrical  length,  so  that 
it  can  serve  as  a  director.  Thus,  a  two- 
element  Yagi-Uda  antenna  array  is 
achieved. 

m.  RESULTS 

The  slot,  the  reactive  FET  circuit 
and  the  feed  line  of  the  external  source 
are  fabricated  on  RT/duroid  5870  by 
ROGERS  with  er=2.33  and  31  mils  in 
thickness,  and  GaAs  MESFET 
transistors  manufactured  by  NEC 
(NE76184A)  are  used  as  the  active 
devices  in  the  circuit.  The  capacitors 
used  are  2.4  pF. 

The  measured  result  of  Sll  of 
the  reactive  FET  circuit  is  shown  in 
Figure  4.  The  tuning  voltages  are  OV  to 
-1.54V  for  the  gate  voltage  (Vgs)  and  OV 
to  0.36V  for  the  drain  voltage  (Vds).  The 
variation  of  the  resonant  frequency 
mainly  depends  on  the  change  of  Vgs. 
Since  the  reactive  FET  has  the  reactive 
element  on  the  source  terminal,  the 
operation  is  not  quite  stable  as  Vds 
increases.  Under  small  Vds,  however,  it 
is  possible  to  obtain  wide  tuning  range. 


A  tuning  range  of  8%  in  resonant 
frequency  of  the  reactive  FET  circuit  is 
obtained.  With  this,  we  are  able  to 
change  the  electrical  length  of  the  slot  by 
approximately  the  same  amount.  The 
resulting  H-plane  radiation  pattern  of  the 
two-element  Yagi-Uda  antenna  array  is 
shown  in  Figure  5.  With  a  DC  bias  of 
Vgs=- 1.535  V  and  Vds=0.35  V,  a  front- 
to-back  ratio  of  10  dB  is  observed. 
Without  DC  bias,  only  a  slight 
directivity  is  observed,  and  the  measured 
front-to-back  ratio  is  3  dB. 

IV.  CONCLUSION 

In  this  work,  a  two-element 
Yagi-Uda  antenna  array  with  two 
physically  equal-length  elements  is 
designed.  Reactive  FET  circuits  are 
utilized  to  tune  the  electrical  length  of 
one  of  the  two  slots.  A  tunable  two- 
element  Yagi-Uda  antenna  array  is 
demonstrated,  and  a  10  dB  front-to-back 
ratio  is  obtained.  This  work  provides  a 
novel  idea  that  the  mainbeam  of  Yagi- 
Uda  antenna  arrays  can  be  made 
switchable  by  the  utilization  of  reactive 
FET  circuits. 
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Figure  4.  Shiftiiig  of  Sll  of  the  reactive  FET 
circuit  when  DC  bias  is  imposed. 
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Figure  5.  Radiation  Pattern  in  H-plane 
with  DC  bias;  Vgs=-1.535  V, 
Vds=0.35V. 


Table  1.  Comparison  of  front>to-back 
ratio  with  and  without  bias. 


Without  bias 

3dB 

With  bias 

10  dB 
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sistently  produce  lower  compression  currents 
than  time-domain  measurements  for  all  tested 
photodetectors.  The  difference  can  be  ex¬ 
plained  by  the  delayed  evolution  of  the  space- 
charge  field  and  by  the  fact  that  saturation  is 
frequency  dependent.  High-frequency  signals 
become  saturated  at  much  lower  levels  than 
low-frequency  signals.  Frequency-domain 
measurements  are  usually  performed  at  high 
frequencies,  whereas  time-domain  measure¬ 
ments  represent  the  average  saturation  level 
over  a  broad  frequency  range.  We  report  de¬ 
tails  of  measurements,  including  the  bias  de¬ 
pendence,  simulation  results,  buildup  time  of 
the  space-charge  field,  and  comparison  be¬ 
tween  the  time-  and  frequency-domain  mea¬ 
surements. 
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Long-wavelength  velocity-matched 
distributed  photodetectors 
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High-power,  high-speed  photodetectors  are  a 
key  component  in  microwave  fiber-optic  links, 
optoelectronic  generation  of  microwaves,  and 
millimeter  waves.  ^  Previously,  we  have  re¬ 
ported  a  GaAs/AJGaAs  velocity-matched  dis¬ 
tributed  photodetector  (VMDP)  that  operates 
at  860  nm  and  demonstrated  its  potential  for 
high-saturation  photocurrents.^  For  applica¬ 
tions  in  rf  photonic  systems,  however,  InP- 
based  long- wavelength  photodetectors  operat¬ 
ing  at  1.3  or  1.55  jxm  are  required.  In  this 
paper,  we  report  on  the  experimental  results  of 
InGaAs/InAlGaAs/InP  VMDP. 


CThM3  Fig.  2.  Contour  plot  of  the  funda¬ 
mental  mode  field  amplitude  in  the  passive  opti¬ 
cal  waveguide  using  BPM  simulation.  See  Fig.  1 
for  layer  structure. 


The  schematic  structure  of  the  VMDP  is  illus¬ 
trated  in  Fig.  1.  Active  metal-semiconductor- 
metal  (MSM)  photodiodes  are  periodically  dis¬ 
tributed  on  top  of  a  passive  optical  waveguide. 
Optical  signal  is  evanescently  coupled  from  the 
passive  waveguide  to  the  active  MSM  photo¬ 
diodes.  Photocurrents  generated  from  each 
MSM  photodiode  are  added  in  phase  through  a 
50-n  coplanar  strip  (CPS)  microwave  transmis¬ 
sion  line  that  is  velocity  matched  to  the  optical 
waveguide.  The  MSM  photodiodes  serve  two 
functions  here:  they  generate  photocurrents  and 
provide  the  periodic  capacitance  loading  needed 
for  velocity  matching.  The  VMDP  design  allows 
the  passive  waveguide,  the  active  photodiodes, 
and  the  microwave  coplanar  strips  to  be  inde¬ 
pendently  optimized.  In  Fig.  1,  the  active  MSM 
photodiode  consists  of  an  InGaAs  absorption 
layer,  an  InGaAs/InAlAs  graded  superlattice  lay¬ 
ers,  and  an  InAlAs  Schottky-barrier  enhance¬ 
ment  layer.  At  10-V  bias,  a  dark  current  of  1 90  pA 
is  measured  for  an  1 1  X  48-p.m^  MSM  photo¬ 
diode.  For  a  1.2-mm-long  VMDP  with  12  11  X 
48-p.m^  MSM  photodiodes,  a  dark  current  of  25 
nA  is  measured  at  10-V  bias. 

Figure  2  shows  the  contour  plot  of  the  fun¬ 
damental  mode  field  amplitude  in  the  passive 
waveguide  using  beam  propagation  method 
(BPM)  simulation.  The  VMDP  is  designed  so 
that  only  the  fundamental  mode  exists  in  the 
optical  waveguide  and  the  active  photodiode 
regions.  Simulation  results  with  the  BPM 
method  indicate  a  4%  scattering  loss  and  7% 
absorption  loss  per  photodiode.  Scattering  loss 
can  be  reduced  to  1%  per  diode  by  optimizing 
the  geometric  structure  of  the  active  photo- 


Frequency  (GHz) 

CThM3  Fig.  3.  Measured  frequency  response 
of  long- wavelength  VMDP  (12  photodiodes,  1.2 
mm  long). 


diode.  The  device  under  test  is  1,2  mm  long 
and  consists  of  12  MSM  photodiodes.  In  this 
experiment,  optical  lithography  is  employed  to 
pattern  the  MSM  photodiodes.  Both  the  fin¬ 
gers  and  the  spacing  are  1  p.m  wide.  After 
antireflection  (AR)  coating,  a  quantum  effi¬ 
ciency  of  34%  was  measured  (this  includes  the 
coupling  loss  from  optical  fiber).  High- 
saturation  current  measurement  with  a  high- 
power  erbium-doped  fiber  amplifier  is  in 
progress. 

The  frequency  response  of  the  VMDP  is 
characterized  by  the  optical  heterodyne 
method.^-5  The  system  consists  of  two  external 
cavity  tunable  lasers  at  1.55  p.m;  the  frequency 
of  each  laser  can  be  tuned  in  1-GHz  steps.  The 
optical  signals  are  combined  by  a  3-dB  coupler 
and  coupled  to  the  VMDP  through  a  fiber 
pickup  head.  The  microwave  signal  generated 
by  optical  mixing  in  the  VMDP  is  collected  at 
the  output  end  of  the  CPS  by  a  50-GHz  picop- 
robe  (GGB  Industries),  which  is  connected  to 
an  rf-power  sensor  and  monitored  by  an  rf- 
power  meter.  The  calibrated  frequency  re¬ 
sponse  of  long- wavelength  VMDP  is  shown  in 
Fig.  3.  A  3-dB  bandwidth  frequency  of  18  GHz 
is  measured.  By  scaling  down  the  MSM  to  0. 1 
|xm  scale,  a  frequency  response  >100  GHz  is 
expected. 

In  summary,  we  have  experimentally  demon¬ 
strated  a  high-speed,  high-power  long- 
wavelength  VMDP  using  optical  lithography.  A 
3-dB  bandwidth  of  18-GHz  and  an  overall  quan¬ 
tum  efficiency  of  34%  have  been  achieved. 
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CThM3  Fig.  1.  Schematic  structure  of  long- wavelength  velocity- matched  distributed  photodetector 
(VMDP). 
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the  Fourier  transform  of  the  EO-signal  is  shown  in  Fig.  3.  To  cor¬ 
rect  the  frequency  response  for  the  effect  of  finite  duration  of  the 
optical  sampling  and  excitation  pulse,  the  Fourier  transforms  of 
their  measured  autocorrelation  traces  are  used.  A  further  effect  on 
the  system  response  results  from  the  finite  transit  time  of  the  sam¬ 
pling  pulse  in  the  electro-optically  active  InP  substrate.  Its 
response  characteristic  is  calculated  by  accounting  for  the  spatial 
distribution  of  the  electric  field  of  the  CPW  line,  assuming  quasi- 
TEM  behaviour  of  the  propagating  signal.  The  resulting  frequency 
response  of  the  EO-sampling  system  as  well  as  the  corrected 
response  of  the  detector  are  shown  in  Fig.  3.  A  3dB  bandwidth  of 
70GHz  results  from  the  corrected  response  characteristic.  The 
comparison  of  the  corrected  response  with  the  RC-linuted 
response  (also  shown  in  Fig.  3),  which  is  obtained  from  the  capac¬ 
itance  (29fF)  and  series  resistance  (5fi)  of  the  detector  when  con¬ 
nected  to  a  son- transmission  line,  shows  that,  in  the  frequency 
range  120  and  170GHz,  the  response  is  dominated  by  the  RC- 
effect.  Moreover,  the  comparison  indicates  that  the  transit-time 
limited  bandwidth  of  the  devices  exceeds  150GHz.  The  corrected 
response  slightly  exceeds  the  RC-limit  at  very  high  frequencies  in 
Fig.  3.  This  may  be  attributed  to  the  imcertainty  of  the  correction 
procedure,  which  is  estimated  to  be  ±0.5dB. 


Fig.  3  Frequency  response  characteristics 

■  corrected  response 
•  uncorrected  response 
O  response  of  measuring  system 
□  RC  limited  response 


Typically,  the  dark  currents  of  the  devices  range  between  1  and 
2nA  at  2V  bias.  The  CW  quantum  efficiency  of  uncoated  devices 
at  1.55  and  1.3Mm  is  found  to  be  3.5  and  10%,  respectively.  The 
application  of  a  silicon  nitride  AR  coating  leads  to  an  enhanced 
efficiency  of  6.7  and  17%,  respectively.  Note  that  the  beha,viour  of 
the  efficiency  does  not  obey  the  geometric-optics  approximation. 
Rather,  its  treatment  requires  the  application  of  exact  diffraction 
theory  [4].  The  polarisation  dependence  of  the  quantum  efficiency 
is  found  to  be  <  0.2dB.  Moreover,  the  AR  coating  gives  rise  to  a 
larger  capacitance  of  40fF  due  to  the  increase  in  the  effective  die¬ 
lectric  constant.  Thus,  preserving  the  bandwidth  of  70GHz  with 
AR  coated  devices  will  require  a  reduction  of  the  active  area  by 
-30%,  As  recently  demonstrated  in  conjunction  with  InGaAs 
MSM  detectors  [6],  a  drastic  improvement  of  the  quantum  effi¬ 
ciency  without  a  significant  loss  to  bandwidth  performance  can  be 
achieved  by  rear  illumination  in  combination  with  resonant  cavity 
enhanced  (RCE)  absorption.  Such  an  RCE  configuration  is 
expected  to  result  in  a  quantum  efficiency  of -40%  at  1.55jjm. 

Conclusion:  Front-side  illuminated  14|jm  active  area  diameter 
InGaAs  photodetectors  for  high-speed  and  polarisation-insensitive 
operation  have  been  fabricated.  The  devices  possess  0.2|jm  feature 
size  finger  electrodes  with  a  semicircular  shape.  A  pulse  response 
time  of  3.8ps  at  2V  bias  has  been  measured  by  EO-sampling  at 
1.55pm  wavelength.  The  corrected  frequency  response  exhibits  a 
3dB  bandwidth  of  70GHz. 
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Long  wavelength  velocity-matched 
distributed  photodetectors  for  RF  fibre  optic 
links 

T.  Chau,  L.  Fan,  D.T.K.  Tong,  S.  Mathai,  M.C.  Wu, 
D.L.  Sivco  and  A.Y.  Cho 

An  InP-based  long  wavelength  velocity-matched  distributed 
photodetector  with  metal-semiconductor-metal  photodiodes  is 
experimentally  demonstrated.  A  3dB  bandwidth  of  18GHz  and 
an  external  quantum  efficiency  of  0.42  A/W  have  been  achieved. 

Introduction:  High  power,  high  frequency  photodetectors  are  a  key 
component  for  high  performance  microwave  fibre  optic  links  [1  -  3]. 
High  optical  power  in  externally  modulated  links  can  greatly 
enhance  the  link  gain,  signal-to-noise  ratio,  and  the  spurious-free 
dynamic  range  [4].  It  is  also  very  useful  for  millimeter-wave  gener¬ 
ation  by  photomixing.  Previously,  we  have  reported  a  GaAs/ 
AlGaAs  velocity-matched  distributed  photodetector  (VMDP) 
operating  at  860 nm  wavelength,  and  demonstrated  its  potential 
for  high  saturation  photocurrent  [1].  For  applications  in  RF  phot¬ 
onic  systems,  however,  InP-based  long-wavelength  photodetectors 
operating  at  1.3  or  1.55  pm  are  required.  The  successful  fabrication 
of  a  long-wavelength  ^  reported  in  [5].  In  this  Letter, 

we  report  on  the  experimental  results  and  RF  performance  of  the 
InGaAs/InAlGaAsAnP  VMDP. 

Device  structure  and  fabrication:  The  schematic  structure  of  the 
VMDP  is  illustrated  in  Fig.  la.  Active  metal-semiconductor-metal 
(MSM)  photodiodes  are  periodically  distributed  on  top  of  a  pas¬ 
sive  optical  waveguide.  The  optical  signal  is  evanescently  coupled 
from  the  passive  waveguide  to  the  active  MSM  photodiod^. 
Photocurrents  generated  from  the  MSM  photodiodes  are  added  in 
phase  through  a  500  coplanar  strip  (CI^)  microwave  transmis¬ 
sion  line  that  is  velocity-matched  to  the  optical  waveguide.  The 
MSM  photodiodes  serve  two  functions:  generating  photocurrents 
and  providing  the  periodic  capacitance  loading  need^  for  velocity 
matching.  The  VMDP  design  allows  the  passive  waveguide,  the 
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active  photodiodes,  and  the  microwave  coplanar  strips  to  be  inde¬ 
pendently  optimised.  The  active  MSM  photodiodes  consist  of  an 
InGaAs  absorption  layer,  InGaAs/InAlAs  graded  superlattice  lay¬ 
ers  an  InAlAs  Schottky-barrier  enhancement  layer,  and  a  Ti/Au 
contact  with  200A/2000A  thickness.  The  fabrication  proc^  is  as 
follows:  first,  Ti/Au  interdigitated  fingers  with  1.5pm  width  and 
spacing  are  patterned  using  standard  photolithography  and 
liftoff  techniques,  then  active  mesas  for  individual  photodiodes  ^ 
defined  by  wet  etching  down  to  the  InAlAs  upper  cladding  II.  Fig. 
\b  shows  the  cross-section  of  the  VMDP  after  mesa  etching.  A 
passive  ridge  waveguide  with  a  ridge  height  of  0.1pm  is  employed 
to  connect  the  photodiodes.  After  mesa  and  waveguide  etching, 
the  CPS  microwave  transmission  line  is  patterned  using  a  standard 
liftoff  technique. 
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Fig.  1  Schematic  structure  of  long-wavelength  VMDP  and  cross-section 
after  mesa  etching 


fl  Schematic  structure  of  VMDP 
b  Schematic  cross-section  of  VMDP  after  mesa  etching 


=( 

— — 

rwal 

Fig.  2  Contour  plot  of  fundamental  mode  field  amplitude  in  photodiode 
region  using  BPM  simulation 

See  Fig.  I  for  detailed  layer  structure 


The  beam  propagation  method  (BPM)  is  used  to  design  and 
simulate  the  optical  performance  of  the  VMDP.  A  large-core  opti¬ 
cal  waveguide  is  employed  to  reduce  the  coupling  loss  between  the 
passive  waveguide  and  the  photodiode  region  as  well  as  to  reduce 
the  optical  power  density  in  the  absorption  region.  Only  the  fun¬ 
damental  mode  exists  in  both  the  passive  wavegvude  and  the  pho¬ 
todiode  regions.  Fig.  2  shows  the  optical  field  distribution  of  the 


VMDP  in  the  photodiode  region.  Most  of  the  optical  energy  is 
concentrated  in  the  waveguide  core.  The  optical  absorption  and 
the  coupling  loss  per  photodiode  for  the  VMDP  with  11pm  long 
photodiodes  and  89pm  spacing  between  photodiodes  are  estimated 
to  be  6.2%  and  1.8%,  respectively,  by  the  BPM  simulation.  The 
coupling  loss  can  be  further  improved  by  reducing  the  width  of 
the  active  mesa. 


Device  characteristics!  The  finished  devices  are  cleaved  and 
mounted  on  copper  heat  sinks  before  testing.  The  VMDP  exhibits 
a  very  low  dark  current:  190pA  at  lOV  bias  for  individual  photo¬ 
diodes  (11  X  48Mm2),  and  25 nA  for  the  1.2nmi  long  VMDP  with 
12  photodiodes.  The  excess  dark  current  in  the  VMDP  is  attrib¬ 
uted  to  the  leakage  current  under  the  CPS  electrodes,  which  h^  a 
much  larger  area  than  the  active  photodiodes.  This  can  be  elimi¬ 
nated  by  placing  the  electrodes  on  thin  dielectric.  The  external 
quantum  efficiency  is  measured  to  be  0.42A/W  after  anti-reflection 
(AR)  coating.  The  dominant  loss  comes  from  the  coupling  from 
the  optical  fibre  to  the  VMDP.  All  devices  exhibit  linear  respon- 
sivity  up  to  12mA  of  photocurrent. 


Fig.  3  Measured  frequency  response  of  long  wavelength  VMDP 
12  photodiodes,  total  length  =  1,2  mm 


The  microwave  performance  of  the  VMDP  is  measured  by  an 
HP  8510C  network  analyser.  The  microwave  return  loss  (Su)  is 
<-^22dB  from  0.1  to  40GHz.  The  characteristic  impedance  of  the 
VMDP  is  matched  very  well  to  500  (within  4%)  for  the  entire  fre¬ 
quency  range.  The  frequency  response  of  the  VMDP  is  character¬ 
ised  by  the  optical  heterodyne  method  [6,  7]  two  external 
cavity  tunable  lasers  at  1.55pm,  The  optical  signals  are  combined 
by  a  3dB  coupler,  and  coupled  to  the  VMDP  through  a  fibre 
pickup  head.  The  output  microwave  is  collected  by  a  SOGHz  pico- 
probe  (GGB  Industries)  and  monitored  by  an  RF  power  meter. 
The  calibrated  frequency  response  of  the  VMDP  is  shown  in  Fig. 
3.  A  3dB  bandvwdth  of  18GHz  is  measured.  The  bandwidth  is 
limited  by  the  transit  time  of  the  MSM  photodiodes.  By  scaling 
dovm  the  MSM  to  a  0.1  pm  scale,  much  higher  bandwidth 
(>  lOOGHz)  can  be  obtained. 

The  maximum  photocurrent  of  12mA  is  currently  limited  by 
thermal  damage  to  the  MSM  photodiode.  The  device  remains  lin¬ 
ear  up  to  the  maximum  photocurrent.  Inspection  of  the  failed 
device  by  energy  dispersive  X-ray  analysis  (EDX)  indicates  that 
gold  diffiision  into  the  semiconductor  is  the  main  failure  mecha¬ 
nism.  Improving  the  diffusion  barrier  of  the  Schottky  contact 
should  lead  to  even  higher  photocurrent. 

Conclusion:  We  have  experimentally  demonstrated  a  long  wave¬ 
length  velocity -matched  distributed  photodetector  with  twelve 
metal-semiconductor-metal  photodiodes.  A  3dB  bandwidth  o\ 
18GHz  and  an  external  quantum  efficiency  of  0.42A/W  have  beer 
achieved. 
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Optoelectronic  clock  recovery  circuit  using 
resonant  tunnelling  diode  and  uni-travelling- 
carrier  photodiode 

K.  Murata,  K.  Sano,  T.  Akeyoshi,  N.  Shimizu,  E.  Sano, 
M.  Yamamoto  and  T.  Ishibashi 

An  optoelectronic  clock  recovery  circuit  is  reported  that 
monolithically  integrates  a  resonant  tuimelling  diode  and  a  uni- 
travelling-carrier  photodiode.  The  integrated  circuit  extracts  an 
electrical  lL55GHz  dock  signal  from  11.55Gbit/s  RZ  optical 
data  streams  in  a  wide  locking  range  with  low  power  dissipation. 
Furthermore,  the  extraction  of  a  subharmonic  clock  from  23.1 
and  46.2Gbit^s  input  data  streams  is  also  confirmed. 

Introduction:  Because  of  the  growth  in  multimedia  services,  ultra- 
high-speed  optical  receivers  will  be  indispensable  for  future  back¬ 
bone  networks.  To  realise  a  small  high-speed  receiver  circuit  with 
low  power  consumption,  an  optoelectronic  circuit  using  resonant 
tunnelling  diodes  (RTDs),  and  a  new  type  of  wideband,  high-satu¬ 
ration  power  photodetector  called  the  uni-travelling-carrier  photo¬ 
diode  (UTC-PD)  [1]  is  a  promising  candidate  because  of  its 
inherent  high-speed  operation.  For  example,  an  RTDAJTC-PD 
integrated  circuit  (IC)  demultiplexing  an  80Gbit/s  optical  signal 
into  40  Gbit/s  electrical  signals  at  an  extremely  low  power  of 
7.75mW  has  been  demonstrated  [2].  A  clock  recovery  circuit  is 
another  key  component  for  realising  the  receiver  circuit.  Several 
kinds  of  optical  injection-locked  RTD  oscillator  [3  -  5]  have  been 
reported,  and  phase  locking  to  optical  input  has  b^n  demon¬ 
strated,  However,  clock  extraction  from  an  optical  data  signal, 
such  as  a  pseudo-random  bit  stream  (PRBS),  and  subharmonic 
clock  extraction  have  not  been  reported. 

This  Letter  describes  an  optoelectronic  clock  recovery  circuit 
that  monolithically  integrates  an  RTD  and  a  UTC-PD  on  an  InP 
substrate.  The  principle  of  synchronisation  is  optical  injection¬ 
locking  of  a  free-running  oscillator.  The  fabricated  circuit  success¬ 
fully  extracted  an  electrical  11. 55  GHz  clock  signal  from  an 


1 1.55  Gbit/s  RZ  optical  input  data  stream  in  a  wide  locking  range 
with  low  power  dissipation.  Furthermore,  the  IC  exhibited  sub- 
harmonic  clock  extraction  from  23.1  Gbit/s  and  46.2Gbit/s  input 
data  streams. 


Fig.  1  Circuit  diagram  of  clock  recovery  circuit 


Circuit  configuration:  Fig.  1  shows  a  circuit  diagram  of  the  pro¬ 
posed  clock  recovery  circuit.  The  dicuit  consists  of  an  oscillator, 
which  is  constructed  with  an  RTD  and  a  transmission  line,  and  a 
UTC-PD.  In  the  RTD  oscillator,  the  RTD  is  biased  in  the  nega¬ 
tive  differential  resistance  (NDR)  region,  and  the  collector  is  con¬ 
nected  to  the  transmission  line  whose  other  terminal  is  connected 
to  the  ground.  The  essential  synchronisation  principle  is  injection¬ 
locking  of  the  RTD  oscillator  using  the  photocurrent  generated  by 
the  UTC-PD.  Here,  the  use  of  the  UTC-PD  is  important  to  main¬ 
tain  fast  photoresponse  at  low  bias  voltages  [6]  corresponding  to 
the  NDR  region  of  the  RTD.  The  self-oscillation  of  the  RTD 
oscillator  is  theoretically  analysed  in  [7],  and  the  oscillating  fre¬ 
quency  is  inversely  proportional  to  the  propagation  delay  time 
of  the  transmission  line.  In  the  present  circuit,  the  transmission 
line  is  a  coplanar  wave  guide  monolithically  fabricated  on  the  InP 
substrate.  TTie  physical  length  and  characteristic  impedance  Zq  of 
the  transmission  line  were  1150pm  and  SOD,  respectively.  The 
active  areas  of  the  RTD  and  UTC-PD  were  6  and  iopm^,  respec¬ 
tively.  The  bias  voltage  V  was  +0.43  V;  the  chip  size  was  1.9  x 
0.5mm. 


Experimental  results:  The  IC  was  fabricated  by  monolithically 
integrating  an  InGaAs/AlAs/InAs  RTD  and  an  InP/InGaAs  UTC- 
PD  on  an  semi-insulating  InP  substrate  [2,  8].  The  peak  and  valley 
current  of  the  RTD  were  7.4mA  at  0.35  V  and  0.6niA  at  0.7  V, 
respectively.  The  responsivity  of  the  UTC-PD  was  0.26A/W  at  a 
wavelength  of  1.55pm,  and  the  3dB  bandwidth  around  the  bias 
voltage  Kwas  --  80GHz. 


oaft'3 

Fig.  3  Locking  range  of  clock  recovery  circuit 

#  fundamental 
■  first-order  harmonic 
A  third-order  harmonic 


The  IC  was  tested  on  a  wafer.  The  input  RZ  optical  data  stream 
(k  =  1.55  pm)  was  generated  by  an  electro-optic  pulse  pattern  gen¬ 
erator  that  can  output  a  10-80 Gbit/s  opticd  data  stream  [9].  The 
pulsewidth  of  the  optical  data  input  was  <  lOps,  and  the  data 
sequence  was  PN  2^^-!  PRBS.  The  input  optical  signal  illuminated 
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Distributed  Photodetectors  For  RF  Fiber  Optic  Links 
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Abstract 

Improved  performance  of  InP-based  long  wavelength  velocity-matched  distributed  photodetector  (VMDP)  with 
metal-semiconductor-metal  photodiodes  is  experimentally  demonstrated.  A  3-dB  bandwidth  of  13  GHz  and  an 
external  quantum  efficiency  of  0.57  A/W  have  been  achieved. 

1.  Introduction 

High  power,  high  frequency  photodetector  is  a  key  component  for  high  performance  microwave  fiber  optic 
links  [1-3],  High  optical  power  in  externally  modulated  links  can  greatly  enhance  the  link  gain,  signal-to-noise  ratio, 
and  spurious-free  dynamic  range  [4].  In  conventional  photodetectors,  there  is  a  trade-off  between  the  saturation 
photocurrent  and  the  device  bandwidth.  High  power  device  requires  a  laige  absorption  volume,  thus  the  device  is 
usually  large,  resulting  in  large  RC  time  or  long  carrier  transit  time  that  limits  the  device  bandwidth.  In  contrast, 
conventional  high-speed  photodetectors  with  small  RC  time  and  small  transit  time  usually  have  low  saturation 
photocurrents.  The  velocity-matched  distributed  photodetectors  (VMDP)  is  a  novel  type  of  travelling  wave 
photodetector  which  can  achieve  both  high  bandwidth  and  high  saturation  power.  It  was  first  proposed  in  1993  [6], 
and  has  been  successfully  demonstrated  experimentally  [1,5].  The  successful  fabrication  of  long-wavelength  VMDP 
for  use  in  1.3  or  1.55  |im  RF  Photonics  systems  was  first  reported  in  [5].  In  this  paper,  we  report  on  the  performance 
of  improved  long  wavelength  VMDP  with  new  design  and  fabrication  procedures.  A  3-dB  bandwidth  of  13  GHz  and 
an  external  quantum  efficiency  of  0.57  A/W  have  been  achieved. 

II.  Device  Structure  and  Fabrication 

The  schematic  structure  of  the  VMDP  is  illustrated  in  Figure  1.  A  passive  (non-absorbing)  optical  waveguide 
is  used  to  serially  connect  an  array  of  periodically  spaced  metal-semiconductor-metal  (MSM)  photodiodes.  Light  in 
the  optical  waveguide  is  evanescently  coupled  to  the  MSM  photodiodes.  The  photocmrents  are  added  in  phase  and 
collected  by  a  50G  coplanar  strips  (CPS)  microwave  transmission  line.  The  active  photodiodes  are  designed  to  have 
small  optical  confinement  factor  to  keep  them  bellow  saturation  imder  high  optical  illumination.  The  bandwidth  of 
the  VMDP  is  limited  by  that  of  the  individual  photodiode,  and  the  residual  velocity  mismatch.  Since  photocurrents 
are  collected  from  many  photodiodes  along  the  transmission  line,  the  individual  photodiode  does  not  need  to  have 
high  quantum  efficiency  and  therefore  can  be  made  small  and  fast.  The  MSM  photodiodes  serve  two  functions: 
generating  photocmrents  as  well  as  providing  the  periodic  capacitance  loading  needed  for  velocity  matching.  The 
VMDP  design  allows  the  passive  waveguide,  the  active  photodiodes,  and  the  microwave  cofrtanar  strips  to  be 
independently  optimized. 


0-7803-4936-9/98/$  10.00  ©  1998  IEEE 
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The  principle  of  the  VMDP  has  been  discussed  in  more  detail  in  [1].  We  have  made  three  major  changes:  (1) 
mesa  width  reduction  to  reduce  optical  coupling  loss  from  optical  waveguide  to  active  MSM  photodiodes;  (2)  nitride 


passivation  on  mesa  sidewall  and  underneath  the  large  CPS  electrodes  to  reduce  dark  current  and  improve  device 
reliability;  and  (3)  use  of  platinum  in  Schottky  metal  contacts  to  prevent  gold  diffusion  at  high  power  operation 

MSU 


Figure  1 .  Schematic  structure  of  long- wavelength  Figure  2.  Contour  plot  of  the  fundamental  mode  profile 

velocity-matched  distributed  photodetector  (VMDP)  in  the  photodiode  region  using  BPM 


Beam  propagation  method  (BPM)  is  used  to  design  and  simulate  the  optical  performance  of  the  VMDP.  Large- 
core  optical  waveguide  is  employed  to  reduce  the  coupling  loss  between  the  passive  waveguide  and  the  photodiode 
region  as  well  as  to  reduce  the  optical  power  density  in  the  absorption  region.  Only  fundamental  mode  exists  in  both 
the  passive  waveguide  and  the  photodiode  regions.  Figure  2  shows  the  optical  field  distribution  of  VMDP  in  the 
photodiode.  Most  of  the  optical  energy  concentrates  in  the  waveguide  core.  The  optical  absorption  and  the  coupling 
loss  per  photodiode  for  the  VMDP  shown  in  Fig.  1  are  estimated  to  be  8.3%  and  1.6%,  respectively,  by  the  BPM 
simulation. 
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Figures.  Cross  section  of  VMDP  wafer  Figure  4.  SEM  picture  of  a  single  MSM  photodiode 

after  mesa  etching  of  VMDP 


The  active  MSM  photodiodes  consists  of  InGaAs  absorption  layer,  InGaAs/InAlAs  graded  superlattice  layers; 
InAlAs  Schottky-barrier  enhancement  layer,  and  interdigitated  fingers.  The  fabrication  process  is  as  followed:  first, 
metal  alignment  maricers  are  patterned  on  the  substrate  for  subsequent  process.  Next,  active  mesas  for  photodiodes 
are  defined  by  wet  etching  down  to  the  InAlAs  Upper  Cladding  II.  Then  optical  ridge  waveguide  with  ridge  height 
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of  0.1|im  is  fonned  by  wet  etching.  After  mesa  and  waveguide  etching,  a  Si3N4  passivation  layer  is  deposited  to 
protect  the  mesa  edges.  The  nitride  directly  on  top  of  the  active  mesas  are  opened  for  Schottky  contacts. 
Interdigitated  Ti/Pt/Au  fingere  with  Ipm  finger  width  and  Ipm  finger  spacing  are  patterned  in  the  open  windows 
using  optical  lithography.  Finally,  the  CPS  microwave  transmission  line  is  fabricated  using  standard  liftoff 
technique.  The  cross  section  of  VMDP  wafer  after  mesa  etching  is  shown  in  Figiuu  3.  Figure  4  shows  the  scanning 
electron  micrograph  (SEM)  of  a  single  MSM  photodiode  in  the  VMDP. 

III.  Device  Characteristics 

All  devices  under  test  are  mounted  on  copper  heat  sinks.  The  temperature  is  set  at  19  ®C  using  a  temperature 
controller.  The  VMDP  exhibits  very  low  dark  current:  8.3  nA  at  lOV  bias  for  a  1-mm-long  VMDP  with  13 
photodiodes.  The  dark  current  has  been  reduced  by  10  times  compared  to  the  VMDP  in  [5].  It  is  attributed  to  the 
nitride  passivation  which  prevents  the  leakage  currents  through  the  sidewalls  and  the  CPS  electrodes.  The  external 
quantum  efficiency  is  measured  to  be  0.4  A/W.  With  anti-refiection  (AR)  coating,  it  could  reach  0.57  A/W.  The 
dominant  loss  comes  fi’om  the  coupling  fi’om  optical  fiber  to  the  VMDP.  Figure  5  shows  the  DC  responsivity  of  the 
new  device.  The  responsivity  of  the  VMDP  in  [5]  is  also  shown  for  comparison. 


Input  Optical  Power  (mW) 


Figures.  DC  responsivity  of  VMDP 
(without  AR  coating) 


Figure  6.  Schematic  of  optical  heterodyne  system  setup  for 
fi*quency  response  characterization  of  VMDP 


The  frequency  response  of  the  VMDP  is  characterized  by  optical  heterodyne  method  [7,8].  The  schematic  of 
the  experimental  setup  is  shown  in  Figure  6.  The  system  consists  of  two  external  cavity  tunable  lasers  at  1 .55|im,  the 
frequency  of  each  laser  can  be  tuned  in  IGHz  step.  The  optical  signals  are  combined  by  a  3dB  coupler,  and  coupled 
to  the  VMDP  using  a  fiber  pickup  head.  The  microwave  signal  generated  by  optical  mixing  in  the  VMDP  is 
collected  at  the  output  end  of  the  CPS  by  a  50GHz  picoprobe  (GGB  Industries),  which  is  connected  to  an  RF  power 
sensor  and  monitored  by  an  RF  power  meter.  The  calibrated  fiuquency  response  of  long  wavelength  VMDP  is 
shown  in  Figure  7.  At  10  Volt  bias,  a  3dB  bandwidth  fiuquency  of  13GHz  is  measured.  By  scaling  down  the  MSM 
to  deep  sub-micron  scale,  much  higher  bandwidth  (>100  GHz)  is  expected. 
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Figure?.  Measured  frequency  response  of  a  long  wavelength  VMDP 
with  10  MSM  photodiodes  at  different  DC  bias  voltages 

IV.  Conclusion 

In  summary,  we  have  experimentally  demonstrated  the  improved  performance  of  the  long  wavelength  velocity- 
matched  distributed  photodetector  (VMDP).  A  3dB  bandwidth  of  13  GEb  and  a  responsivity  of  0.57  A/W  have  been 
achieved. 

•  This  project  is  supported  by  TRW,  ONR  MURI  on  RF  Photonics,  NRAO,  Tracer  and  HRL  under  UC  MICRO 
andJSEP. 
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Abstract 

A  novel  velocity-matched  distributed  balanced  photodetector  with  a  50Q  coplanar  waveguide 
output  transmission  line  has  been  experimentally  demonstrated  in  the  InP/InGaAs  material  system. 
Distributed  absorption  and  velocity  matching  are  employed  to  achieve  high  saturation  photocurrent.  A 
common  mode  rejection  ratio  greater  than  27  dB  has  been  achieved.  The  RF  link  experiment  conducted  at 
6.48  GHz  shows  that  the  laser  intensity  noise  has  been  suppressed  by  more  than  17  dB. 
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Abstract—  A  novel  velocity-matched  distributed  balanced 
photodetector  with  a  50Q  coplanar  waveguide  output 
transmission  line  has  been  experimentally  demonstrated  in  the 
InP/InGaAs  material  system.  Distributed  absorption  and 
velocity  matching  are  employed  to  achieve  high  saturation 
photocurrent.  A  common  mode  rejection  ratio  greater  than  27 
dB  has  been  achieved.  The  RF  link  experiment  conducted  at 
6.48  GHz  shows  that  the  laser  intensity  noise  has  been 
suppressed  by  more  than  17  dB. 

Index  Terms —  RF  photonics,  microwave  photonics, 
noise  suppression,  analog  fiber  optic  links,  balanced 
photodetectors,  optical  receivers 

L  INTRODUCTION 

Balanced  photodetectors  are  of  great  interest  to 
analog  fiber  optic  links  because  they  can  suppress  laser 
relative  intensity  noise  (RIN)  and  amplified  spontaneous 
emission  noise  (ASE)  from  erbium-doped  fiber  amplifiers 
(EDFA)  [1].  Because  balanced  photodetectors  can  achieve 
shot  noise-limited  link  performance,  we  can  continue  to 
improve  the  noise  figure  and  the  spurious-free  dynamic 
range  (SFDR)  of  externally  modulated  links  by  increasing 
the  power  of  the  optical  carrier.  Therefore,  balanced 
photodetectors  with  broad  bandwidth  and  high  saturation 
photocurrents  are  particularly  important  for  analog  fiber 
optic  link  applications.  Though  discrete  balanced 
photodetectors  with  high  saturation  power  have  been 
reported,  their  bandwidth  is  limited  [2].  Most  of  the 
reported  integrated  balanced  receivers  suffer  from  low 
saturation  power  and  are  not  suitable  for  analog  links  [3-5], 
Previously,  we  have  reported  a  velocity-matched  distributed 
photodetector  (VMDP)  with  a  peak  saturation  photocurrent 
of  56  mA  and  a  3-dB  bandwidth  of  49  GHz  [6].  Recently, 
InP-based  long  wavelength  VMDP  has  also  been  reported 
[7].  Compared  with  other  photodetector  structures,  the 
VMDP  is  more  suitable  for  implementing  balanced 
photodetection  since  it  has  separate  optical  and  microwave 
waveguides.  In  this  paper,  we  propose  and  demonstrate  a 
novel,  monolithic  distributed  balanced  photodetector  that 
can  simultaneously  achieve  high  saturation  photocurrent 
and  large  bandwidth.  A  common  mode  rejection  ratio  of  27 
dB  and  a  noise  suppression  of  17  dB  have  been 
experimentally  demonstrated. 


n.  DESIGN 

Figure  1  shows  the  principle  and  the  schematic 
structure  of  the  distributed  balanced  photodetector.  It 
consists  of  two  input  optical  waveguides,  two  arrays  of 
high-speed  metal-semiconductor-metal  (MSM)  photodiodes 
distributed  along  the  optical  waveguides,  and  a  50Q 
coplanar  waveguide  (CPW)  output  transmission  line.  The 
detector  operates  in  balanced  mode  when  a  voltage  bias  is 
applied  between  the  two  ground  electrodes  of  the  CPW. 
The  common-mode  photocurrent  flows  directly  to  the 
bottom  ground  electrode  while  the  difference  photocurrent 
(signal)  flows  to  the  center  conductor.  The  signal  is  then 
collected  by  the  CPW.  The  diodes  are  23  pm  long  and  5  pm 
wide.  The  separation  between  photodiodes  is  150  pm.  The 
MSM  fingers  with  1pm  width  and  1pm  spacing  are 
patterned  by  optical  lithography.  The  overlap  of  the  MSM 
fingers  is  10.5  pm.  The  central  conductor  of  the  CPW  has  a 
width  of  55  pm  and  the  separation  between  the  central 
conductor  and  the  ground  conductors  is  85  pm.  Without  the 
photodiodes,  the  velocity  of  the  CPW  is  about  31.8%  faster 
than  the  light  velocity  in  the  optical  waveguide.  The 
photodiode  arrays  provide  the  periodic  capacitance  loading 
to  slow  down  the  microwave  velocity.  By  adjusting  the 
length  and  separation  of  photodiodes,  velocity  matching 
between  the  CPW  and  the  optical  waveguides  is  achieved. 
The  impedance  of  CPW  is  also  matched  to  50Q. 

The  distributed  balanced  photodetector  inherits  the 
basic  advantages  of  the  VMDP,  namely,  high  saturation 
photocurrent,  high  quantum  efficiency,  and  large 
bandwidth.  It  should  be  noted  that  even  though  only  the 
difference  current  (AC  signal)  is  collected  in  the  balanced 
photodetector,  the  DC  light  is  still  absorbed  in  the 
photodiodes.  As  a  result,  high  DC  saturation  photocurrent  is 
required  for  the  distributed  balanced  photodetectors.  By 
coupling  only  a  small  fraction  of  light  from  the  passive 
waveguide  to  each  individual  photodiodes,  the  photodiodes 
are  kept  below  saturation  even  under  intense  optical  input. 
Though  longer  absorption  length  is  required,  the  bandwidth 
of  the  distributed  balanced  photodetector  remains  high 
because  of  the  velocity  matching.  The  linearity  of  the 
detector  is  also  improved  by  distributed  absorption  because 
the  photo  carrier  density  is  reduced  in  the  active  region. 

Optical  Waveguide 

The  optical  waveguide  consists  of  the  following:  a 
200-nm-thick  Ino.52Alo.37Gao.11  As  lower  cladding  layer,  a 
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500-nm-thick  Ino.52Alo.178Gao.302 As  core  region,  a  200-nni- 
thick  Ino.52Alo.37Gao.11  As  first  upper  cladding  layer,  and  a 
thin  Ino.52Alo.48As  second  upper  cladding  layer.  The  150- 
nm-thick  absorption  region  is  located  on  top  of  the 
waveguide  for  evanescent  coupling.  Since  the  Schottky 
barrier  height  of  most  metals  on  InGaAs  is  typically 
between  0.2-0.3  eV,  an  Ino.52Alo.48As  cap  layer  is  used  to 
increase  the  Schottky  barrier  height  and  therefore  decrease 
the  dark  current  of  the  photodiodes  [8,9].  A  graded  layer  is 
incorporated  in  the  structure  to  reduce  the  minority  carrier 
trapping  at  the  InAlAs-InGaAs  band  edge  discontinuity. 

A  scalar  three-dimensional  beam  propagation 
method  (BPM)  was  used  to  simulate  the  optical  properties 
of  the  balanced  VMDP.  The  absorption  per  photodiode  and 
the  optical  coupling  loss  between  photodiodes  are  estimated 
to  be  8.8%  and  3.2%,  respectively.  Our  simulation  shows 
that  an  AC  quantum  efficiency  of  42%  can  be  achieved  for 
balanced  receiver  with  10  pairs  of  photodiodes.  Since  the 
two  optical  waveguides  are  140  pm  apart,  no  optical 
coupling  between  the  waveguides  is  expected.  This  is 
confirmed  by  the  BPM  simulation.  Figure  2  shows  the 
optical  intensity  profile  and  the  power  distribution  along  the 
waveguides  when  a  fundamental  mode  is  launched  into  the 
first  waveguide  (WGl)  at  the  input.  The  waveguide  powers 
were  calculated  using  computational  windows  surrounding 
the  individual  waveguides.  Indeed,  no  power  coupling  to 
WG2  is  observed.  The  slight  decay  in  WGl  (0.195 
dB/mm)  is  attributed  to  loss  to  slab  mode. 

Modeling  of  Microwave  Transmission  Line 

The  impedance  and  the  phase  velocity  of  the  CPW 
are  calculated  using  the  equivalent  circuit  model  described 
in  Ref.  [10].  The  length  of  the  photodiodes  and  the 
separation  between  them  are  adjusted  to  achieve 
simultaneous  velocity  matching  and  impedance  matching. 
Since  the  separation  between  the  central  conductor  and  the 
ground  electrodes  (85  pm)  are  much  smaller  than  the 
wavelength  of  the  RF  signal  (about  two  order  of  magnitude 
smaller  at  50  GHz),  quasi  static  analysis  is  reasonably 
accurate  [11].  The  capacitors  and  resistors  of  each 
photodiodes  are  considered  as  lumped  elements  in  our 
quasi-static  simulation.  After  optimizing  the  receiver 
structure,  a  full-wave  analysis  was  performed  to  verify  the 
design.  We  found  that  the  quasi-static  results  agree  very 
well  with  the  full  wave  analysis  for  frequency  below  100 
GHz.  The  period  in  our  device  corresponds  to  a  cut-off 
frequency  of  300  GHz,  well  above  our  expected  frequency 
of  operation.  Therefore,  the  dispersion  due  to  the 
periodicity  of  the  structure  is  negligible  for  frequencies 
below  100  GHz  [12], 

The  overall  performance  of  the  balanced  VMDP  is 
simulated  using  the  model  in  Ref.  [13],  which  includes  the 
transit-time  frequency  response  of  the  photodiodes,  the  loss 
and  dispersion  of  the  microwave  transmission  line,  the 
optical  coupling  loss  between  the  photodiodes,  and  the 
residue  velocity  mismatch  between  the  CPW  and  the  optical 
waveguides. 


in.  FABRICATION 

Fabrication  of  the  receiver  started  with  removing 
the  InGaAs  layer  except  in  the  active  areas  of  the 
photodiodes.  Ridge  waveguides  with  100  nm  ridge  height 
were  formed  by  wet  chemical  etching.  The  active  regions  of 
the  photodiodes  were  defined  by  opening  6x23-pm^ 
windows  on  the  150-nm-thick  silicon  nitride  (Si3N4)  film 
deposited  by  plasma-enhanced  chemical  vapor  deposition 
(PECVD).  Buffer  HF  was  used  to  open  the  windows.  The 
Ti-Au  electrodes  and  contact  pads  were  then  delineated  by 
standard  lift-off  processes.  The  tips  of  the  MSM  fingers  are 
placed  on  top  of  the  Si3N4  to  suppress  the  soft  breakdown 
and  enable  the  MSM  diodes  to  operate  over  a  wider  range 
of  bias  voltages  [14].  A  350-nm-thick  coplanar  waveguide 
was  formed  by  standard  lift-off  process  to  connect  the 
distributed  balanced  photodetectors.  Finally  the  balanced 
detector  structure  is  cleaved  and  mounted  on  copper  heat 
sinks.  By  measuring  the  forward  current-voltage 
characteristics,  the  barrier  height  of  the  metal- 
semiconductor  junction  was  estimated  to  be  0.57  eV. 

IV.  EXPERIMENTAL  RESULTS 

The  balanced  VMDP  exhibits  very  good  electrical 
and  optical  characteristics.  The  dark  current  is  measured  to 
be  28  pA/cm^  at  10  V  bias;  the  lowest  reported  for 
InAlAs/InGaAs  MSM  photodiodes  (Figure  3).  At  the 
operating  voltage  of  4  V,  the  total  dark  current  of  a 
balanced  receiver  with  5  pairs  of  photodiodes  is  1.5  nA.  We 
used  a  pair  of  lensed  fibers  to  couple  light  into  the 
photodetector.  Figure  3  also  shows  the  measured 
responsivity  of  the  photodetector  as  a  function  of  bias 
voltage.  The  average  DC  responsivity  was  measured  to  be 
0.45  AAV  at  8  V  bias.  Responsivity  as  high  as  0.6  AAV  has 
been  observed  in  some  devices.  With  anti-reflection 
coating,  the  average  responsivity  can  be  increased  to  0.64 
AAV.  The  DC  photocurrent  at  8  V  versus  the  input  optical 
power  is  shown  in  Figure  4.  The  photocurrent  of  each 
VMDP  remains  linear  up  to  12  mA  (27  mW  of  optical 
power).  We  did  not  increase  the  optical  power  further  to 
avoid  damaging  the  devices. 

An  HP  85 IOC  network  analyzer  was  used  to 
measure  the  characteristic  impedance  and  the  microwave 
return  loss  (Sn)  of  the  balanced  receiver.  Figure  5  shows 
the  measured  S  parameters  of  the  receiver.  The  Sn  is  below 
-30  dB  from  45  MHz  to  40  GHz,  indicating  that  the 
impedance  of  the  detector  is  very  well  matched  to  50  Q. 
The  loaded  CPW  has  very  low  insertion  loss.  The  measured 
S12  shows  a  drop  of  only  0.6  dB  from  45  MHz  to  40  GHz. 

The  frequency  response  of  the  balanced  VMDP 
was  first  characterized  with  light  coupled  to  one  waveguide 
only.  Figure  6  shows  the  fi-equency  response  of  the 
photodetector.  Using  the  optical  heterodyne  technique  with 
two  external  cavity  tunable  semiconductor  lasers  at  1.55 
pm,  the  3-dB  bandwidth  was  found  to  be  16  GHz  for  both 
photodetector  arrays.  The  bandwidth  is  currently  limited  by 
the  carrier  transit  time  of  the  MSM  photodiodes.  Since  the 
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bandwidth  of  our  capacitance  loaded  CPW  is  much  greater 
than  40  GHz,  the  bandwidth  of  the  balanced  VMDP  can  be 
increased  by  scaling  down  the  MSM  photodiodes. 
Theoretical  simulation  indicates  that  bandwidth  >100  GHz 
is  achievable. 


Figure  7  depicts  the  experimental  setup  for 
balanced  detection.  A  distributed  feedback  (DFB)  laser  with 
1542-nm  wavelength  and  0  dBm  output  power  is  employed 
as  optical  source.  It  is  amplified  by  an  EDFA  and  then 
filtered  by  an  optical  bandpass  filter  with  2-nm  bandwidth. 
The  microwave  signal  was  modulated  onto  the  optical 
carrier  by  an  X-coupled  Mach-Zehnder  modulator  (MZM), 
which  produces  two  complimentary  outputs  for  the 
balanced  VMDP.  The  outputs  are  coupled  to  the  balanced 
VMDP  by  two  lensed  fibers.  To  maximize  the  signal 
enhancement  and  noise  cancellation,  it  is  important  to 
match  the  amplitudes  and  phases  of  the  two  detected 
microwave  signals.  In  our  experiment,  a  variable  attenuator 
was  used  to  match  the  amplitudes  of  the  photocurrents. 
Typical  balance  between  the  VMDP  was  within  2%  of  the 
total  photocurrent  reading.  We  also  employed  a  variable 
optical  delay  line  to  match  the  microwave  phases  of  the 
detected  signals.  For  balanced  detection,  the  RF  signals 
need  to  be  exactly  180°  out  of  phase.  The  maximum  delay 
in  our  setup  is  50  mm,  which  allows  us  to  match  the  phases 
of  RF  signals  above  6  GHz. 


Balanced  detection  is  achieved  by  applying  a  bias 
of  8  V  between  the  two  ground  electrodes  of  the  CPW.  A 
custom-made  high-frequency  probe  with  an  integrated  De¬ 
blocking  capacitor  on  one  ground  probe  is  used  to  collect 
the  microwave  output  signal.  Figure  8  shows  the  biasing 
scheme.  The  DC  photocurrents  are  measured  through  two 
ammeters:  the  difference  photocurrent  {Idiff)  is  monitored 
through  the  bias-T  connected  to  the  probe,  and  the  common 
mode  photocurrent  {icohd  is  monitored  between  the  ground 
electrodes  of  the  CPW.  Figure  9  shows  the  common  mode 
rejection  ratio  (CMRR),  defined 
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photocurrent.  Very  high  CMRR  (>  27  dB)  is  measured  for  a 
wide  range  of  photocurrent  from  30  nA  to  12  mA.  This  is 
attributed  to  the  well-matched  characteristics  of  the 
photodiodes  in  our  monolithic  balanced  detectors. 


To  verify  the  balanced  detection,  we  modulated 
the  optical  input  at  8  GHz  and  varied  the  delay  of  the  fiber 
optic  delay  line  from  0  to  50  mm.  Figure  10  shows  the 
detected  signal  versus  the  delay.  When  the  received  RF 
signals  have  a  0°  phase  difference,  the  RF  output  is 
cancelled  out  in  balanced  mode.  When  the  detected  RF 
signals  are  180®  out  of  phase,  the  small  signal  photocurrents 
from  the  two  photodiode  arrays  are  added  in  phase  and 
maximum  RF  signal  is  observed.  The  extinction  ratio  is 
more  than  44  dB. 


cancellation  ratio  of  our  device,  we  compare  the  noise 
spectra  of  a  DFB  laser  measured  by  our  photodetector  in  the 
unbalanced  and  the  balanced  modes.  The  DFB  laser 
operates  in  CW  condition  and  has  a  RIN  peak  at  6.5  GHz. 
The  top  trace  in  Figure  11  shows  the  noise  spectra 
measured  near  its  RIN  peak  when  only  one  waveguide  is 
illuminated.  We  confirm  that  the  noise  is  dominated  by  the 
RIN  of  the  DFB  laser.  When  the  optical  input  power  is 
doubled,  the  noise  floor  increases  by  6  dB.  We  then  biased 
the  photodetector  in  the  balanced  mode  and  coupled  the 
input  to  both  the  optical  waveguides.  The  bottom  trace  in 
Fig.  10  shows  the  noise  spectra  detected  by  the  balanced 
detector.  Suppression  of  the  RIN  by  as  much  as  15  dB  is 
observed.  The  maximum  suppression  occurs  at  6.48  GHz, 
at  which  frequency  the  phases  of  the  detected  noise  matches 
exactly.  Since  optical  path  lengths  of  the  two  input  fibers 
are  not  equal  (due  to  the  insertion  of  the  fiber  delay  line), 
the  cancellation  of  the  laser  RIN  is  not  as  effective  when 
the  frequency  deviates  from  6.48  GHz.  With  equal  fiber 
lengths  for  both  inputs,  uniform  noise  cancellation  is 
expected  for  all  the  frequencies. 

Figure  12  shows  the  RF  spectra  of  the  output  from 
the  balanced  VMDP  in  the  unbalanced  (only  one  waveguide 
is  illuminated)  and  the  balanced  mode.  Suppression  of  the 
noise  floor  by  17  dB  has  been  observed  in  the  balanced 
mode  over  a  wide  frequency  range  from  6  to  15  GHz.  This 
is  equivalent  to  23  dB  improvement  of  noise  if  two  outputs 
were  considered.  The  signal  is  also  enhanced  by  6  dB.  For 
frequencies  below  6  GHz,  the  fiber  length  for  one  RF 
period  is  longer  than  50  mm,  exceeding  the  maximum  delay 
length  in  our  current  setup.  We  are  currently  working  on  a 
new  setup  with  two  fiber  optic  delay  lines  to  achieve 
exactly  equal  optical  path  lengths  for  broad  band  operation 
without  any  tuning  of  the  delay  line.  This  will  allow  us  to 
measure  the  AC  responsivity  of  the  device  in  the  balanced 
mode. 


V.  CONCLUSION 

We  have  successfully  designed,  fabricated,  and 
experimentally  demonstrated  a  balanced  velocity-matched 
distributed  photodetector  (VMDP)  with  both  impedance 
and  velocity  matching.  The  device  exhibits  a  very  low  dark 
current  (1.5  nA  for  the  balanced  VMDP  with  10  pairs  of 
photodiodes)  and  a  high  external  quantum  efficiency  (0.64 
AAV).  The  relative  intensity  noise  (RIN)  of  a 
semiconductor  distributed  feedback  laser  has  been 
suppressed  by  17  dB,  and  the  RF  signal  has  been  enhanced 
by  6  dB.  To  the  authors’  knowledge,  this  is  the  first  report 
of  monolithic  distributed  balanced  photodetectors.  The 
experimental  results  indicate  that  the  distributed  balanced 
photodetector  will  have  a  major  impact  on  most  RF 
photonic  systems. 


One  key  advantage  of  the  balanced  photodetector 
is  its  ability  to  cancel  out  the  laser  RIN.  To  evaluate  the 
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Figure  1.  (a)  Principle  and  (b)  schematic  structure  of  the  distributed  balanced  photodetector.  The  inset  shows  the 
active  region  with  an  MSM  photodiode. 
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Figure  2.  (a)  Simulated  optical  field  profile  and  (b) 
optical  power  distribution  along  the  waveguides. 


Figure  4.  DC  photocurrent  versus  input  optical 
power  (without  AR  coating).  The  DC 
photocurrent  remains  linear  up  to  12  mA. 
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Figure  3.  The  dark  current  (left)  and  the  responsivity 
(right)  versus  bias  voltage  for  individual  VMDP  (one 
array  of  photodiodes)  in  the  balanced  detector.  At  the 
operational  bias  (4  V),  the  dark  current  is  0. 15  nA. 


Figure  5.  The  Measured  S  parameters  of  the  CPW. 


Figure  6,  Frequency  response  of  the  distributed 
balanced  photodetector  with  only  one  waveguide 
illuminated. 
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Figure  7.  Experimental  setup  for  balanced  detection.  The 
complimentary  input  signals  are  produced  by  the  X 
coupled  MZ  modulator. 
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Figure  9.  Common  mode  rejection  ratio  (CMRR) 
versus  photocurrent  for  the  distributed  balanced 
receiver.  The  high  CMRR  results  from  the  closely 
matched  photodiode  characteristics  in  our  monolithic 
device. 


Capacitor  for 
Balanced  Biasing 


Figure  8.  The  biasing  circuit  for  balanced  detection. 


Delay  Line  Reading  (mm) 

Figure  10.  The  peak  RF  power  detected  in  the 
balanced  mode  versus  the  time  delay  of  one  input  light. 
Periodic  variation  is  observed.  The  period  is  measured 
to  be  38.75  mm  at  8  GHz  RF  signal,  which  is  very 
close  to  the  theoretical  value.  The  Extinction  ratio  is 
greater  than  44  dB. 
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Figure  11.  The  noise  spectra  of  a  DFB  laser  measured 
by  the  distributed  balanced  photodetector  in 
unbalanced  mode  (top  trace)  and  balanced  mode 
(bottom  trace).  Maximum  noise  cancellation  occurs  at 
the  frequency  at  which  the  RF  noise  is  180°  out  of 
phase. 
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Figure  12.  The  RF  spectra  of  the  modulated  signal  detected  in  (a)  unbalanced  and  (b)  balanced  mode.  The  peak  of  the 
balanced  output  is  6-dB  higher  than  that  of  the  unbalanced  signal.  Noise  suppression  in  excess  of  17-dB  is  achieved. 
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ABSTRACT 

We  report  on  the  noise  suppression  properties  of  a  novel  distributed  balanced  photodetector  that  we 
designed  and  fabricated.  Large  improvement  in  signal-to-noise  ratio  (SNR)  has  been  observed  over  a 
wide  range  of  received  optical  power  and  RF  signal  phase  mismatch. 
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Balanced  photodetectors  are  of  great  interest  to  analog  fiber  optic  links  because  they  can  suppress 
relative  intensity  noise  (RIN)  of  lasers  and  amplified  spontaneous  emission  noise  (ASE)  from  erbium- 
doped  fiber  amplifiers  (EDFA).  When  used  in  conjunction  with  an  external  modulator  with 
complementary  outputs,  shot  noise-limited  system  performance  can  be  achieved  [1].  In  order  to  fully 
exploit  the  advantages  of  the  balanced  systems,  balanced  photodetectors  with  high  saturation 
photocurrents  and  broad  bandwidth  are  required.  Monolithically  integrated  balanced  photodetectors  offer 
superior  performance  (broader  bandwidth,  better  matching  of  photodiodes)  and  reduced  packaging  cost. 
However,  most  of  the  integrated  balanced  photodetectors  reported  to  date  have  low  saturation 
photocurrents  and  are  not  suitable  for  analog  applications. 

Previously,  we  have  reported  a  novel  distributed  balanced  photodetector  [2].  It  integrates  two 
velocity  matched  distributed  photodetectors  (VMDP)  with  a  microwave  coplanar  waveguide  (CPW).  This 
new  device  inherits  the  basic  advantages  of  the  VMDP,  namely,  broad  bandwidth  and  high  saturation 
photocurrent  [3].  In  this  paper,  we  reported  on  the  noise  suppression  properties  of  the  monolithic  balanced 
distributed  photodetector.  A  high  signal-to-noise  (SNR)  and  large  suppression  of  the  laser  RIN  over  a 
broad  range  of  input  powers  has  been  observed  in  the  RF  photonic  link  employing  the  balanced  VMDP. 

Figure  1  shows  the  schematic  of  the  balanced  distributed  photodetector.  It  consists  of  two  input 
optical  waveguides,  two  arrays  of  high-speed  metal-semiconductor-metal  (MSM)  photodiodes  distributed 
along  the  optical  waveguides,  and  a  50Q  coplanar  waveguide  (CPW)  output  transmission  line.  By  biasing 
the  photodiodes  in  balanced  mode,  only  the  difference  photocurrents  are  collected  by  the  CPW.  The 
balanced  VMDP  showed  a  RIN  suppression  of  more  than  17  dB  with  a  very  good  common  mode 
rejection  ratio  (CMRR)  >  27  dB. 

An  RF  photonic  link  with  a  dual-output  complementary  Mach-Zehnder  modulator  (MZM)  has 
been  set  up  to  characterize  the  noise  properties  of  the  distributed  balanced  receiver.  The  RF  carrier  is 
fixed  at  6.5  GHz.  Figure  2  shows  the  SNR  of  the  link  versus  the  received  optical  power  for  both  the 
distributed  balanced  receiver  and  the  reference  receiver  with  a  single  detector.  The  SNR  for  the  single 
detector  receiver  is  almost  constant  with  increasing  optical  power,  indicating  that  the  receiver  noise  is 
dominated  by  the  RIN.  In  contrast,  the  SNR  for  the  balanced  receiver  increases  monotonically  with 
optical  power  upto  15  mW.  For  optical  power  greater  than  15  mW,  the  SNR  still  increases  monotonically 
with  optical  power  at  a  smaller  rate.  This  suggests  the  presence  of  some  residue  RIN.  Comparing  the  two 
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receivers,  it  is  also  noted  that  the  SNR  for  the  balanced  receiver  is  23  dB  higher  than  the  single  detector 
receiver. 

The  RF  signals  detected  by  the  balanced  detector  should  be  exactly  at  180“  out  of  phase  for 
optimum  noise  suppression.  This  requires  the  lengths  of  the  fiber  from  the  complimentary  MZM  to  the 
detector  to  be  exactly  the  same  m  length.  In  practical  applications,  the  optical  path  length  will  drift 
slightly  due  to  environmental  changes  and  it  is  important  to  imderstand  the  impact  of  phase  mismatch. 
Figure  3  plots  the  variation  in  the  signal  peak  and  the  amoimt  of  noise  suppression  when  the  phase 
difference  of  the  RF  signal  at  6.5  GHz  deviates  from  180°.  The  measured  data  shows  that  even  with  a 
phase  variation  of  100°,  our  device  can  suppress  more  than  9  dB  of  RIN  noise  whereas  the  signal  peak 
reduction  was  less  than  4  dB.  As  balanced  detector  amplifies  the  RF  signal  by  6dB,  with  a  phase 
deviation  of  100°,  the  signal  is  still  higher  than  that  of  an  unbalanced  detector.  For  the  RF  signal  at  6.5 
GHz,  100°  phase  variation  between  two  fibers  corresponds  to  ~8.5nim  of  fiber  length. 

In  conclusion,  we  measured  the  SNR  and  phase  tolerance  of  the  input  RF  signals  of  the  high 
power  distributed  balanced  receiver.  This  monolithic  balanced  receiver  has  the  potential  to  dramatically 
improve  the  performance  of  RF  photonic  links. 
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Figure  1.  Principle  and  schematic  structure  of  the  distributed  balanced  photodetector.  We 
cascaded  multiple  balanced  photodetector  pairs  in  series  along  a  coplanar  waveguide  to  increase 
saturation  photocurrent. 
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Figure  2.  Measured  signal-to-noise  ratio  (SNR)  of  the  balanced  receiver.  The  lower  curve  shows 
the  SNR  for  the  receiver  in  unbalanced  configuration. 
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Phase  Difference  (Degrees) 


Figure  3.  Total  amount  of  noise  suppression  and  the  deviation  of  the  signal  peak  versus  phase 
deviation  of  the  RF  signal  from  180°.  Even  with  a  phase  variation  of  100°,  the  receiver  can 
suppress  more  than  9  dB  of  RIN  noise  whereas  the  signal  peak  reduction  was  less  than  4dB. 
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ABSTRACT 

A  novel  velocity-matched  distributed  balanced  photodetector  with  a  50Q  coplanar 
waveguide  output  transmission  line  has  been  experimentally  demonstrated  in  the  InP/InGaAs 
material  system.  Distributed  absorption  and  velocity  matching  are  employed  to  achieve  high 
saturation  photocurrent.  A  common  mode  rejection  ratio  of  27  dB  has  been  achieved.  The  RF 
link  experiment  conducted  at  6.48  GHz  shows  that  the  laser  intensity  noise  has  been  suppressed 
by  more  than  17  dB. 
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I.  INTRODUCTION 

Balanced  photodetectors  are  of  great  interest  to  analog  fiber  optic  links  because  they  can 
suppress  laser  relative  intensity  noise  (RIN)  and  amplified  spontaneous  emission  noise  (ASE) 
firom  erbium-doped  fiber  amplifiers  (EDFA)  [1].  Because  balanced  photodetectors  can  achieve 
shot  noise-limited  link  performance,  we  can  continue  to  improve  the  noise  figure  and  the 
spurious-firee  dynamic  range  (SFDR)  of  externally  modulated  links  by  increasing  the  power  of 
the  optical  carrier.  Therefore,  balanced  photodetectors  with  broad  bandwidth  and  high  saturation 
photocurrents  are  particularly  important  for  analog  fiber  optic  link  applications.  Though  discrete 
balanced  photodetectors  with  high  saturation  power  have  been  reported,  their  bandwidth  is 
limited  [2].  Most  of  the  reported  integrated  balanced  receivers  suffer  from  low  saturation  power 
and  are  not  suitable  for  analog  links.  [3-5]  Previously,  we  have  reported  a  velocity-matched 
distributed  photodetector  (VMDP)  with  a  peak  saturation  photocurrent  of  56  mA  and  a  3-dB 
bandwidth  of  49  GHz  [6].  Recently,  InP -based  long  wavelength  VMDP  has  also  been  reported 
[7].  Compared  with  other  photodetector  structures,  the  VMDP  is  more  suitable  for  implementing 
the  balanced  photodetection  since  it  has  separate  optical  and  microwave  waveguides.  In  this 
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letter,  we  propose  and  demonstrate  a  novel,  monolithic  distributed  balanced  photodetector  that 
can  simultaneously  achieve  high  saturation  photocurrent  and  large  bandwidth.  A  common  mode 
rejection  ratio  of  27  dB  and  a  noise  suppression  of  17dB  have  been  experimentally  demonstrated. 


II.  DESIGN  AND  FABRICATION 

Figure  1  shows  the  schematic  of  the  distributed  balanced  photodetector.  It  consists  of  two 
input  optical  waveguides,  two  arrays  of  high-speed  metal-semiconductor-metal  (MSM) 
photodiodes  distributed  along  the  optical  waveguides,  and  a  500  coplanar  waveguide  (CPW) 
ou^ut  transmission  line.  The  diodes  are  23  pm  long  and  5  pm  wide.  The  MSM  fingers  with  1  pm 
width  and  1pm  spacing  are  patterned  by  optical  lithography.  The  central  conductor  of  the  CPW 
has  a  width  of  55  pm  and  the  separation  between  the  central  conductor  and  the  ground 
conductors  is  85  pm.  The  10.5  pm  overlap  of  the  MSM  fingers  is  designed  to  provide  the 
required  capacitive  loading  for  velocity  matching.  Unlike  previously  reported  slow-wave  CPW 
[8]  that  ignored  the  resistance  of  the  MSM  photodiode  fingers,  finite  metal  thickness  and 
transmission  line  discontinuities  in  their  quasi-static  calculations,  we  used  a  full-wave  analysis  to 
achieve  the  velocity  and  impedance  matching. 

The  distributed  balanced  photodetector  inherits  the  basic  advantages  of  the  VMDP, 
namely,  high  saturation  photocurrent,  high  quantum  efficiency  and  large  bandwidth.  It  should  be 
noted  that  even  though  only  the  difference  current  (AC  signal)  is  collected  in  the  balanced 
photodetector,  the  DC  light  is  still  absorbed  in  the  photodiodes.  As  a  result,  high  DC  saturation 
photocurrent  is  required  for  the  distributed  balanced  photodetectors.  By  coupling  only  a  small 
firaction  of  light  from  the  passive  waveguide  to  each  individual  photodiode,  the  photodiodes  are 
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kept  below  saturation  even  under  intense  optical  input.  The  bandwidth  of  the  distributed 
balanced  photodetector  remains  high  because  of  the  velocity  matching. 

Beam  propagation  method  (BPM)  simulation  was  used  to  investigate  the  optical 
properties  of  the  balanced  VMDP.  The  simulation  results  indicate  that  there  is  no  crosstalk 
between  the  parallel  optical  waveguides  that  are  140  pm  apart.  The  simulation  also  shows  that  an 
input  beam  with  TM  polarization  contributes  to  the  optimum  performance  and  the  absorption  per 
photodiode  and  the  optical  coupling  loss  between  photodiodes  are  shown  to  be  8.8%  and  3.2%, 
respectively  for  the  TM  polarization.  Our  calculation  shows  that  the  AC  quantum  efficiency  of 
42%  can  be  achieved  for  balanced  receiver  with  10  pairs  of  photodiodes.  The  optical  waveguide 
consists  of  the  following:  a  200-nm-thick  Ino.52Alo.37Gao.11  As  lower  cladding  layer,  a  500-nm- 
thick  Ino.52Alo.i78Gao.302As  core  region,  a  200-nm-thick  Ino.52Alo.37Gao.11  As  first  upper  cladding 
layer,  and  a  thin  Ino.52Alo.48 As  second  upper  cladding  layer.  The  150-nm-thick  absorption  region 
is  located  on  top  of  the  waveguide  for  evanescent  coupling.  Since  the  Schottky  barrier  height  of 
most  metals  on  the  InGaAs  is  typically  between  0.2-0.3  eV,  an  Ino.52Alo.48As  cap  layer  is  used  to 
increase  the  Schottky  barrier  height  and  therefore  decrease  the  dark  current  in  the  photodiodes 
[9,10].  A  graded  layer  is  incorporated  in  the  structure  to  reduce  the  minority  carrier  trapping  at 
the  InAlAs-InGaAs  band  edge  discontinuity. 

Fabrication  of  the  receiver  started  with  removing  the  InGaAs  layer  except  in  the  active 
areas  of  the  photodiodes.  Ridge  waveguides  with  100  nm  ridge  height  were  formed  by  wet 
chemical  etching.  The  active  regions  of  the  photodiodes  were  defined  by  opening  6pm  x23  pm 
windows  on  the  150-nm-thick  silicon  nitride  (Si3N4)  film  deposited  by  plasma-enhanced 
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chemical  vapor  deposition  (PECVD).  Buffered  HF  was  used  to  remove  the  Si3N4  in  the  windows. 
The  Ti-Au  electrodes  and  contact  pads  of  the  MSM  photodiodes  were  then  delineated  by 
standard  lift-off  processes.  The  tips  of  the  MSM  fingers  are  placed  on  top  of  Si3N4  to  suppress 
the  soft  breakdown  and  enables  the  MSM  diodes  to  operate  over  a  wider  range  of  bias  voltage 
[1 1].  A  350-nm-thick  coplanar  waveguide  was  formed  by  standard  lift-off  process  to  connect  the 
distributed  balanced  photodetectors.  By  measuring  the  forward  current-voltage  characteristics, 
the  barrier  height  of  the  metal-semiconductor  junction  was  estimated  to  be  0.57  eV. 


III.  MEASUREMENTS  AND  DISCUSSIONS 
The  balanced  VMDP  exhibits  very  good  electrical  and  optical  characteristics.  The  dark 
current  is  measured  to  be  28pA/cm^  at  lOV  bias,  the  lowest  reported  for  InAlAs/InGaAs  MSM 
photodiodes.  By  coupling  light  directly  fi*om  a  lensed  fiber,  the  average  DC  responsivity  (15 
devices  tested)  was  measured  to  be  0.45AAV  at  8V  bias.  Responsivity  as  high  as  0.6AAV  has 
been  observed  in  some  devices.  The  photo  response  of  a  laser  beam  with  TM  polarization  is 
measured  to  be  ~3  dB  higher  than  that  of  TE  polarization.  The  coupling  efficiency  of  the  lensed 
fibers  in  our  setup  was  calculated  to  be  ~50  %.  With  antirefiection  (AR)  coating  to  the  balanced 
receiver  facet  (30%  Fresnel  loss),  the  average  responsivity  can  be  increased  to  0.64 A/W.  The 
responsivity  can  be  further  improved  by  optimizing  the  coupling  efficiency  of  the  lensed  fiber  as 
well  as  reducing  the  coupling  loss  between  the  passive  and  active  waveguide  regions  by  better 
controlling  the  etching  steps  in  the  fabrication. 

An  HP  85 IOC  network  analyzer  was  used  to  measure  the  characteristic  impedance  and  the 
microwave  return  loss  (Sn)  of  the  balanced  receiver.  The  characteristic  impedance  of  the 
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receiver  is  very  well  matched  to  50Q.  The  Sn  is  as  low  as  -30  dB  from  45MHz  to  40GHz.  The 
CPW  has  very  low  insertion  loss.  The  measured  S12  shows  a  drop  of  only  0.6dB  from  45MHz  to 
40  GHz. 

The  frequency  response  of  the  balanced  VMDP  was  first  characterized  with  light  coupled 
to  one  waveguide  only.  Using  the  optical  heterodyne  technique  with  two  external  cavity  tunable 
semiconductor  lasers  at  1.55  pm,  the  3-dB  bandwidth  was  found  to  be  16  GHz  for  both 
photodetector  arrays.  The  bandwidth  is  currently  limited  by  the  carrier  transit  time  of  the  MSM 
photodiodes.  Since  our  bandwidth  of  the  capacitance  loaded  CPW  is  much  greater  than  40  GHz, 
the  bandwidth  of  the  balanced  VMDP  can  be  increased  by  scaling  down  the  MSM  photodiodes. 
Theoretical  simulation  indicates  that  bandwidth  >100  GHz  is  achievable. 

Figure  2(a)  depicts  the  experimental  setup  for  balanced  detection.  A  distributed  feedback 
(DFB)  laser  with  1542-nm  wavelength  and  0  dBm  output  power  is  employed  as  optical  source.  It 
is  amplified  by  an  EDFA  and  then  filtered  by  an  optical  bandpass  filter  with  2-nm  bandwidth. 
The  microwave  signal  was  modulated  onto  the  optical  carrier  by  an  X-coupled  Mach-Zehnder 
modulator  (MZM),  which  produces  two  complementary  ouqjuts  for  the  balanced  VMDP.  The 
outputs  are  coupled  to  the  balanced  VMDP  by  two  lensed  fibers.  To  maximize  the  signal 
enhancement  and  noise  cancellation,  it  is  important  to  match  the  amplitudes  and  phases  of  the 
two  detected  microwave  signals.  In  our  experiment,  a  variable  attenuator  was  used  to  match  the 
amplitudes  of  the  photocurrents.  Typical  balance  between  the  VMDP  was  within  2%  of  the  total 
photocurrent  reading,  which  was  limited  to  4  mA  in  our  current  experiments.  Since  the  accuracy 
of  cutting  optical  fiber  lengths  is  about  1  mm,  a  variable  delay  line  was  used  to  ensure  equal 
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optical  path  lengths  between  the  two  fibers,  which  is  necessary  to  produce  the  broadband  180" 
out-of-phase  RF  signals  for  balanced  VMDP.  The  output  from  the  balanced  VMDP  was 
measured  using  a  spectrum  analyzer.  A  common  mode  rejection  ratio  of  -27dB  has  been 
achieved. 

Balanced  detection  is  achieved  by  applying  a  bias  of  8V  between  the  two  ground 
electrodes  of  the  CPW.  A  custom-made  high-frequency  probe  with  an  integrated  DC-bloeking 
eapacitor  on  one  ground  probe  is  used  to  collect  the  microwave  output  signal.  The  DC  difference 
photocurrent  is  monitored  through  the  bias-T  connected  to  the  probe.  To  verify  the  balanced 
detection,  we  modulated  the  optieal  input  at  8-GHz  and  varied  the  delay  of  the  fiber  optic  delay 
line  from  0  to  50  mm.  Figure  2(b)  shows  the  RF  signal  versus  the  delay.  When  the  RF  signals 
have  0"  phase  difference,  the  AC  output  is  cancelled  in  balanced  mode.  At  180°  phase  difference, 
the  RF  signal  is  amplified,  the  extinction  ration  is  more  than  44  dB. 

Figure  3  shows  the  RF  speetrum  of  the  output  from  the  balanced  VMDP  in  the 
unbalanced  (  only  one  waveguide  is  illuminated  )  and  the  balanced  mode.  Suppression  of  the 
noise  floor  by  17-dB  has  been  observed  in  the  balanced  mode  over  a  wide  frequeney  range  from 
6  to  15  GHz.  The  signal  is  also  enhanced  by  6-dB.  For  bandwidths  lower  than  6  GHz,  the  fiber 
length  for  one  RF  period  is  longer  than  50  mm,  the  maximum  delay  length  in  our  eurrent  setup. 
We  are  currently  working  on  a  new  setup  with  two  fiber  optic  delay  lines  to  aehieve  exactly 
equal  optical  path  lengths  for  broad  band  operation  without  any  tuning  of  the  delay  line.  This 
will  allow  us  to  measure  the  broadband  charaeteristics  of  the  device  in  the  balanced  mode. 
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IV.  CONCLUSION 


We  have  successfully  designed,  fabricated,  and  experimentally  demonstrated  a  balanced 
velocity-matched  distributed  photodetector  (VMDP)  with  both  impedance  and  velocity 
matching.  The  device  exhibits  a  very  low  dark  current  and  a  high  external  quantum  efficiency  of 
0.64  AAV.  Using  the  balanced  VMDP,  the  relative  laser  intensity  noise  (RIN)  and  the  erbium- 
doped  fiber  amplifier  (EDFA)-added  noise  have  been  suppressed  by  17-dB,  and  the  RF  signal 
has  been  enhanced  by  6-dB.  This  was  the  first  integration  of  balanced  detectors  with  high-power, 
high-speed  distributed  photodetectors  for  high  performance  RF  photonic  links. 
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Figure  Captions; 


Fig.  1  (a)  Schematic  structure  of  the  distributed  balanced  photodetector.  The  inset  shows  the  active 

region  with  an  MSM  photodiode,  (b)  The  photodiodes  connected  to  the  optical 
waveguide  in  balanced  detection  mode. 

Fig.  2  (a)  Balanced  detection  setup.  The  complimentary  input  signals  are  produced  by  the  X-coupled 
MZ  modulator,  (b)  Peak  power  of  AC  signal  in  balanced  mode  versus  the  delay.  For  RF  signal 
at  8  GHz,  the  period  was  measured  to  be  38.75  mm,  which  is  very  close  to  the  theoretical 
value.  The  power  difference  between  0°  and  1 80°  phase  difference  is  more  than  44  dB. 

Fig.  3  The  AC  signal  in  unbalanced  and  balanced  detection  modes  are  plotted.  The  balanced  output  is 
almost  6-dB  higher  in  amplitude  than  the  unbalanced  signal.  Noise  suppression  of  greater 
than  17-dB  is  achieved. 


11 


Figure  1(a) 


Distributed  Balanced  Photodetectors  for  High  Performance 

RF  Photonic  Links 


TuD2 


M.  S.  Islam,  T.  Chau,  A.  Nespola*,  A.  R.  Rollinger,  S.  Mathai,  W.  R.  Deal, 

T.  Itoh  and  M.  C.  Wu 

UCLA,  Electrical  Engineering  Department,  66-147D  Engineering  IV,  405  Hilgard  Ave.  Los  Angeles,  CA  90095-1594 
^  Tel;  (310)  825-6859,  Fax:  (310)  825-6954,  email;  wu@ee.ucla.edu 

Politecnico  di  Torino,  Dipartimento  di  Elettronica,  Corso  Duca  degli  Abruzzi  24, 1-100129  Torino.  Italy 


Balanced  photodetectors  are  of  great  interest  to  analog  fiber  optic  links  because  they  can  suppress  laser 
relative  intensity  noise  (RIN)  and  amplified  spontaneous  emission  noise  (ASE)  from  erbium-doped  fiber  amplifiers 
(EDFA)  [1].  Because  balanced  photodetectors  can  achieve  shot  noise-limited  link  performance,  the  noise  figure  and 
spurious-free  dynamic  range  of  the  link  continue  to  be  improved  by  increasing  the  power  of  the  optical  carrier. 
Therefore,  balanced  photodetectors  with  broad  bandwidth  and  high  saturation  photocurrents  are  particularly 
important  for  analog  fiber  optic  link  applications.  Previously,  we  have  reported  a  velocity-matched  distributed 
photodetector  (VMDP)  with  a  peak  saturation  photocurrent  of  56  mA  and  a  3-dB  bandwidth  of  49  GHz  [2]. 
Compared  with  other  photodetector  structures,  the  VMDP  is  more  suitable  for  implementing  the  balanced 
photodetection  since  it  has  separate  optical  and  microwave  waveguides.  Here,  we  demonstrate  a  novel  distributed 
balanced  photodetector  that  can  achieve  high  saturation  photocurrent  and  large  bandwidth  simultaneously.  A  3-dB 
frequency  of  16  GHz  and  a  responsivity  of  0.64  AAV  have  been  achieved. 


Figure  1.  Schematic  structure  of  the  distributed  balanced  photodetector.  The  inset  shows  the  active  region  with 
an  MSM  photodiode.  Figine  on  the  right  shows  the  photodiodes  connected  to  the  optical  waveguide  in  balanced 
detection  mode. 


Figure  1  shows  the  schematic  of  the  balanced  distributed  photodetector.  It  consists  of  two  input  optical 
waveguides,  two  arrays  of  high-speed  metal-semiconductor-metal  (MSM)  photodiodes  distributed  along  the  optical 
waveguides,  and  a  50Q  coplanar  waveguide  (CPW)  output  transmission  line.  Beam  propagation  simulation  results 
indicate  no  crosstalk  between  the  parallel  optical  waveguides,  and  the  absorption  per  photodiode  and  the  optical 
coupling  loss  between  photodiodes  are  8.8%  and  3.2%,  respectively.  The  diodes  are  23  pm  long  and  5  pm  wide.  The 
MSM  fingers  are  patterned  by  optical  lithography.  The  central  conductor  of  the  CPW  has  a  width  of  55  pm  and  the 
separation  between  the  central  conductor  and  the  ground  conductors  is  85  pm.  The  finger  overlap  of  10.5  pm  is 
designed  to  provide  the  required  capacitance  loading  for  velocity  matching.  The  matching  ensures  broadband 
operation.  Unlike  previously  reported  slow-wave  CPW  [3]  that  ignored  the  resistance  of  the  MSM  photodiode 
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out-of-phase  RF  signals  for  balanced  VMDP.  The  output  from  the  balanced  VMDP  was  measured  using  a  spectrum 
analyzer. 


DFB  Laser 


EDFA  Band  Pass 
Filter 


X-Coupled 
Mach-Zehnder 


RF  input 


Spectrum 

Analyzer 


RF 

Amplifier 


Variable 

Balanced  Detector  Delay  Line 
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Figure  3.  Balanced  detection  setup.  The  complimentary  input  signals  are  produced  by  the  X-coupled  MZ 
modulator. 


The  device  was  biased  in  balanced  mode  by  biasing  one  of  the  ground  electrodes  of  the  CPW  at  8 Volts. 
We  used  a  custom  made  probe  with  integrated  DC-blocking  capacitor  on  the  ground  that  connects  to  the  CPW 
ground  biased  at  higher  voltage.  To  verify  the  balanced  detection,  we  modulated  the  optical  input  at  8GH2  and  tuned 
the  fiber  optic  delay  line  from  0  to  5mm  delay.  Figure  4  shows  the  RF  signal  verses  the  delay.  When  the  RF  signals 
have  0®  phase  difference,  the  AC  output  is  cancelled  in  balanced  mode.  At  180®  phase  difference,  the  RF  signal  is 
amplified,  the  extinction  ration  is  more  than  44dB. 


period  was  measured  to  be  38.75  mm,  which  is  very  close  to  the  theoretical  value.  The  power  difference 
between  0  and  180®  phase  difference  is  more  than  44-‘dB.  (b)  The  same  plot  in  linear  scale. 
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Figure  5  shows  the  RF  spectrum  of  the  output  from  the  balanced  VMDP  in  the  non-balanced  (  only  one 
waveguide  is  illuminated  )  and  balanced  mode.  Suppression  of  the  noise  floor  by  12-dB  has  been  observed  in  the 
balanced  mode  over  a  wide  frequency  range  from  6  to  15  GHz.  The  signal  is  also  enhanced  by  6-dB.  Typical 
balance  between  the  VMDP  was  within  2%  of  the  total  photocurrent  reading,  which  was  limited  to  2  mA  in  our 
preliminary  experiments.  For  bandwidths  lower  than  6  GHz,  the  fiber  length  for  one  RF  period  is  longer  than  5  mm, 
the  maximum  delay  length  in  our  current  setup.  We  are  currently  working  on  a  new  setup  with  two  fiber  delay  lines 
to  achieve  exactly  equal  optical  path  lengths  for  broad  band  operation  without  any  tuning  of  the  delay  line.  This  will 
allow  us  to  measure  the  AC  responsivity  of  the  device  in  balanced  mode. 
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Figure  5.  The  AC  signal  in  non-balanced  (a)  and  balanced  (b)  detection  modes  are  plotted  separately.  The 
balanced  output  is  almost  6-dB  higher  in  amplitude  than  the  non-balanced  signal.  Preliminary  results  show  a 
noise  suppression  level  greater  than  12-dB. 


In  conclusion,  a  balanced  velocity-matched  distributed  photodetector  (VMDP)  with  both  impedance  and 
velocity  matching  has  been  designed  and  successfully  fabricated.  Preliminary  experimental  results  shows  that  the 
relative  intensity  noise  has  been  suppressed  by  12-dB,  and  the  signal  has  been  enhanced  by  6-dB.  This  was  the  first 
integration  of  high  power  and  high  speed  balanced  detectors  with  a  slow  microwave  transmission  line  to  be  used  in 
high  performance  RF  photonic  links. 
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